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IMPORTANT NOTE 


This volume was written by and for engineers and scientists who are 
concerned with the analysis and synthesis of piloted aircraft flight 
control systems. The Bureau of Aeronautics undertook the sponsorship of 
this project when it became apparent that many significant advances were 
being made in this extremely technical field and that the presentation 
and dissemination of information concerning such advances would be of 
benefit to the Services, to the airframe companies, and to the individ- 
uals concerned. 

A contract for collecting, codifying, and presenting this scattered 
material was awarded to Northrop Aircraft, Inc., and the present basic 
volume represents the results of these efforts. 

The need for such a volume as this is obvious to those working in 
the field. It is equally apparent that the rapid changes and refine- 
ments in the techniques used make it essential that new material be 
added as it becomes available. The best way of maintaining and improving 
the usefulness of this volume is therefore by frequent revisions to keep 
it as complete and as up~to-date as possible. 

For these reasons, the Bureau of Aeronautics solicits suggestions 
for revisions and additions from those who make use of the volume. In 
some cases, these suggestions might be simply that the wording of a 
paragraph be changed for clarification; in other cases, whole sections 


outlining new techniques might be submitted. 


Each suggestion will be acknowledged and will receive careful study. 
For those which are approved, revision pages will be prepared and dis- 
tributed. Each of these will contain notations as necessary to give full 
credit to the person and organization responsible. 

Thies cooperation on the part of the readers of this volume is vital. 
Suggestions forwarded to the Chief, Bureau of Aeronautics (Attention 
AE~612), Washington 25, D. C., will be most welcome. 


Le M. Chattler 

Head, Actuating & Flight Controls Systems Section 
Airborne Equipment Division 

Bureau of Aeronautics 


PREFACE 


This volume, The Hydraulic System, has been written under BuAer 
Contract NOas 51-514(c) to present to those concerned with the problems 
of designing integrated aircraft control systems certain basic informa- 
tion regarding hydraulic systems used to position aircraft control sur- 
faces. 

The purpose of the volume is to develop mathematical models of 
typical aircraft hydraulic control systems. The analyses are used to 
show the effects of various system parameters on stability where aero- 
elastic behavior of the airframe is not a prominent factor. Portions 
of the analyses are used to outline possible flutter systems including 
the aeroelastic airplane, the hydraulic actuator, and other servomech- 
anisms of flight controller loops. 

It is mathematically expedient, though somewhat artificial, to 
treat the flutter problem as a separate entity. This separation corre- 
sponds to current practice in the study of the control of aircraft. The 
material of this volume is therefore presented in two parts. 

Part I, composed of eight chapters, is primarily concerned with 
derivation of the controlling equations of typical actuators and with 
stability of typical systems embracing an actuator and a simplified 
control surface. The problem of instability is given heavy emphasis, 
and several methods for overcoming this difficulty are presented. 
Details of design have been considered only where the reader may be 
aided in a physical understanding of the behavior of typical systems. 
The hydraulic actuator configurations which are discussed have been 


chosen as most nearly representing current design practice, 
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Part II of the volume is composed of two chapters. The first deals 


, oo 


with the conceptual modification of the classical flutter system to 
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account for the effect on flutter of a complex control system embracing 


, 


several servomechanisms. A second chapter presents certain techniques 
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for effecting a solution of the complete servo-flutter system. 
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PART I 
HYDRAULIC ACTUATING SYSTEMS 
~ CHAPTER I 
INTRODUCTION 


Part I of this volume is concerned specifically with the problems and 
tie: Seneral design criteria associated with hydraulically powered aircraft 
surface controls, The early sections deal with general considerations 
relating to this subject, and the later chapters discuss more specific 
applications of particular systems. It should be emphasized that through- 
out the volume the basic analytical approach rests upon the mathematical 
models of the transfer functions of interest. The effects of parameter 
variations upon these transfer functions and upon over-all performance are 
éteciawed in the final chapters. 

A pyamary consideration in the design of any type of servomechanisn, 
or of control systems tn general, is the ease with which the flow of power 
is controlled. Among other things which a servo does, it operates as a 


power amplifier, and therefore the facility with which the power may be 


metered into some dynamical system determines to a large degree the choice 
of various components comprising the system. Nearly all the difficulties 


encountered in the design of any device arise from a misapplication of the 


most fundamental engineering concepts. In the particular case of a hydraulic 


control device, certain design shortcomings may be attributed to improper 
applications of such basic concepts as fluid flow and system dynamics, and 
to other simpler considerations, mechanical in nature, 

With these introductory thoughts in mind, the more general consider- 


ations entailed in hydraulics analysis will now be presented. 
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CHAPTER Il 
GENERAL CONSIDERATIONS 
SECTION 1 = THE NATURE OF THE HYDRAULIC MEDIUM 


The hydraulic medium is considered to be a (relatively) incompressible 
fluid. In addition it is considered to be completely continuous and homo= 
geneous. This implies that, in the hydraulic system, there are no disconti- 
nuities, such as cavitation, free air, etc. The hydraulic medium is thus 
capable of completely filling any cavity or chamber. This hydraulic mediun, 
which will be referred to hereafter as a "liquid," is capable of withstanding 
compressive loading only. The column of liquid comprising the hydraulic system 
can be considered as a fluid piston. Therefore a displacement of a volume of 
liquid anywhere in the hydraulic system will result in an equivalent displace- 
ment elsewhere in the system. 

Consider a system made up of two movable pistons connected by means of 


external plumbing. 


Working Fluid 
under Pressure 
Ve 





A 
_ Output Load Kay or 
= Aa 2 4 Ap: 


Figure II-l,. Simple Hydraulic System 


The chambers formed by the two pistons and the bores in which they operate 
and also the interconnecting plumbing are completely filled with liquid. 


Disregarding the transient state, if one piston is displaced, ,the, entire 
TT..1 
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column of liquid undergoes an identical volumetric displacement. The linear 
displacement of the second piston is equal to that of the first piston modified 
by a proportionality factor which is the ratio of the areas of the two pistons 
exposed to the hydraulic fluid, Ap JAs « If the output piston operates against 
a load, the liquid contained in the system is subjected to a compressive loading. 
The intensity of this loading is the magnitude of the output load divided by 

the area of the output member as seen by the fluid under compression. The 
intensity of the compressive stress so generated is referred to as "pressure," 
which is generally expressed in pounds per square inch. If the velocities 
encountered in a system are not high, this compressive stress is uniform through- 
out the system. 

The product of the volumetric flow rate (corresponding to the output piston 
velocity) and the generated pressure yields a time rate of work performed, or 
power output. In the hydraulic system described above, if the load on the out- 
put member is such that a pressure of 1714 pounds per square inch is induced 
throughout the system, and if the input member is moved so that the fluid is 
displaced at the rate of one gallon per minute, the system is transmitting one 
horsepower. | 

SECTION 2 = GENERATION OF HYDRAULIC POWER 

Since the main function of the hydraulic control system is to meter 
hydraulic power into a dynamical system, a brief discussion of how this power 
is generated is pertinent. The source of the hydraulic power is a fluid 
punp of some type, usually driven by the main propulsion unit of the aircraft, 

In general, two types of pumps are available: these are classified as 


hydrokinetic and hydrostatic, 
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The hydrokinetic type pump involves the handling of relatively large vol- 
umes of fluid at fairly low pressures. An example of this type of hydraulic 
power generator is the well-knowm centrifugal punp. In a device of this sort, 
large volumes of fluid are received in a central chamber. This fluid is accel- 
erated in a radial direction und attains high velocity. The momentum thus 
imparted to the fluid is delivered to the high pressure outlet of the pump as 
the fluid is decelerated in the outlet chamber, Since this type of fluid 
power generator is not normally used in aircraft hydraulic systems, only pass= 
ing mention of the device is made here, 

Because minimization of the weight und volume of components used in aire 
craft systems is of prine inportance, the hyurostatic, or positive displacement, 
type pump is almost universally used, ‘he simplest type of pump (see Figure 


II-2) is a simple reciprocating piston. 


Return 


Piston Reciprocated 
by Prime Mover 





Its 
Suction ‘J nr 


CS Grog eng os @& 


Sump 
Figure II-2, Simple Reciprocating Piston 
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The piston is connected to suction and discharge ports by means of appropriate 
valving such that on the suction part of the stroke the chamber formed by the 
piston and the bore in which it operates is completely filled with fluid taken 
from the suction line. During the reverse, or discharge, stroke this volume of 
fluid is delivered to the discharge line at a somewhat higher pressure. If the 
piston is driven continuously in an approximately sinusoidal manner by means 
of a crank mechanism, a simple positive displacement pump is the result. 

A device of this type is suitable for use with a compressible fluid; an 
example is its use as an air compressor. However, this kind of pump cannot 
be used as a hydraulic power generator Since tHe discharge rate throughout the 
cycle is not constant, or even approxfmately so. In order to provide a satis- 
factory hydraulic power generator, a mult#plicity of such pistons must be pro- 
vided. In addition, the angular phasing of the pistons should be such that 
each piston performs its suction and discharge functions in a fixed phase 
relutionship to the piston which follows it. To minimize pump output ripple, 
an uneven number of pistons is employed. 

Valve Plate (Stationary) 


H Cylinder Barrel 
Piston (One of, say, nine ) 


D s ANY; — Radial Bearing 
SS Thrust Bearing 
{\ 


Pump Housing fi 








WN Shaft Driven by 
S Prime Mover 

Universal Joints 
socket Ring 


Figure II-3(a). Positive Displacement. Hydraulic Pwnp 
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Stationary 





Note: Socket ring thrust and radial bearing not shown. 


Figure II-3(b). Positive Displacement Iiyaraulic Pump 


A typical configuration (see Figure II-3) consists of, say, seven or nine 
pistons (A) operating in the bores of a rotatable cylinder barrel (B). The 
pistons are connected to a rotatable socket ring (C) by means of short connect- 
ing rods (D) and ball and socket joints (E). The socket ring and the cylinder 
barrel are synchronized by means of a smal] comnecting shaft (F) and universal 
joints (G) as required. Consider that the socket ring is rotated by some 
prime mover such as the main propulsion plant. Therefore, the cylinder barrel 
and the pistons also rotate about a longitudinal axis. If the socket ring is 
displaced angularly about a lransverse axis normal to the axis of rotation, 
continuous reciprocation of the pistons relative to the cylinder barrel results 
as the socket ring is driven by the prime mover. Because the bores in which 
the pistons operate are arranged parallel to the axis of rotation of the cylin- 
der barrel and are located at the same radial distance from the axis of rotation, 
and since the angular spacing of the bores is equal in each case, the desired 
fixed phase relationship between the individual pistons is achieved. The 


remaining necessary feature is the valving means which alternately connects 
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the pistons, as they reciprocate, to the suction and discharge lines of the 
hydraulic system. The cylinder barrel is provided with a flat surface normal 
to the axis of its rotation. The chambers of variable volume, formed by the 
piston heads and the parallel bores in the cylinder barrel, communicate with 
the flat surface of the base of the cylinder barrel by means of drilled holes 
(H) of reduced diameter. These holes, of course, move through a circular path 
as the cylinder barrel rotates in synchronism with te socket ring. Rotation 
of the cylinder barrel through 180° moves some reference piston into the cyl- 
inder barrel, thereb, reducing the volume of fluid within the cylinder. Dur- 
ing he remaining half cycle, the piston travels in the reverse direction, 
drawing fluid into the cylinder. All that remains to be provided is means 

to isolate these two regions of operation from each other and to connect them 
respectively to the discharge and suction lines of the hydraulic System. 

The flat base of the cylinder barrel rotates against a flat surface which 
forms a valve plate (J) fixed to the pump housing (K). The valve plate is 
equipped with two kidney shaped ports (L). These two ports are so arranged 
that one of them communicates through the holes (H) with the chambers in which . 
the pistons are advancing and for which the volume of the chambers is there- 
fore decreasing; and the other port gives access to the chambers in which the 
volume is increasing. 

At the two dead center positions, corresponding to the point where any | 
piston reverses direction, commutating lands separate the two ports of the 
valve plate. The width of these two commutating lands is approximately equal 
to the diameter of the holes in the flat surface of the cylinder barrel. The 
two ports of the valve plate are connected directly with the suction and dis- 


charge lines of the hydraulic system. The configuration of the cylinder bores, 
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the reduced diameter holes, and the semi-circular ports of the valve plate are 
arranged by means of careful design so that the system is in a state of hydrau- 
lic balance. Thus internal pressure has no tendency to disturb the intimate 
surface contact of the rotating cylinder barrel with the stationary valve plate 

If the angular relationship of the socket ring to the cylinder barrel 
about the transverse axis is fixed, the pump is a ccnstant displacement type. 
If this angular relationship can be changed during the normal operation of the 
pump, it becomes a variable displacement type. Both are widely used in air- 
craft installations. | 

The type of equipment just described when driven by a prime mover is a 
hydraulic pump. The same basic unit when supplied by a source of hydraulic 
fluid under pressure can be used as a motor. The comtination of a variable 
displacement pump driving a constant displacement motor constitutes the basic 
configuration of a positive displacement hydraulic transmission. This device 


is widely used and is to be considered in some detail in this chapter. 
SECTION 3 = FUNCTION OF THE HYDRAULIC SERVOMECHANISM 


The function of the hydraulic servomechanism is to control the flow of 
hydraulic power. This power is in the form of a volumetric flow rate of a 
(relatively) incompressible fluid operating, generally, against an appreciable 
load-induced pressure. This volumetric flow is used to actuate some type of 
servo motor. The servo motor or actuator is used to control the position, 
velocity, or acceleration of some dynamical system which is made up of lumped 
or distributed elemerits of inertia, damping, and elasticity. Therefore, the 
dynamical nature of the load and the attendant coupling of the dynamics of the 
load to those of the servo actuator with its associated control elements, 


present a problem of some complexity. Generally the power level of the output 
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of the hydraulic servo is many times greater than the power level of the input 
signal. Thus the hydraulic servomechanism operates as a power amplifier 

Because the unit is actuated by the existence of an error, generally 
appearing as a valve displacement, a feedback loop and an error sensing element 
are necessary parts of the system. The hydraulic servomechanism must control 
the position, velocity, and acceleration (and sometimes higher devivatives) of 
the output control ember in a stable fashion throughout the ewmtire frequency 
spectrum. Marginal instability exhibited by a condition of incipient hunting 
is sufficient reason for rejecting any specific servo design, 

The functions and performance requirements of the hydraulic servomechanism 
are of course identical to those of any other type of servomechanism. In the 
functional consideration of the hydraulic servomechanism it is necessary to 
study those particular features which make the use of hydraulic components par- 
ticularly attractive to the controls designer. Those features which are worthy 
of consideration are enumerated as follows: 

1. Size 

2. Weight. 

3. Power level 

4. Speed of response (torque to inertia ratio) 

5. Reliability 

6, Cost 

7. Ease of control 

The first four items of this list are those in which the hydraulic servo 
has an advantage over other types through a certain range of power levels. The 
data in Table II-1 indicate the important features of several commercially 


available hydraulic motors adaptable directly for use as servo actuators. Most 
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VICKERS INC. MOTOR NO. 


MF-2/,-3906- | MF~36-3908- |MF-3909- 
30BC-2 30Z=-2 30Z~2 

ears [= [oe [om 
Pressure 


Maximum R,.P.M. 
Displacement in? / rev. man 0.240 0.367 
Torque at Maximum in.-1b. 175 

Pressure 
Horsepower at Maximum 

Pressure & R.P.M. 

1b.in.sec* 1.09 x 10 4.95 x 10° 

mimeo sol 

0.015 0.038 0.058 
rad, / sec*%! 410,000 230,000 171,000 


Table II-l. Rotary Hydraulic Motor Data 


MOTOR 





CHARACTERISTICS 









10.4 


Moment of Inertia 1.02 x 10°* 





Volume of 011 Under 
Compression 


Approximate Volumetric 
_Efficiency 


Weight lbs, 





Torque to Pressure 
Ratio 


Torque to Inertia 
Ratio 


el) 


Hydraulic Compliance* 


*Based on a bulk modulus for the hydraulic fluid of 250,000 psi. 
(Compliance is reciprocal of "stiffness.") 
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of the features listed are from the manufacturers' catalog data. Some of the 
others are calculated from these data, 

For a high performance servo, the most important item on the list is the 
torque to inertia ratio. This ratio is often given as a figure of merit for a 
Servo motor and is a direct measure of the available acceleration, Of almost 
equally great importance is the low weight (and corresponding small size) of 
any of the motors listed for the given power rating. Comparison of the items 
of this list with those for an equivalent electric servo actuator clearly dem- 
onstrates the superiority of the hydraulic servo, 

The cost of hydraulic apiponeite is generally fairly high because of the 
necessary precision of the manufacturing techniques. ven so, the cost of 
hydraulic components compares favorably with that of equivalent components of 
other types. This is especially true in the field of aircraft control systems, 

In general, the reliability of hydraulic components is excellent. The 
flexibility and ease of control further enhance the over~all excellence of the 
hydraulic servo. It is believed that the ease with which hydraulic power can 
be controlled will become evident later in this chapter when the control of 


hydraulic power is more specifically considered. 
SECTION 4 - TYPES OF HYDRAULIC SERVOMOTORS 


Hydraulic servomotors can be classified under two generic headings; they 
are the linear and the rotary types. The former of these is the direct acting 
hydraulically controlled piston actuator. This type of actuator generally 
operates with the chambers at each side of the piston head pressurized. (An 
exception to this is the use of this kind of actuator in the so-called "open 


center" system.) The piston may be of the balanced area type in which the 
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areas exposed to the respective pressures at the two sides of the piston head 
are approximately equal. With such an actuator, the pressures on both sides 
of the piston head, and the resulting fluid flows, are precisely controlled by 


complementary servo elements. (See Figure II-4(a) ) 






yo 
im, 





(a) Balanced area piston type actuator. 
(Pressures ( and # at each side of piston head 
both controlled by valving means.) 


Constant = 
High Pressure 14 A 





(b) Differential area piston type actuator. 
(Only one pressure, 4 , controlled by valving ineans.) 


Figure II-4. Types of Servo Actuator 


An alternate type of linear servo actuator is the piston type in which 
the ratio of the areas of either side of the actuator piston head is approx= 
imately two to one. (See Figure II-4(b).) In such a device, the small area 
chamber is ported directly to the constwt pressure system supply. The pres- 
sure and the flow rate to the larger area chamber control the actuator per- 


formance. 
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The rotary type actuator assumes a variety of forms. Practically any type 
of positive displacement hydraulic pump can be adapted for use as a rotary 
servo actuator. These include axial and radial multiple piston actuators as 
well as gear and vane type actuators. 

In general, any hydraulic servo actuator, linear or rotary, can be con- 
trolled either by the use of conventional valving or by the use of a servo con= 


trolled vuriable displacement pump. 
SECTION 5 ~ POSITIVE DISPLACEMENT HYDRAULIC TRANSMISSION 


The positive displacement hydraulic transmission consists basically of a 
variable displacement pump and a fixed displace;nent motor hydraulically inter- 
connected in such a way that the volumetric output of the pump drives the 
hydraulic motor. The variable displacement pump is driven at a constant speed 
by some type of prime mover. The angular velocity of the motor is proportional 
to the volumetric flow rate to it; the latter is proportional to the angle of 
tilt of the socket ring in the pump if the device is of the same tye as that 
described earlier in this chapter. Therefore, servo control of the output 
sember, that is, of the fixed displacement motor, is achieved by control of the 
tilt angle of the variable displacement pump. Generally the variable displace- 
ment pump is referred to as the "A-end," and the fixed displacement motor is 
referred to as the "B-end." The simple system is as shown: 


A-end Tilt Control (X) 





©, 
Wy We, 
VVNUVaPnrpdagyy © 
oe 
aJ To 
ffosssasse-nfled Load 
Drive Motor A~end (Pump) B-end (Motor) 


Figure II-5., Hydraulic Transmission 
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The basic equation involved can be expressed as: 
(II-1) Qn" y* & +O. 


where 

Qp> is the total volumetric flow rate from the pump 

, is the flow rate consumed by motor rotation 

Q, is the leakage flow rate 

Q. is the compressibility "flow rate" 
It is possible to define each of the quantities above in terms of known system 
characteristics .* 

> varies from zero to a maximum in stepless increments as controlled 


by the A-end stroke (X). X varies from zero to + 1.0. Therefore 
(11-2)  Qp = Xow, (in? / sec.) 


where 
L is the volumetric displacement of the A-end at full stroke (X“ 1.19) per 
unit of angular displucement (in? / radian). 
4); is the (constant) angular velocity at which the A-end is driven. 


Also: 
(II-3) Q,= 4A (ind / sec.) 


where: 
L is the system leakage factor of the complete transmission, (This 
is derivable from catalog curves of volumetric efficiency versus pres- 
sure drop across the motor. This is very nearly a linear function.) 
*The following derivation of (II-6) is essentially that given in 


Brown, G. S., and Campbell, D. P., Principles of Servomechanisms, John 
Wiley & Sons, New York, 1948. 
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Ff} is the pressure drop (psi) across the motor. (This is entirely a load 


induced pressure. ) 


Further: 
-% GR (3 
(II-4) O& N dt (in? / sec.) 
where: 


% (is the total volume of liquid under compression in the active system. 
(This is a const a:t and does not vary with motor position.) 


/V is the bulk modulus of the liquid. (See Chapter III, Section 2.) 


Finally: 
(II-5) Ow = Lay Gy (in? / sec.) 
where: 


d,, 18 the volumetric displacement of the B-end per unit angular dis- 
placement. 

G}, is the angular velocity (variable) of the B-end. 

Combining (II-?), (II-3), (II-4), and (TI-5): 


% d 
(11-6)  Xdpa = dah + LB + —6 oa 


At this point, it is convenient to introduce the quantity Z, , which is 
referred to as "load impedance." 


By definition, 
yi 
Z,=— 
LZ Qn 
therefore, from (II-5): 
RR = 2, a, 


Substituting this value of / into Equation (II-6): 
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(11-7) Xdyaly = Ayah + 2, dy oh (Z + és] 


The Laplace transform is introduced at this point to simplify presentation of 
the equations, 


The open loop transfer function of the hydraulic transmission is then: 


Ga) al. uly / 
x (3) ain 


nn ra 
d AC ¢ % s) 
For the relationship between the angular output position (@,) of the B-end 
and the control input (X) , there results: 


(1-9) 3G) = — 
1+ (l 4 ¥e 5) 
where 
K = Met 
Lp, 


Several simple cases of specific types of loading are considered: 
Case (a) 
The case wherean inertia element is rigidly coupled to the B-end output 


shaft, 





The load torque, 


T = Ss* 
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Since 
Pp 
z= 
Z Om 
P 
Z, = f 
‘ sq 
Also, the torque, 
T= ch 


where ¢ (in?) is the motor torque per unit pressure (inch pounds per psi 
active system pressure differential). This value is often given as part of 
the manufacturer's model designation. 


Therefore, in this case, 
2 
z= JSS OQ, 
cd,s&@ 
or 


SS 
cd, 





Z,(5) = 


and the system (open loop) transfer function becomes 


lm (cat Jls + 1% ot 
N 
or 
2) K ; 
_ 9 Or age pe = agen nga 
(II-10a) (5) s (Le eee @ 1} 
Case (b) 


The case where an inertia element is coupled to the B-end output shaft 


by means of an elastic member and with output damping, 
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As before, 
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The load torque, 


T= k(9,-9,) = vs*6, + Bs 


L 


fron which 
s k & 
. J/s?+ Bs +k 


Then the torcue, 


k &o 


2 
T=(Us + Bs): Bs tk 


Since 
Zz, - Jory ue 
cd,,$ Op 


k(/s + B) 


and the systom transfer function becomes 


d(/s*» Bs tk 
aay Sele) - BO arn 


cd, (/s*+ Bs + k)+ klLl/s+B)+ lust» Bs) 


ae. @ (Js? + Bs + k) 
~ $ 


JUk%\ + f Jk L Ea / BAL J 
(vB) *(8* ce * ca? (* * ca 


II~17 








Section 6 


Case (c) 
This system is the same as (b) but with zero damping (B=O) , 


Lo + / 


(II-12) <2(s) = — @ 2a) 


k ca,N Cd, 


In a similar manner, load dynamics of any degree of complexity can be 
included in the over-all system transfer function. Study of the transfer func- 
tion above yields the design requirements of the servo device which controls 


the stroke (X) of the A-end. 
SECTION 6 = SEPARATION OF LOAD DYNAMICS 


In the preceding section, the concept of load impedance was introduced. 
This was done in order to establish a means whereby the load dynamics of a 
system may be separated from the dynamics of the unloaded servo device. Al- 
though it is not necessary to separate the load dynamics in this manner, doing 
so sometimes yields a better physical insight. The techniques involved may 
appeal to certain engineers, and this section is presented for that reason. 

In general the dynamics of the unloaded servo element are fairly easy to 
define. If no other means are available, a series of simple bench tests will 
suffice, in many cases, to establish the important parameters applicable to 
the unloaded servo element. If the dynamics of the load can be handled as a 
separate entity, the over-all problem of determining system behavior becomes 


somewhat simplified. 
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In order to clarify this idea, the biock diagram below is considered. 


fe re eK OO eee — 





] — KsGs——, 
] 
e, ,D=1 UNLOADED | i ,D>? | ee 
i. SERVO 
le 


Figure II-6., Load Block Diagram 

In the diagram the transfer function of the unlosded servo system is des- 
ignated by A,G, . For a condition of no load the system output (a) is identi- 
cal with @,, since the contribution of the load dynamics is zero. For the 
system éperatine under an appreciable load, the unloaded system output (a) 
differs from the actual output (9) by the amount (b) supplied by the load 
dynamics. It is well tc keep in mind that the diagram does not necessarily 
represent any actual physical array of hardware. Parts of the diagram may be 
purely fictitious, and such an arrangement is used here only to provide a means 
of understanding the over-all problem. 


if the diagram is accepted as such, the following relation:hips are appar- 


ent. 
(I) G@-=a-b 

b is subtracted from @ by the fictitious differential device D~2, 
(II) E=-=O-& 


E can be considered the instantaneous error and is formed by comparing 
system output (@,) to the input ( OG; ) at the differential device D-l. 
(III) a = (KG )E 


KgG,; is the open loo, transfer function of the unloaded system. 
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(IV) b= (1, G)& 
Ay & is the "transfer function" of the load dynamics. 
(V) & = (K; 6, )E 
A, G, is the actual open loop systom transfer function under load. 


Combining (I), (III), and (IV), 
& = (Ky G5 )E - (KG, )@ 


and 


S _ C=) 
E Ks Gs /+ KG 
io 


Since “F is by definition Kk; Gr, 


(II-13) K>G, = Kz ares) 

This last expression relates the over-all system transfer function to the 
transfer function of the unloaded system and that of the load dynamics as a 
separate entity. This may be compared with the general expression for the over= 
all system transfer function of the positive displacement hydraulic transmission, 


Equation (II-9), which is repeated below as: 





(II-14) SP (s)= bal f = 
J; + Z, ( + és) 
In this expression the unloaded system transfer function is that part of 


the right hand side not enclosed in the brackets. Therefore, in this case: 
K 

$ 

This indicates that the unloaded positive displacement hydraulic transmission 


KG, * 


acts as a perfect integrator, as indeed it does, That is, for a fixed position 


of A-end tilt, a proportional B-end output velocity results, 
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Further comparison of Equations (II-13) and (II-14) shows that the transfer 


function of the load dynamics is: 


(II-15) KG, = 2, (0 + % } 
N 

The performance of a well-designed servo system employing a positive 
displacement hydraulic transmission will be such that in the range of fre- 
quencies covered by the. input signal, the device will operate very nearly like 
a perfect integrator. In the higher frequency range where the effects of the 
load dynamics become more pronounced, there must be no evidence of incipient 
instability. In addition, of course, the device must be capable of meeting 
whatever performance specifications are applicable, 

These last comments apply, in a general sense, to the valve controlled 
actuator servo device. This type of servo element is discussed in the follow- 
ing section. An attempt is made to show the similarity of the valve controlled 
servo, which is basically a non-linear device, to the positive displacement 


type, which is essentially a linear device, 
SECTION 7 = VALVE CONTROLLED SERVO ELEMENT 


The valve controlled hydraulic actuator is probably more widely used as 
a servo element in aircraft control surface actuation than any other device, 
A typical configuration employing a closed center four-way control valve with 


unity feedback is shown in Figure II-7,. 
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Figure II-7, Valve Controlled Hydraulic Actuator 
In the figure, & is system supply pressure, and & is system return pres~ 
sure. For the purpose of the present discussion, Pp is assumed to be zero and 
fa constant. 
Control of the actuator is accomplished by simult‘neously controlling the 
rolunetric flow to and from the two chambers adjacent to the two sides of the 
}. ston nead, The flow to the cylinder through an inlet orifice is given by 


the basic relationship: 


(1-16) Q*Ca/2g” 


(This equation can be readily derived from Bernoulli's equation. See, for 
example, Vennard, J. K., Elementary Fluid Mechanics, John Wiley & Sons, New 
York, 1940.) 
where 

C is the orifice coefficient 

@ is the orifice area 

uy is the specific weight of the fluid 

Ap is'the pressure drop across the orifice 


g is the acceleration due to gravity 
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Assuming an essentially closed-center valve configuration, i.e., negligible 
leakage flow past the valve lands at a centered valve position, the flow 
relationship above also expresses the volumetric time rate of cylinder dis- 
placement. Further, the pressure drop,4p , across the orifice becomes 

a simple function of the system pressure and the load induced pressure 
RP (P= P-f) . Since the inlet orifice area is a direct function of 
the instantaneous system error £ (or valve displacement from neutral), it 
follows that the effective flow tending to move the cylinder is a function 


of the three variables just named, i.e., 
Q= FL, 2, R) 


If it is assumed that the total flow Q consists of a steady flow.Qr and 
a perturbation, g , about this steady valuasQ= Ong 3 then the perturhed 
flow can be approximated as the first-order Taylor's series expansion of Gg 
as a function of £,/, and: 


| dQ JQ dQ 
a a" yeh * Spat TBA 


but since the system supply pressure has been assumed constant: 


where, by definition, 





(II-19) C= a : S = — S 


The perturbation flow rate Gg may also be thought of as the sum of two 
components: 
(1) An incompressible component, 9, , which causes motion of the cyl- 


inder relative to the piston. 
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(2) A compressible component, g, » due to the compliance of the fluid 
within the cylinder; this component causes no actuator motion. 


Therefore 
g = Z 4 J. (leakage flow is assumed negligible) 


where 
G = Asx, 


and is that part of the flow which results in actuator motion; 


also 
4 


ee 
G. ~ WN OF 
and is that part of the flow through the valve orifice which does not result 
in actuator motion. 


Therefore , 


Gy = AsX, + TSR = GE- GA 


or 


A 
(II-20) CoE = AsXo (cp ? Xs)a 


where 
A is the area of the actuator piston 
¥” is the effective volume under compression 
Xo is the actuator output position 
re 


is the system instantaneous error 

2, is the load induced pressure | 

N is the bulk modulus of the fluid 

All the variables above are perturbation quantities. The detailed treat- 


ment of this part of the description is presented elsewhere in this volume. 
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The foregoing makes no attempt to justify completely the validity of the state. 
ments made. The reader is referred to Chapter III of this volume which treats 
this subject in greater detail, Specifically, the reader will need to under- 
stand what is meant by "the effective volume, rr." 
The concept of load impedance is again applied to compare this type of 
system with the positive displacement hydraulic transmission. 
By definition, 
ang 
O 
Since 
g = Asx, ; A, = 2, AsX, 


Substituting this value of % into Equation (II-20) gives: 
y’ 
(II-21) CoE = ASX_ + Z,AsXo[-s - >) 


Writing this expression in the form of the open loop transfér function 


results in: 


This may now be compared with the equivalent expression for the positive 
displacement hydraulic tranamissions: 


(11-23) 2(s) = 4 ae] 


/+z were 


The two expressions are: identical in form, with the load dynamics iso~- 
lated in each case. However, some very important differences exist. In 
(II~23), all the system parameters as given are constant and are fairly easy 


to evaluate, 
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These constants are: 


dy, 

Gp uty system gain term 
Ly 
% volume of oi] under compression 


l leakage factor (which gives rise to system damping) 


However, in the case of the valve controlled actuator device, none of the 
equivalent terms is a constant except over a very small operating range. Over 
any appreciable operating range, the problem is to assign reasonably accurate 
average values to Ce , Cp , and & Z Actually in practice it is fairly easy 
to establish a realistic value for C;. For conservative design the walue of 
¥ ‘tor the piston at the midpoint of the stroke may be used since this is the 
point of maximum hydraulic compliance, 

The most elusive parameter is the effective damping factor, Co . No simple, 
direct method of evaluating this quantity, short of actually testing an exist- 
ing device, is known to be completely satisfactory at this writing. 

The reason for presenting the positive displacement hydraulic transmission 
and the valve controlled actuator in the manner in which this was done was to 
bring out the rather remarkable similarity of the two devices. Actual fre- 
quency response tests carried out on examples of the two types operating against 
similar loading and at the same power level yield very similar results. This 
is true in spite of the fact that in one case the system is linear throughout, 


and in the other case the system is very definitely non-linear, 
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CHAPTER 
ANALYSIS OF THE GENERALIZED HYDRAULIC SERVO ACTUATOR 
| SECTION 1 - INTRODUCTION 


In this chapter, an analysis of the generalised hydraulic servo actuator 
shown in Figure III-1 is presented. 
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Figure III-l1. The Generalized Hydraulic Servo Actuator 


Since the input, or control member, is subject to only a portion of a 
static load on the control surface, the mechanism shown in the figure is a 


power boost actuator; that is, the hydraulic system provides only a portion of 
the energy necessary to move the output, or controlled member, against the 


imposed load. The rest of the necessary energy must be supplied by another 


source at the control end of the actuator. In an aircraft application, this 
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additional mechanism is either the human pilot or a flight controller servo motor, 


The energy requirements for actuating aircraft control surfaces have increased 


steadily through the years, requiring the hydraulic system to provide an increas=- 
ing proportion of the energy supplied. The limiting type — that in which the 
control member is not subject to any appreciable static load — is called the 
fully=powered actuator, and is of great importance at the present time. There 
remains, however, considerable interest in the power boost type. The mechanism 
of Figure III-l1 has been selected so that the results of its analysis may be 
specialized to describe the behavior of a large percentage of the hydraulic servo 
actuators in current use. 

The slide valve in Figure III-l is of the conventional four-port type and 
is mounted integrally with the hydraulic cylinder. A more general actuator in 
which the valve housing is mounted on a linkage rather than directly attached 
to the cylinder would require additional flexible tubing, and the follow-up 
linkage would be likely to have considerable backlash. Because of these dis- 
advantages, the more complicated system is of limited importance, The actuator 
of Figure III~l may be regarded as typical of practice. 

To accomplish its purpose of effectively duplicating input displacement in 
output displacement regardless of load, a positional actuator must have one or 
more feedback loops. The presence of the feedback makes an actuator a servo —- 
mechanism and therefore makes an investigation of its stability a critical part 
of the analysis. The mechanism of Figure III-l has more than one feedback loop, 
but the principal feedback is due to the integral mounting of valve and cylinder, 

The valve may in general be one of the types broadly classified as “open 
center," "closed center," and "open center —closed center"; that is, the valve 


spool in a central position may fall short of covering the operating ports 
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leading to the cylinder, it may overlap these ports, or it may come very near 
to covering the ports exactly. The range of operation of a slide valve will in 
general include a central region of open center type behavior and adjacent 
regions of substantially closed center behavior. In the region of open center 
type behavior, the flow from an operating port to the cylinder, i.e., the flow 
tending to move the cylinder, is the difference of two flows — that from the 
pressure port past a valve spool land into the operating port and that from 
the operating port past the valve spool land to the sump port. 

A hydraulic servo actuator which includes a predominantly open center type 
valve is characterized by considerable position error when operating under load. 
On the other hand, with an open center type valve, the fluid on each side of 
the cylinder is always under compression, minimizing the entrainment of air and 
the resultant reduction of oil rigidity. The "dead band" characteristic of & 
completely closed center type valve (one with overlap) is regarded as undesira- 
ble, and most present practice is a compromise, utilizing valves which are 
essentially of the closed center type, but which have a small central region 
of open center type behavior. 

The analysis to be presented will be applicable to both a region of closed 
center type valve operation and a range of open center behavior, although the 
result of its application in the latter case will be only a good approximation, 
becoming less and less accurate as the range of open center behavior is made 
more extensive. 

The mechanism under discussion is inherently nonlinear. Known mathemati- 
cal techniques do not at present provide means for solving, in other than numeri- 
cal form, the true nonlinear problem. However, a linearized treatment is suffi- 
cient basis for designing to meet most requirements, and the following analysis 
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is in terms of small increments, or "perturbations," about an operating point 


described by a fixed set of values of independent variables. 
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plan of attack is as follows: 

Derivation of flow equations of the hydraulic servo. 

Derivation, under fairly general hypotheses, of the force and linkage 
equations of the servo. 

Combination of (1) and (2) to secure a literal description of over-all 
servo behavior. 

following assumptions govern the analysis: 

Linearity exists about operating points; i.e., the relationship between 
any two variables may be regarded as proportional within a sufficiently 
small range of values of the variables, 

All masses and dampings can be replaced by lumped parameters, 

The valve is essentially symmetrical; i.e., the valve spool and valve 
housing may assume such a relative position that all areas open to 

flow are approximately equal in magnitude and geometrically similar, 
Both valve spool and valve housing are rigid. 

The valve is designed in such a way that the rate of change of momentum 
of the fluid as it passes through the valve is negligible resulting in 
the effective elimination of centering and decentering forces on the 
valve; or, equivalently, the numerical values assigned to an input cou- 
pling spring and a valve damper (see Figure III-1) are altered by the 
proper amounts to account for the valve forces, 

The fluid in the cylinder is always compressed to the extent that cave 
itation is negligible. 


The piston is double-ended. 
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8, Flows through ports connecting valve and cylinder are substantially 
equivalent, 
9. The cylinder is rigid; or, equivalently, allowance is made for cylinder 
flexibility when considering oil compressibility. 
10. There is no leakage past the piston. 


SECTION 2 ~ FLOW EQUATIONS OF THE ACTUATOR 


(a) VALVE ERROR 

Referring to Figure III-l1, it is evident that the flows from operating 
ports of the valve to the portions of the cylinder separated by the piston are 
primarily dependent on the position of the valve relative to its housing. It 
is desirable to introduce a variable to represent this relative position. In 
the following, this variable is called the "valve error." Since both the valve 
spool and valve housing are assumed rigid, no deflections of these structures 
need enter the definition of valve error. 

If X, is the displacement of an arbitrary point of the valve spool, and 
if X, is displacement of an arbitrary point of the cylinder, both measured 
from a point on a fixed structure, a valve error —& is defined as to form by 
the relationship £ ™ X,-X,. It remains only to choose a relative position of 
valve and cylinder which will be associated with a zero value of the error. A 
convenient choice for this purpose is usually the relative position which would 
be determined in test by capping the two operating ports with pressure gages, 
making supply and sump connections, and moving the valve spool relative to the 
housing until the operating port pressures match. For a symmetrical valve, a 
relative position determined in this way will be one for which all areas open 


to flow between valve and cylinder are equal in magnitude and geometrically 
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similar, and for which the common value of operating port pressures is one-half 
the sum of supply and sump pressures. 
Let the perturbations of X,, X, , andf& bea, ,x, , and E, respectively; 


then the error perturbation is given by 


(III-1) E€= x,- x, 





Supply 
Pressure 
Ph 
Valve Slide 
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Cylinder 
Piston 


Figure III-2. Flow, Pressure, and Displacement Definitions 


(b) FLOW FROM THE VALVE 

Figure (III-2) shows flows Q, and Q, , from operating ports into the cyl- 
inder. The flow Q, is of necessity a function of the supply pressure ? , the 
sump pressure R , the valve error £ , and the pressure / in the cylinder vol- 
ume on the right side of the piston. Similarly, Q, is a function of — , 


R , & , and the pressure / in the cylinder volume to the left of the piston. 
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If the range of open center type operation of the valve is not extensive, there 
will be little interaction between the hydraulic actuator and either the supply 
or the sump systems; for this reason f and f are regarded as constants in 
the following analysis. This leaves Q, and Q, as functions of two variables: 
Q, is determined by P and F ; and ah » by P andé. 

To construct a simple, practical analysis of the actuator, it is desirable 
to reduce the number of quantities on which the flows depend. In particular it 
is preferable to restrict the analysis so as to permit the use of a single vol- 
umetric flow variable Q instead of the two flows Q, and Q, » and also a sin- 
gle load induced pressure variable / instead of the two pressures / and 
The first requirement is met by restricting the mathematical model of the actu- 
ator so that the flows Q and Q, are equal in magnitude: 


(IlI-2) Q,2=-Q, = Q 


For the elimination of one of the pressure variables, the variable / is 


defined as the pressure differential across the piston: 
(III-3) p23 P-P 


A rigorous correlation on the basis of the expressions (III=2) and (III-3) 
is not possible in general, and it is a matter for inquiry to determine what 
limitations the use of these expressions places on the actuator analysis. These 
limitations will now be illustrated for the case of a completely symmetrical valve, 
i.e., a valve in which all) areas open to flow between valve and cylinder are equal 
in magnitude and geometrically similar at zero valve error, 

The functional dependence of the flow on relevant parametera is obtained by 
use of the fundamental orifice equation for substantially incompressible flow: 


(III=-4) Q = Ca ia ila 
Gy 
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where 

Q2 is the orifice area 

AP is the pressure differential across the orifice 

g is the acceleration due to gravity 

C is the orifice coefficient 

ay is a mean specific weight of the fluid 
In applying the orifice equation to a slide valve, the area is a function only 
of valve error. However, the functional dependence of this area on the error 
is rather complex. In an open center range of valve behavior, the valve slide 
forme four orifices with the ports. Flow through each of these orifices may be 
expressed, in accordance with the orifice equation given above, as the product 
of a function of valve error and the square root of a pressure difference, 

For a symmetrical valve, the four functions of valve error must be equal 
in magnitude at equal values of the error, and must be such that the two pairs 
of functions are symmetric with respect to zero valve error, The net flows Q, 
and Q, from valve to cylinder (see Figure III-2) are each the difference of 


two orifice flows, and expressions for these quantities may be written as: 


cams) 0, = WENR™ FR’ - VENP- B® 


(a6) = YHENB™ B- WEN -B 


where the asterisks indicate a fixed value, and YE) is a function of the gene 
eral type shown in Figure III-3. 
If the net flows Q, and Q, are each taken as zero in (III-5) and (III-6), 


the result is: 


(I1I-7) Pte = ere 


III-8 


aa ee We ee Oe ee WO Oe WD 


i. 


Section 2 


from which, using the definition (III-3), 


* 
(III-8) ps (Re+ B+ kh 


2 


and 


(III-9) R= 


Figure III-3. General Form of the Error Coefficient of 
Flow Through a Single Orifice 
If (III-2) is assumed valid, and the error is large enough in a positive sense 
so that Y/-F) is effectively zero, (III-7) is obtained, and (III-8) and (III-9) 
follow. Again, if the error is large enough in a negative sense so that W (E) 
ie effectively zero, (III-7), (III-8), (III-9) again apply. 

From the above remarks, it appears that if (III-2) is valid and the valve 
is symmetrical, a correlation in? , rather than in f and , is possible 
in regions of closed center type valve operation and also in that part of the 
open center range of valve error where Q= 0 . It may be shown that for inter- 
vening ranges of valve error a correlation in / alone is not possible, but 


that an analysis on the basis of (III-2) and (III-3) will yield an excellent 
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approximation when the region of open center type operation is not extensive, 
The single flow equation for the case of the symmetrical valve is obtained 
from (III-2), (I1I-5), (III-6), (III-8), and (III-~9): 


ce a a a 


For a valve having a marked degree of asymmetry of behavior, a correlation 
in? becomes dubious. This is the motivation for Assumption 3 given in the 
introduction to this chapter. 

All of the following analysis assumes the approximate validity of a corre-~ 
lation of flow behavior in terms of a single flow © and the pressure difference 


fh . An equation in perturbed quantities may then be written as: 


> 
(men) g- 226+ 2p 
L 





Figure III-4 shows typical variations of net flow Q with pressure differ- 
ential P , valve error £ being held constant; and of Q with & , /P being 


held constant. 






t P. = Constant 
! Q E = Constant + | R=0 
- Q 
| Q Ay 
: Po i [a 
E=0 
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Type Operation 


(a) Variation of Net Flow from Valve to (b) Variation of Net Flow from Valve 
Cylinder with Pressure Differential to Cylinder with Valve Error 
Across the Piston 


Figure III-/,. Typical Variations of Net Flow 


III-10 


ee 


ie 


L_, 


— oJ wl cl 


} 


bk Sc 


J 


| 


— 


Section 2 


The figure also shows a typical operating point, the fixed values of the vari- 
ables defining the operating point being Q* ss ui ’ p* » and indicates possible 
values of the perturbed quantities. The partial derivatives in (III-11) are 
the slopes of the curves evaluated at the operating point, 

In Figure III-4, note that JQ/9£ is always positive and that 9Q/dR 
is always negative. It will be advantageous in the development that follows to 


make the normal signs of these quantities apparent; that is 
(III-12) gG= CoE - Cop 


where 





The term C,, is analogous to the slope of the torque-speed curve of a shunt 
motor and gives rise to a similar damping action. The term Cy, is analogous to 
the slope of the speed-field current curve of a shunt motor and produces a similar 
gain term in the over-all system. 

Cy» varies from a very small positive value to infinity as a function of 
fF and & , the very small positive value occurring when the valve error is sero 
and no pressure differential exists across the piston, and the infinite value 


occurs at a stalled cylinder condition. 


(c) FLOW INTO THE CYLINDER 

It now becomes necessary to relate the flow from the valve to cylinder 
motion, piston motion, and pressure differential across the piston. In devel- 
Oping the equations for cylinder flow, the compressibility of the fluid will 


be taken into account. 
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The total mass of fluid to one side of the piston at any time ¢ is 
(I1I-13) m=p97 


where 
0 is the volume occupied by the fluid mass 
fp is the fluid density at the time ¢ , assumed independent of position 
within the volume 


Then 
dm. (52% ,y~# 
(III-14) ri oes Lore 


and the volumetric flow corresponding to this mass flow rate is 


(III-15) OQ. ==. = 


The compressibility of the fluid is characterized by the expression 


dP dP 
eae = Ya ¥|m=constant “~ deo 


where 
N is the bulk modulus of the fluid 
fo is the fluid density 
P is the fluid pressure 

By inserting this relation into (III-15): 


dr, % dP 


eee Ee es ae 


(III-17) applies to each of the regions of the cylinder separated by the 


piston (see Figure III-2). 





AY, x AP 
(III~18) Q, —{ + N, r 
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If it is assumed that the instantaneous flow into one side of the cylinder 
is equal to the instantaneous flow out of the other side, 4i.e., Q=-Q,2 Q, 
and if it is further assumed that the bulk modulus NV is constant, and the same 
magnitude for both regions, (III-18) and (III-19) becomes 


(111-20) q- 2%, Uae. (th, & dB) 


t Ndt/ 
Since 
dy. aR 
dt at 


it follows that 


S 


df 
= geen, tee al 
(I-21) Yt = - BT 


Addition of the quantity ( % dP / dt) to both sides of (III-21) gives: 





ap. wt 4i,_5 2 _% 
(TI-22) at FE ae (PB) y+ de 
where 


Pp 2=P-8, the pressure differential across the piston 
Further, the time rate of change of the volume on the high pressure side 


of the cylinder is 


AY, 


a -_ 
(111-23) 7 = A-(X. xX, ) 


where 
X, is the cylinder displacement relative to structure 
X, is the piston displacement relative to structure (non-rigid piston) 


A is the cylinder area 
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By substituting the values of dP/dt and d x /dt from (III-22) and (III-23) into 
(III-20): 





a //%% \aP 
(III-24) Q2=A , (Xo X,) + if ) 








RACE 
or 
(III-25) Q = A 2 7 (Xo- Xp) + a 
where 
y= ree y? the effective oil volume* 


It should be pointed out that this analysis is strictly applicable only 
when 2 = 2%, or when the ports are blocked, since these are the only conditions 
where dy /dt = -di,/dt « However, the equation will represent an excellent 
approximation for any valve and cylinder where the flow through the valve ports 
due to compressibility is negligible, i.e., for any system where the valve 
decouples the effects of flow in the lines from flow into and out of the cylin- 
der. 

(II1-25) is written in terms of total quantities, i.e., operating point 
values plus perturbations. The similar relation in terms of perturbed quanti- 


ties is of the form: 


a(x xp) (7° “, ar") ad 1 (ap, 
(T1I-26) gy = A athena aes Gamma”? Gee dt Ar’ 


ix 
where 9» X%)» Xp» 7, » and Ayare perturbations, and 7” and (dP /dat)* 
are values at the operating point. 

ls 
The nonlinear term (4 ¥/N\ da /dt) is of higher order and may be ignored 


in a linear analysis, It is further assumed that (dP Jd t ye is relatively small 


* The physical representation of a column of oil as a spring also involves 
this same effective oil volume. 
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enough eo that the term(//N\dP/dt)“ AY can be neglected, With these sim 
plifications, the flow equation in perturbed quantities becomes: 


d(x. - Xp) ‘ Yan, 


III- = oe. 
(111-27) g = A dt N “at 


or, by designating s as the derivative operator, (d/dt) : 
abe 
(IlI-28) g = As(x,- xp) + oP 


To facilitate combining the flow and force equations, it is helpful to 
rewrite the coefficient x, N in a mre usable form. Since y"4s an area, 
and VV a bulk modulus, it may readily be verified that the quantity YUN has 
the dimensions of an area squared over lb/in., that is, of an area squared over 
a linear spring constant. The relevant area in this case is the net piston 
area, A. If the equivalent spring constant is represented by the symbol 4‘, , 
(III-28) may then be rewritten as: 


At 
(III-29) g = As(x,- x») + rede 


se 


where k, = y 


Thus the effect of compressibility of the oil appears as an equivalent spring 
rate, 
By eliminating the perturbed flow g from (III-12) and (III~29) there 


results: 
At 
(III-30) CeE = As(x,- Ap) + 4s ? Ce v2) 
This equation and the definition of the valve error, (III-1), are the basic 
relationships employed in the analysis which follows. 
For the sake of simplicity, the preceding derivation of flow relationships 


has been limited in such a way that it rigorously a, lies only to an actuator 
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with a symmetrical, closed center valve, with equal volumes on each side of the 
cylinder, and with equal net areas on each side of the piston which operates 
through small displacements. The analysis is relatively insensitive to deviations 
from these conditions, however; and it scope of applicability is quite wide. 

It should be mentioned that there are systems for which the analysis as 
presented here is not a sufficiently accurate description to permit useful appli- 
cation. The following are examples of conditions under which this statement 
may be true. If the valve is definitely an open center type, it is possible to 
have interaction between the actuator, and the tubing and pressure regulators 
of the supply system. Even if this interaction is taken into consideration for 
an open center valve, the analysis as given may lead to insufficiently accurate 
results because of the pronounced effect of neutral leakage. Space considerations 
sometimes demand that the piston have no blind rod or that the volumes on each 
side of the piston be markedly different for other reasons. There are some 
hydraulic fluids in use which display a constancy of bulk modulus only at fairly 
high pressures, The higher operating pressures can be obtained by purposely 
designing the valve to be asymmetrical. It is also possible to compensate for 
unequal net areas on the two sides of the piston or for preload (paired actuators) 
by designing the valve to be asymnetrical. In the case of unequal volumes on 
either side of the piston, this cannot be done, however. With any of the devi-e 
ations from complete symmetry mentioned, non-linear behavior is accentuated; 
this leads to chopped and offset waves in the frequency response. In almost any 
case, it is possible to extract a linearized solution which predicts critical 
frequencies and stability, using in general a complete double set of parameters 
Cy,Ce,4,A. This section has provided the initial basis for such a gener- 
alized development; however, it has not been carried through because of the 
specialized nature of its application. The analysis as presented is adequate 


for most design. 
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SECTION 3 = FORCE AND LINKAGE EQUATIONS OF THE ACTUATOR 


(a) DISTRIBUTION OF PARAMETERS 

In the development of the actuator analysis, it is-convenient to regard 
the hydraulic servo actuator as being terminated by two linear springs, 4; and 
k, » Linear motions of the ends of these springs represent inputs to the servo. 
Similarly, forces collinear with these motions are inputs. However, there must 
be an over-all dynamic balance of the forces acting on the system, so that one 
of them is not independent, and the linearized servo behavior may be described 
in terms of only three variables, One of the three variables must be a force, 
or convertible to a force. With these remarks in mind, the force and linkage 
relationships of the servo are assumed to be those implied by the schematic in 


Figure III-5. <A comparison with Figure III-1 is suggested. 
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Figure III~5. Identification of Lumped Parameters 
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The vertical rod of length /7 is assumed to be rigid. The mass \/-, at 
the free end of the vertical rod, is a lumped equivalent of the mass of this 
rod, of the slide valve, of part of the piston rod, and of part of the cable or 
pushrod between the controlling element and the actuator. 8. is a viscous damp- 
ing coefficient; it represents the effect of the damping in that part of the 
system associated with M, , due to velocity relative to the fixed structure, 

The magnitude ™ does not enter explicitly into the actuator equations; 
hence the radii from the fixed point of the vertical rod to the pivot points of 
actuator and valve are designated as the products of m by the dimensionless 
ratios 4 and q. The small error introduced by the fact that the motions of 
these pivot points are angular rather than linear is neglected. 

& is a damping coefficient representing the effect of relative velocity 
between valve slide and cylinder. Mp 4s the equivalent mass of the piston, and 
M. the equivalent cylinder mass. B. is a damping coefficient representing the 
effect of relative velocity between piston and cylinder. 

A force Ap, on the piston and cylinder acts as shown in the diagram. 

As a matter of terminology, points at which displacements such as X, and 
X, are measured will be called "nodes." The appropriateness of this terminol- 
ogy may be better appreciated after considering the mechanical network diagrams, 
or "nodal diagrams," used later in this chapter and throughout the remainder of 
the volume. For a more complete discussion of such diagrams, M. F. Gardner and 
J. L. Barnes, Transients in Linear Systems (John Wiley & Sons, N.Y., 1942), may 
be consulted, 

The nodes x; and X, may be chosen so as to facilitate analysis of a larger 
system of which the servo actuator is only a part. It is to be noted that the 
selection of these nodes determines what are to be considered the coupling 


springs A; and A, . If numerical results of the hydraulic servo analysis are 
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to be supplied to another group, to be used in further, more comprehensive ana- 
lyses, these nodes and consequently the springs should be definitely identified 
in the physical system; doing this imposes no restriction on the way in which 
adjacent sub-systems are analysed. If, on the other hand, a wholly literal analy- 
sis of a more comprehensive system is to be performed, Xx; and X, may well be 
chosen at the terminal nodes of adjacent sub-systems. 

In general, it will be necessary to increase the values of the quantities 
Me » 8: » and M; indicated in the preceding discussion to take into account 
distributed masses and dampings in coupling members. 

The forces fy and f; each arise from inertias, dampings, and applied 
external forces associated with one or more degrees of freedom of an adjacent 
system, 

(b) FORCE AND LINKAGE EQUATIONS 

The equations of the force system of Figure JII-5 are obtained by writing 

out the expressions for force and moment balances in a straightforward way; 


only the results are given here: 
(III-31) Fy thx, = ApXs 


(III-33) (M.s* + Bs + k.)x,~&s x, -(B,s tk JE-AQR = kay 


(11-34) (Mps"*+ Bs + kp) Xp $Bsx,-hbx, +Ag = O 


(III-35)  (M, s*+ Bs + ky + bhp) x, - bkpXp + Q(AS th JE = kzx; 


(III-36) %, = AX, 
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SECTION 4 = COMBINATION OF FLOW, FORCE, AND LINKAGE EQUATIONS 


(4) SIMPLIFICATION OF THE EQUATIONS 

All of the equations defining behavior of the hydraulic servo actuator have 
now been derived: (III-1) defining the valve error, (III-30) describing flow 
behavior, and (III-31) through (III-36) expressing force and linkage relation- 
ships. In all, there are eight equations with ten variables  ,/; , %) » Xp>5 
ams» %pe ss » and E . All these equations are linear with coeffi- 
cients which are algebraic functions of the derivative operator S., The over~all 
servo behavior can then be described by a single literal equation involving only 
three variables, When the equations are written in determinantal form, the coef- 
ficients in 9 may be manipulated in the same way as constants. However, if the 
equations are treated in the present form, the geometrical factors and flow coef- 
ficients make the reduction difficult. If the equations are changed to a form 
which can be easily expressed in a nodal diagram (mechanical-network diagram), 
the reduction is greatly facilitated. 

An important simplification in the force and linkage equations (III-31) 
through (III-36) can be effected by relating the effects of masses, dampers, and 
springs acting at a radius om , to those of equivalent quantities acting at 
some other radius 6 ; for example, the equivalent mass, |/ ‘ » acting at a radius 
&m is equivalent to the original mass, MW , at a radius om multiplied by 
o(/p* ;i.e.,M= M(<'/B*). This simplification also applies to the 8's and 
k's e On the other hand, the displacement ~ at a radius pm equivalent to 2 
atxmis f x foc ‘ 

The force and linkage equations (III-31) through (III-36) may then be sim- 


plified by the substitutions: 





Xm- = DX» , Arp = bX 
_ M: B; ky 
Mp b? 9 B;, 3 re) 9 Kip 2 “pe 
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Also, the node x, is easily eliminated by combining Equations (III~-1) and 
(III-36). With these simplifications the set of equations (III-31) through 
(III-36) take the simpler forms: 


(III-37) £ + kx, = hax, 


(III-38) ey Kip Xap ™ Kip Rep 
(11-39) (Ms' + Ast k,)x,-45%m- (Bs th JE- AR = koX%, 


(III-40) (Mps*+ As + ky)x, - Bx, - kpXm. + AD, = O 
(III-A) (Myst Bg + kip + Kp) Xmp ~ KpXp * = (Bs thE = hip Xp 


(1IT-42) € = - Kme~ Ke 


The flaw equation (III-30) may be rewritten as: 





& 
[(Z)e + 46] 


It may be readily checked that AC, / Cy has the dimensions of a spring constant 
£ 

(1b/in.) and that A /Cy has those of a viscous damping coefficient, These 

equivalent parameters of a force system are designated in the following as ky 


and & , respectively. Equation (III~-43) then becomes: 





(III-44) (By s # A)(52) 8, 3(x,- Xp») - kg E = O 


where 
PY) 
= ACe _ _, 9F L of 
hy Ang PAS 
dF 
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The nature of the "flow damping" B. can be seen by rotating Figure III-4a 
through 90° clockwise, The "flow spring constant" 4, is the slope of the cross- 


plot shown in Figure III-6. 


m | 


B= (Ph —Ps ) Q =Constant 
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Type Operation 


Figure III~6. Variation of Pressure Differential Across the Piston 
with Valve Error 


The definition of the parameters 4, and &. , expressing flow relation- 
ships, implies the presence of a physically important node associated with them. 
This node is located somewhere within the body of oil in the cylinder and can 
thus never be experimentally observed, but it is important in speeding the reduc- 


tion of the actuator equations. The node x, is defined by the relationship: 


(III-45) x, 





+ Xo 


= Ap, 
k, 


° 
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With the introduction of xa, , (III-44) becomes: 
(III-46) (8, s + ky) X%, -— Bs x%_- kx, -kKE =O 


The force Ap is then given by: 


(IIIW-47) Ap, = ky (x,- %) = BS(Xp- HX) + KE 





Figure III-7. Partial Nodal Diagram Expressing Flow Relationships 


The partial nodal diagram shown in Figure III~7 may be constructed from 
Equation (III-46). The only node completely described by this diagram is the 
node x, ‘This form may be used directly to express the flow behavior irrespec- 
tive of linkage configuration, It is interesting to note that if p, is sero, 
X% and 2%, are equivalent. The displacement x, may then be considered as the 
output of an actuator subject to no external or internal loads. This nodal dia- 
gram is of partiediar use in the physical interpretation of the actuator as a 
force source driving x, through a spring 4,. All of the parameters f£, , S, ; 
and 4, have physical values which can be measured by test. 

The incorporation of (III-46) and (III-47) into the set (III-37) through 
(III-41) results in the following: | 
(III-48) A+ Am, = kody 


i 


(III-49) rs ss kip “mp ™ Kp ip 
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(IIT-50)  (M.s"+ Bg +k, + hy)x, ~ B.S%p ~ hy Hy (BS HK )E ™ KyXy 
(111-51) [Mps"+ (B+ Als + hy]xp- Asn, -Bes%y- bpm tHE = O 
(IITI-52) (M:, g*+ By. 3 + kip + kp) X&pp ~ kpi,» + <(6,s tka )E = Kip Rip 
(ITI-53) (4.5 * ky) x, ~ Be Shp ~hoX% -~KpE = O 


Equations (III-50) through (III-53) are now all nodal in form (see Figure 
III-8), and the reduction is greatly expedited. The only force sources are in 
the valve error or its time derivative, where the valve error is defined: 


(IlI-54) E& = = hme + x, =O 


(b) REDUCTION OF THE EQUATIONS 
The analysis presented in the following chapters depends on the reduction 


of the controlling equations of the actuator in such a manner as to isolate four 


relationships: 
iy Dx; x; ) 
F; SAX x: ) 
(III-55) 


Xo - EE > XK, ) 


& = $4(x: Xe 1 Kg ) 


The first three of the functional forms given above are unique. The fourth 
form will be determined in such manner as to facilitate physical interpretation 
of the results. 


The determinant of the coefficients on the left of (III-50) through (III-54) 


III=25 
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Solving (III-50) through (III-54) for x, and X,, , each in terms of X, 
and x; , and substituting in (III-48) and (III-49), 


(111-56) £2 Ae%elAa@ he Ay) = hip Rip (he Any) 
A, . 


= k;,. Kip (A, = kipAyg)— ko Xs (kipA,s) 


A, 


(I1I-57) 


&|an 


where A;, is the (signed) minor of the element in row ¢ and column J of the 
system determinant, and 4, is the determinant itself. For purposes of dis- 
cussion, it is convenient to have the above equations expressed in terms of 
operators which reduce to unity, or a good approximation thereof, when 3 
equals sero. This is accomplished by dividing numerator and denominator of 
each equation by the product of major actuator springs 4,4,4,4;,Akp» A, 
being omitted because it is small in comparison with A, . The equations may 


now be rewritten: 


(rnsa) fg = Ae 2lO%e) 
Ci 
k 
fs 


(111-59) hw BOX) HAs 
. {Se 
k. 


where 
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(F)4u) 


G2 = 


k, kp ky 


b 
( = \o — Mcp Ass) 
93 * k, ky Kp k 


= Me 
4 ko kph, 
Gs — A, 


Ke 7 k, kp ke k, ki, 


J» 9g» G+ Gs» and Js are polynominals in S$, which reduce to unity 
when S equals zero if the input coupling spring A, ie much smaller than the 
flow spring kK, « 

One additional operator of some use in later work may be derived from the 


system determinant. Equation (III~56) may be written in the form: 


k Ay | 
ol Ke ip Kinks , 
(111-60) & = Ak ee ile db x, - 92 (a%;) 
($s 3.) 
ke kK, 


An expression equivalent to this is 


e 6 e 


(4) 
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where 
a 
/ / b | ab 
—_  — = ~ —— 
Kee k, ki, k, kt 


g. = Ai, Kyi 


and g, is an expression in gs which reduces approximately to unity when Ss 
equals sero. The third required relationship of (III-55) is derived by taking 
the determinant of the system of equations consisting of (III-50), (III-51) 
with the node <,,, eliminated by use of (III-54), and (III-53): 


Ms" *As thik, ~89 2 
~(Bs +S hy) Myst +(e Bs + kp 8, ¢ 
~k, ~B.s Bs + k, 


The solution to this system of equations is given by: 


(Bs ke)B, + [ky Ld hy) Br + k, By, 
| a | B,, 
E+ |rs 





(III-62) x, 2 $2 8s “pt el 
Sn 
AF | 


where A, is the determinant of this system and the by its minors, as before. 
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The fourth functional form in (III-55) remains to be determined. If (III-50), 


(III-51), and (III-53) are added, there results: 
(111-63)  -(B,s +k, )E + Mp3" Xp ~(MpS** kp) Xn ~ Xp) = ke X%y — (MS k,) x, 


The determinant of the system of equations (III-52), (III-54), and (III-63) is: 


a 
> (6s +k.) M:.s"+ Bs + kip ~kp 

/ as O 
—-(B,3 + k,) Mp s* —(Mps* + ky ) 


By solving these three equations for E&: 





2 Cy 2 2. C3, 
) (ax;) - (Ms + k)Cy~ Ce, x, + 8 dbk. Ky 
(111-64) € = +2 ae aoe ee! 
ee 
Kp Kip 


In expanding the determinants of this chapter it is worthwhile to exploit 
certain generally applicable relations between the parameters. Since the equiv- 
alent flow damping 8- is inversely proportional to the valve flow coefficient 
C, » which remains quite small, & is always large in comparison with 8B. , 4, , 
and 4, . Also, the input coupling spring 4, will always be negligible in com- 


parison with the flow spring 4, . 
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B, >> & 


B, >? (a-b) B, 


6, >> B 


A, >> k 


The desired relations of (III-55) are now established. 

The equations derived in this chapter are very general. In Chapter IV, 
these equations are restricted by imposing the load of an aircraft control 
surface on the hydraulic control system. Chapters V and VI consider two speci- 


fic types of hydraulic surface actuator systems. 
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CHAPTER IV 
A GENERALIZED HYDRAULIC CONTROL SYSTEM 
SECTION 1 - SIMPLIFIED CONTROL SURFACE CHARACTERISTICS 


A general *nalysis of control surface behavior is outside the scope of this 
volume, since this would require that the reader have a rather comprehensive 
knowledge of aeroelasticity in order to understand the complex interaction possi- 
ble between non-rigid parent and control surfaces. 

It is usually possible to represent control surface behavior fairly accu- 
rately by including with the actuator analysis a mass M4 and a coefficient of 
viscous damping 8. . ‘These control surface paraneters are referred to the line 
of action of the actuator cylinder, ‘The control surface, regarded as a flexible 
etructure, has some effective spring constant which expresses the flexibility 
between its effective mass center and the surface horn. The coupling spring ke 
of the actuator analysis is conveniently taken as including this surface flexi- 
bility, The displacement *%, then becomes a node at the effective mass center 


of the surface, The nodal diagran is show in Figure IV-l. 
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Figure IV-l, Nodal Diageran of Control Surface 
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The force F of Figure IV-1 is a quantity which is predominantly the aero- 
dynamic moment acting on the control surface, referred to the line of action of 
the cylinder; however, this variable may also include terms expressing a coupling 
between multiple modes of vibration of control] and parent surfaces. The equation 


corresponding to Figure IV-l is 
(IV-1) (Ms + Bs + k,)x,- kx, 2 F 


An equivalent relationship in terms of the force ( existing in the coupling 
spring A, is secured by use of Equation (III-48). 


(Iv-2) (M,s"+As)x,+ fh =F 


It must be remembered here that the spring constant A, is understood to include 
surface flexibility. 


SECTION 2 ~ A GENERALIZED HYDRAULIC CONTROL SYSTEM 
It ie now proposed to consider the system composed of the hydraulic servo 
actuator of Figure III-8 and the simplified control surface of Figure IV-l. The 
final equations of Chapter III will be so combined with (IV-1) and (IV-2) as to 
permit investigation of the behavior of the system as a component, and also of 
the interaction through forces between the surface~actuator system and adjacent 


systems. The results of the combination will be a set of five relationships: 
Fe 2 (x, , X;) 
f= 8 (x; F) 
(Iv-3) x," Z,(% , F) 
x,=#F,(E ,F) 
E= 2,(x,, x. ,F) 
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The first relationship of (IV-3) is of use in the investigation of a gen- 
eral flutter problem in which there is interaction between the control system as 
a servo and the aeroelastic behavior of control surface and parent surface. Parv 
II of this volune is concerned with the statement of the flutter problem where. 
there may be interaction with servos of the aircraft control chain. The rela- 
tionshin for f in terms of xX, and xX; will be expressed in three ways, each 
of which lias some useful application to servo-flutter coupling interaction. 

If the behavior of a system at the input end of the actuator is markedly 
affected by the force exerted on it by the actuator, there is force coupling, 
as well as displacement coupling, at the input end of the hydraulic control sys-~- 
tem. In this case it is necessary to use the second relationship of (IV-3). 

The third, fourth, and fifth relationships of (IV-3) represent the behav- 
ior of the hydraulic control system as a positional servo. The over-all, or 
closed-loop behavior is given by x, (which for convenience of comparison with 
test results is taken as the output of the system) in terms of A; , the displace- 
ment input to the system, and fF , the disturbing force. The fourth and fifth 
expressions of (IV-3) are functional forms of greater utility in the analysis 
than closed-loop expressions. Isolation of the critical variable E , which is 
tia valve error, gives rise to an open-loop transfer function when the disturt--. 
ing force F can be regarded as negligible. 

the relationships of Chapter III which can be combined with (IV-1) and 
(IV-2) to determine the first three relationships of (IV--3) ure (III-58), 


(III-59), and (III-61). They are repeated here for convenience: 


(IV) f= 2%s_— 92 (ax; ) 


(7 


lVv—2 
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(IV-5) 


(Iv-6) 


where 


Fz _ 9Gslaxi) — %4%s 
me) 
ky 
FE) 
k ad é x 92 (ax; ) 
3 red S 
Cea 
k. k. 
a 
/ / / b | / ab 
— Bias pon a eS oS 
k, k. ky Kip k. ky k; 


<2 
and J,» I 9 Sg 9 Gy J, are polynomials in S , each of which reduces to unity 


when S equals zero, if the input coupling spring ke is much smaller than the 


flow spring rate A,. By elimination of / between (IV-2) and (IV-4), 
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By elimination of £ between (IV-2) and (IV-6), 


tak 
sel ° bo 9, (axz) 


(i) J %) 


(Iv~8) F = |(M,s°+ Bs +k.) - 


By rearrangement of (IV~7), 


wr aoltl g 
-_ -- E +(Ms + 9 7 = wie 


i 5 
k, hk. 
By elimination of x, between (IV-5) and (IV-9) 


p| Melee! + (M, ot b,#)%5an:) - bg, F 


_ Ys 
(IV-10) F; E + (M,s** Bs (2)) 


By elimination of ~, between (IV-1) and (IV-9), 


Qe (a st+ ess Vax.) + (Age gs. Kei 9 & cal 
k kee 

(1v=11) x. a k. ke ke Cc Fe 

E + (M,s*+ Bs ra 
¢ 


It is of interest at this point to examine the forms assumed by (IV-9) and 


.(IV-11) when gs is equal to zero: 


F 
TOR 


é 


(IV-lla) Xo ka O 


1V--5 


Scction 2 


= (ax;) + a 


(IV-9a) ; 
¢ 


*s |s =O 
It can be seen from (IV-lla) that if the force - can be neglected, the con- 
trol system is a servo with no steady-state position errors i.e., it has a single 
integration in the open-loop transfer function. The output displacement X, 
attempts to follow the input displacement referred to the line of action of the 
valve slide. With £ left undetermined, there is a position error determined by 
the spring ke from the input node x; to the output node x,. The spring ky» ; 
is the series combination of the flow spring A, and a modified input spring 
k; /ab . If the power boost ratio / /b is infinite, the spring k,, becomes ky, . 
It is important to note that the servo action determines the steady-state behav- 
ior, and Ay and 4, are not involved. Equation ([V-9a) shows tiiat in the steady 
state X,, which is the displacement of the effective mass center of the control 
surface, behaves somewhat as the cylinder displacement X, , except that the 
effective spring determining position error is now A, , the series combination 
of k,, and A, « The position error due to a static force fF is of course 
greater at ~, than at ~,. 
The fourth expression of (IV-3) is obtained by combining Equations (III-50), 
(III-51) with the node ~,,. eliminated by use of (III-54), (III-53), and (IV-1). 


The determinant is 
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The solution may be written: 


(IV-12) A= EE + EF 


where 





b 
ba t ky oy, + (ky - m Mp) 21 + Kye Xz, | 


_ ke Ko kp ke 
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k, Ke Kp Ke 


The fifth and last expression of (IV+3) is obtained by using (III-52), 
(III-54), (III-63), and (IV1). The determinant is 


(Bs + k,) M;.s°+ B..s + k;, ky 
a 
/ b 
2 
-(Bs + ke) M, s* -(Mps+ ky) 7,2 


M,s*+ Ast k, 
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The solution may be expressed in the form: 


(W-13) EE = YO [Y, (ax) - x + %, FI 


E 


where 


ae ke a ~ (3a, + fee | 
k k Kip 
Ve —-. 


$ Ee * LY) Ls 
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Equations (IV-12) and (IV--13) are the basis for the block diagram shown in 
Figure IV-2. 

The quantity F of (IV-11), (IV-12), and (IV-13) in general includes terms 
which express a coupling between multiple modes of vibration of control surface 


and parent surface. If the maynitudes of these coupling terms are the same as, 


IvV-8 
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or less than, that of the true aerodynamic force, the tot*l * composed of 
aerodynamic force and coupling terms can usually be neglected without severely 
affecting the system analysis. If the coupling terms are large, however, the 
system analysis is incomplete, and the results of the actuator analysis must 
be combined with a more sophisticated formulation of the characteristics of 
control surface behavior. Such difficulties are often associated with uncon- 


ventional control surface design, such as all-movable surfaces. 
F F 


Figure IV-2. Block Diagram of the Generalized Hydraulic Control System 


If the force * is considered negligible in Figure IV-2, the open-loop 

transfer function of the positional servo is given by the product of %eg and 

% « From an investigation of the frequency response of the open-loop 
transfer function one can form an excellent idea of the characteristics and 
limitations of various systems. 

The following two chapters are concerned principally with the frequency 
response and stability of two systems. Chapter V deals with the case that the 
power boost ratio //6 is infinite and the input ratio @ is unity. Because 
of the wide use of the "fully-powered" system, portions of the material of 
Chapters III and IV are repeated, so that Chapter V may be read with little 


reference to earlier chapters. In Chapter VI analysis of the generalized 


power boost case is resumed. 
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GHAPTER V 
THE FULLY-POWERED HYDRAULIC CONTROL SYSTisM 


SECTION 1 - THE GENERALIZED SYSTEM 


The purpose of this chapter is to present the equations governing the 
behavior of a fully-powered hydraulic control system and to investigate the 
effects of its various parameters upon frequency response and stability. The 
schematic diagram is shown in Figure V-1, and the corresponding nodal diagram 
is shown in Figure V-2. The form of this diagram is derived from Figure III-5 
by assuming the power boost ratio \/b to be infinite and the control input 
ratio 2 to be unity, and by adding the simplified surface characteristics dis- 
cussed in Chapter IV. 

The viscous damping coefficients 8. and &- represent respectively the 
damping due to the relative velocity of the valve slic. and the cylinder, and 
the damping in the cable and feel device. The mass M; is the effective mass 
of the valve slide, the quadrant, and the cable. The spring constant 4, rep- 
resents the effect of a springing action between the valve slide and the cylin- 
der; and 4; represents the effect of the flexibility between the valve slide 
and a point of the controlling system having a displacement x, , which is the 
actuating quantity of the control system. 

For purposes of classification in this volume, a hydraulic control system 
of the general type shown in Figure III-5 is regarded as fully-powered whenever 
the power boost ratio //b is infinite. Strictly speaking, however, the damp- 
ing 4, and the spring A, cause a deviation from fully-powered operation in 
that a force may be transmitted from the control surface back to an initially 
controlling system, which may be the human pilot or an autopilot. If the behav- 
jur of the pilot or autopilot is markedly affected by a force transmitted from 
the control surface because the value of 4. or of A, or of both is other than 


zcru, there is a possibility that foree coupling, in addition to displacement 
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coupling, exists between the initially controlling system and the hydraulic con=- 
trol system. With the human pilot, force coupling is probably not very important; 
however, with some autopilot installations, the effect might be appreciable. In 
most instances it is sufficient to consider only displacement coupling, and the 


following devclopment is based on that assumption. 
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Figure V-l. Schematic Diagram of the Fully-Powered 
Hydraulic Control System 


The controlling equations are easily written from Figure V-l or Figure V—2: 


(V-1) (M.s°* (B+ BJs +k, * k.]x, - BS %p- B39 %,~ kot, - Ag xy, ™ AD, 
(V.-2) (Mp 3° + Bis 7 Ap ) Xp - Asx, = ~Ap, 

(V--3) (Mis'* Bost ki dry hex = F 

(V-4) [M.5°+(B:4+ &)s the + kx, ~B9x,- kx 7 Ap Xe 
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The flow relations governing the force differential across the piston, Ap, , 
due to the pressure differential D , have been treated in detail in Section 2 
of Chapter III. The equation for A p, is derived from Equation (I-30) 

At 

(V5) CoE = Ag(x,- Xp) n(4s - c,) 
where 

A is the net area of the piston 

A, is the spring rate due to oil compressibility 

E ia the valve error, E = x,- x, 
and C, , Cy are experimentally determined coefficients of the net flow from 
valve to cylinder as a function of valve error E and the pressure differential 
2, across the operating ports, that is, across the piston. When multiplied by 


A/Cp » Equation (V=5) becomes 


(v~6) (4). + kal “2 - (Eee (A) atx X») = O 


3ing, for convenience, equivalent force parameters defined by: 


ky ae ACs 
Cp 
A 
Cp 


the equation (V-6) takes the form 





(es + kyl AZt]+ Bp a(ny~ Xp) ~ he E = O 


The incorporation of the val.ve error € and (V-7) in the set (V-1) through 


(V-4) results in: 


(V-8) (M. s°+ Bs tk.) xo- Bog %p ~ ky xy - x, (92) =(B,st A)E 
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(V-9) (Mp s? + Bs + kp )%> -~-QoSx, thy (42) = O 
°o 





(V-10) (B, + ty 5) * Box, ~ Bsm = kK E 


(V-11) (M, s° + By S + ke) Xs ~ Ke%, = F 
(viz) [Mg s"+ (Be + Bat ke + Bl E + (Mest Bo + hela, = kgm; 
A simple and readily interpreted result is obtained by solving (V-8) through 


(V-11) for x, in terms of € and F and Equation (V-12) for € in terms of ~; 


and x, + The system determinant of the first four equations is: 


M.s'+ Bs +k, -Bs -k, -k. 
-Bs MstBs t kz k, O 
Bs “Bes B.s+k, : O 

-k, O O Ms +B,s + k, 


This set of equations yields the following expression for the cylinder displace- 
ment X,3 


; (4,5 + ke )A, + ky As, | E.+ (-A, A] 
V-13 <= * 


where A is the determinant shown above. The coefficients in (V-~-13) are poly- 
nomials in the derivative operator S. It is convenient for discussion and for 
forming of Bode plots to have the numerator coefficients in such a form that 


sach reduces to unity, or very nearly so, when s is set equal to zero. This 
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is accomplished by dividing numerator and denominator by the product of the 


spring constants, k,. k, k, k, « In symbolic form 
(V-14) Aa = Ee er 


where 


[Pus + ke)A, re As, 


Y= ky ke kp ky ky kp ke 
Li heehe. 
k, ke kp ke 
y = Ke LK hp he J 3 kp ke. 
fete) ‘e) k, kp “I 


When solved for the valve error E€, (V-12) becomes 


k;%; a (M;s" + Bs + ki )%o 


Vel moe 
(V-15) E [M-s*+ (B.+ B.)s ? k; + Kel 


A functional form involving a unity feedback from the output x, is desirable 
to facilitate interpretation of the results. The coefficient of ~, in (V-15) 


is then taken out as a multiplying factor. 





(M; 5" Bs + k;) en 
(V-16) E M:5**(B,* B, )s + k, + k, (M, s +B + k;) x 


If the static gain is isolated in (V-16): 


(V-17) E = Leg 
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where 
M2, & | 
¥ = (- k; 5 = k; tie 
BM tkel]_ Me oe, Bt h 





k; ? ke k; + k, 





(V-14) and (V-17) are readily combined into the functional block diagram shown 
in Figure V-3. 





xy 
Figure V-3. Block Diagram of the Fully-Powered Hydraulic 
Control System 
It is now necessary to consider the quantities involved in the expressions 
for A and A « Certain of the relations between parameters are of such gen=- 


eral validity as to justify basing all discussions on simplified expressions for 
these quantities. Since the equivalent flow damping LB. is inversely propor- 
tional to the small valve flow coefficient C, , & is always large in comparison 
with the other damping coefficients 8. and 4. Also, the control surface mass 
M, will always be large in comparison with the cylinder mass M. and also with 
the piston mass M4,. The input coupling spring A, will always be insignificant 


in comparison with A, . 
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B, >> B, 
A >> & 
ky >> Ke 
My >> Me 
My >> Mp 
If the simplifications indicated by these relations are made, the expressions 


for ); and ) become less complex. 
The numerator of % is 
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To sketch out the frequency response of the fully—powered system, approxi- 
mate factors of the polynomials in Equations (V=18), (V-19), and (V=20) must be 
determined. These factorizations are obtained by taking into account relative 
magnitudes of parameters which are typical of practice. 

An important parameter of the system is the surface mass VM, . To a large 
extent, general design considerations make it necessary that M, be fairly large. 
The cylinder mass M. and the piston mass (/, are always quite small in compari-~ 
son with A4 3 in addition, M, is several times as large as 4,. The piston 
spring Ag and the flow spring A, are each several times as large as the oil 
seine k, » The surface coupling spring k is typically somewhat larger than 
A, but still much smaller than 4, and A, . The cylinder damping & , the 
surface damping 4, , the valve demping &, , and the input damping B; are all 


limited to a low range of values, and for purposes of factorization, they may 
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B, >> B, 
A >> B, 
ky >> Ke 
M, >> Mp 
My >> Mg 
If the simplifications indicated by these relations are made, the expressions 


for ); and ¥ become less complex. 


The numerator of % is 
(Bs + ke)A, As, _ (Me » 8s ). 
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The numerator of 7% is: 
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To sketch out the frequency response of the fully—powered system, approxi- 
mate factors of the polynomials in Equations (V=18), (V-19), and (V-20) must be 
determined. These factorizations are obtained by taking into account relative 
magnitudes of parameters which are typical of practice. 

An important parameter of the system is the surface mass MV, . To a large 
extent, general design considerations make it necessary that M/, be fairly large. 
The cylinder mass M. and the piston mass (/, are always quite small in compari~ 
son with A4 3 in addition, 4, is several times as large as M,. The piston 
spring Ap and the flow spring 4, are each several times as large as the oil 
sonine k, « The surface coupling spring A is typically somewhat larger than 
A, ‘but still much smaller than 4, and 4, . The cylinder damping & , the 
curfuce damping 4, , the valve dimping 4, , and the input damping 8; are all 


Limited to a low range of values, and for purposes of factorization, they may 
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ve regarded as approximately equal in magnitude. The input spring A; is in 
general the spring constant of a long, light cable, and is usually very much 
smaller than 4, , A, , Kp, and A, The input coupling spring A, which can 
be used is limited by the force which it is desirable to impose on the pilot or 
the autopilot, and may be considered to be either of the same magnitude as &; , 
or less in magnitude than /; . 


If 8, and ke are negligible, the numerator of % , (V~18), reduces to: 


=(“%~.2, 5s | e 
(V-21) Mim y le 4 m0) fat an st/ hp gs +/ 
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the quadratic containing M, being effective only at a very high frequency. With 


6, and k, large and Ap several times as large as A, , a good approximation is 


| Msg & 4, Mp _ 2 a ) 
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the last quadratic being characteristically non-minimum phase. 


The numerator of ); , (V-19), has the approximate factors 


~ tpt, Mp ot, ] 
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From (V-20) it is seen that the denominator of Y¥ and & is composed of 
the factor (6./k,)s = (A/C,)s (which expresses the fundamental behavior of the 
control system as an integrating servo) and a sextic which consists of three quad- 


ratics in M,,M, , and My. In investigating this sextic, it is worthwhile to 
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consider first the undamped behavior with 4. and 8, each equal to zero and 
with &. large. With typical magnitudes of the parameters, the approximate criti- 


cal frequencies in this case are: 


(IV-24) 





Typically, wW,“w,<4), Inclusion of the damping terms results in only minor 
modifications to these natural frequencies. From (V-21), (V-22), (V~23), and 
(V-24,), it is seen that only a small amount of damping suffices to cause the 
effective cancellation of Mp quadratics in the numerator and denominator of 


Ye and % . An approximate expression for the cylinder displucement can then 


be written: 
s 2h s’ 2g I [s 
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The frequencies and coefficients of the middle terms of the quadratics in (V-25) 


are approximately as follows: 
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SECTION 2 — THE SYSTEM WITHOUT INPUT COUPLING 
If the damping 8, and the spring A, in Figure V-1 have an appreciable 


magnitude, they effect a force coupling at the input. Section 3 of this chapter 
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discusses how these quantitics can be used to stabilize the servo. If, however, 
these parameters are not delibcrately augmented in design, their effect on sta- 
bility is negligible, and the system where 4. and A, are effectively zero may 
be considered a basic case. Then leg, reduces to unity, and the block diagram 


takes the form shown in Figure V—4. 





Figure V-4, The Fully-Powered System without Input Coupling 


The effect of airload is always in the direction of greater stability and 
is almost always small. Moreover, it is necessary to design a system which is 
satisfactory on the ground or at low airspeed. The force * can therefore be 
neglected. An open-loop transfer function can then be written with the approxi- 


mate factorizations of the precoding section: 
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A typical frequency response corresponding to (II-26) is shown in Figure V-5, 
Inspection of (V-26) and Figure V-5 brings out several important facts. 
The only effective damping in the quadratic containing M. is due to G , the 
damping coefficient pertaining to the relative motion of piston and cylinder. 
The damping ratio 4 is not large, and there is a fairly sharp peak at G), on 
the Bode plot. However, since the amplitude ratio is very low at this point, 
there is no danger of instability due to this peak. The damping ratio, g, P 
will never be large enough to cause a large phase lag in the vicinity of Ud), 


and G),. It is permissible, then, to disregard the quadratic in M_ and write: 
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Section 2 


where the frequencies are given in (V-26). The nature of the stability problem 
and the resultant limitation on performance are now evident. With a, differing 
from @), , principally because the oil spring k, is small, considerable damping 
£ is required to lower the resultant peak and thus prevent instability. The 
crossover frequency «) = k, /B.*G /A principally determines the speed of response 
of the system. It is desirable to make this quantity as high as possible within 
some specified limit of overshoot in transient response. The magnitude of W, , 
the "gain" of the system, which is obtainable in practice is almost always lin- 
ited by an incipient instability due to the separation of u), from G), mentioned 
above. 

The characteristic equation corresponding to the open-loop transfer function 
(V-27) is 


21 [Be Me + [a(t et) eh 
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Application of the Routh criterion to this relation yields the following nece~ 

essary condition for stability: 
M, 88 By B( L) (2 L | 
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A stable control system is not necessarily a satisfactory system. The over= 

shoots in transient response may be too large or the time of response too long 


to meet operating requirements. This is, however, too much a matter of indi- 


vidual design to permit a generalization here, 
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Figure V-5. Open-Loop Transfer Function Without Input Coupling 
and With No Appreciable Airload 
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If the damping ratio £, in (V-27) is sufficiently large, instability 
caused by the separation of @, from @, will not occur. The individual terms 
of 25, fu), thus become of definite interest in design. There is some small damp- 
ing because 4. is finite. The surface damping, 5, , is more effective as a sta- 
bilizing parameter than the cylinder damping, 4, , since the former acts through 
a weaker spring. 

It is important in view of weight and space limitations that requisite sta- 
bilization of the hydraulic control system be carried out with a minimum of 
design complication. Even though & must be increased by a large amount if it 
is to have much effect on stability, the cylinder damping is very important from 
the standpoint of simplicity of designs Considerable cylinder damping is obtain- 
ed in the course of preventing leakage past the pistor by use of O-rings. With 
piston diameter roughly established by general design considerations, 4, can be 
increased by inserting additional rings or by using pairs of rings with the inter- 
mediate region vented to the atmosphere. The latter technique increases the damp~ 
ing by increasing the pressure differential across each ring. 

The surface damping 8, arises from several sources, First, there is some 
parasitic damping at the hinge and surface horn connections. Second, damping 
devices are sometimes inserted at the hinge line. The latter method is not alto~ 
gether successful, since such devices usually take a heavy toll of weight and 
space allotments. A third source of surface damping is aerodynamic. With a fair 
degree of accuracy, the aerodynamic force acting on the control surface of an 
airplane in flight (/f in Figure V-4, or Figure V-5) can be represented as the 
parallel combination of a spring and a damper. The effect of the airload spring 
on control system behavior is not great, and no account of this parameter has 


been taken so far in the analysis. llowever, the aerodynamic damping is simulated 
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as a part of 8, , the surface damping. This portion of &, can be of appreciable 
magnitude, but it cannot be regarded as available for use in stabilizing the sys 
tem, since the servo mst be satisfactory on the ground. 

It is appropriate at this point to discuss briefly the effect on servo 
behavior of the springing action of the airload on the control surface. The 


force F of Figure V-4 is assumed to take the form: 
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If the airload spring A, is included, the fundamental character of the 
servo is disturbed to some extent, as is indicated by the last term in the demon—- 
inator of (V-23). However, this additional term has little effect on servo behav- 


ior because a typical flow spring constant, A, , is from 100 to 1000 times as 
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large as a typical airload spring constant, A, . For all practical purposes, 


the extra term may be omitted. The frequency response then assumes the form: 
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If a system with the airload spring is compared, on the basis of the same 
open-loop gain, with a system where this spring is negligible, the following 
remarks regarding the effect of A, may be made. The damping ratios of both 
numerator and denominator quadratics are virtually unchanged. The frequencies 
G, and Gs; are increased, but their difference, W,—- GZ, , is slightly decreased. 
There is then a slight stabilizing effect due to the airload spring, and thus 
calculations which ignore the airload spring are somewhat conservative as far 
as stability is concerned. A further justification for omitting the airload 
spring in analysis is that it becomes weak at low airspeeds. 


The steady-state servo response, including A, , is given by 


"- ky - As] 
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Since A, is normally much smaller than A, , the position error is small. 

The hydraulic control system must ultimately be considered as a component 
of a larger system including the airframe. In typical practice the hydraulic 
control component is coupled to the more comprehensive system only at very low 
frequencies. It is then possible to view in isolation the question of stability 
of the hydraulic control system. If the hydraulic control system itself is sta- 
ble, its behavior as a component of the larger system can usually be simplified. 

If the numerator and denominator quadratics are deleted from (V-19), the 


open-loop characteristic is that of a perfect integrating servo: 
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The closed-loop behavior is then 
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For purposes of illustration,the values of the time constant /( 8, /k,)= (A Kk) may 
be considered as extending from ,05 to .025, with the latter value running into 
a range where it is sometimes difficult to stabilize the system, 


SECTION 3 = THE EFFECT OF INPUT COUPLING 

The stabilizing means considered in Section 2 are rather limited in effec~ 
tiveness. If the hydraulic control system must be very fast, more effective 
means of stabilization are desirable. The additional means considered in this 
volume use the combination of a force and a relatively weak spring to secure 
compensating displacements at some point of the system. There are, of course, 
numerous variations possible in physical application of this idea, The material 
presented in Section 1 facilitates discussion of some of these possibilities, 

Consider now a case where the valve damping &. or the input coupling spring 
k 


€ 
a system with a high gain valve. The block diagram in Figure V-3 then applies, 


» or both,and also the input damping &- are augmented by physical changes in 


With £ set equal to zero, the operator ¥%, is written from (V-25): 
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The frequency a), of the typically non-minimum phase quadratic in the numer~ 
ator drops sharply with an increase in valve damping 8, 3 however, the decrease 
is offset by a high gaink, /8 . This quadratic can fall inside the denominator 
frequency @), , but it 15 improbable that it should ever reach the neighborhood 
of a), and a), e With: Z; characteristically a small negative danping, there is 
little phase lag at important frequencies, The quantity 7, can then be assumed 


to have the same form as (V=27) of Section 2: 
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The spring constant A, usually expresses the flexibility of a long, light 
cable, and is generally fairly small. The input mass M; is kept low by weight 
considerations, and is normally of about the same magnitude as Mp » the piston 
mass. If the product k, M,; is sufficiently small, the quadratice of (V=38) can 
be split into first order factors by increasing 8 and &,. In this case, an 
equalizer might take the form: 


(V=39) %eg 
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Equation (V-39) suggests several possibilities, which will be briefly 
sketched in the following paragraphs. 
If B; and A, are small, and 6,, /k; is large, 
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In this case, valve damping combines with a weak input spring to form a lag 
factor in the open-loop response, 
If k,. is small, and &, \ ke is large, and & is fairly large, 
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The modifying factor in this equation is a lag-lead. The lead would probably be 
too high to affect stability, but the introduction of 8 makes a lag at low 
frequency more easily realizable. 


If 8, is small and 8; /k; and k, are large, 
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A lead-lag modifier can be realized if the gain of the valve is increased to 
compensate for the decrease in gain in the equalizer. 

The spring constant A; is to some extent a fixed factor in design; how- 
ever, it varies considerably between control systems in different aircraft and 
also between the several control channels in the same aircraft. If modifying 
devices of the types discussod above are used for stabilization, they must be 
"tailor-made" for the specific installation. This is not too serious a drawback, 
because aircraft control systems in general do not permit m:ch standardization, 

The effect of a lag stabilizer on an unstable system is shoum in Figure 


V-6, and the effect of the lead=lag equalizer is shown in Figure V-7. 
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Figure V-6, The Effceet of Valve Damping on Fully-Povwered System 
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CHAPTER VI 
THE POWER BOOST HYDRAULIC CONTROL SYSTEM 


The material of this chapter is a continuation of Chapter IV. It is prima- 
rily concerned with the evaluation andthe approximate factorization for the 
power boost case of the expressions occurring in Equations (IV-12) and (IV-13), 


which are repeated here for convenience: 


(VI-1) x= ZEr+ YF 


(VI-2) e&-= Yeg [%. (ax:)- x, + %-F] 


The corresponding block diagram is shown in Figure IV-2, 

In Chapter V, it has been shown that the effects of the force F are always 
stabilizing and usually emall. in the fully-powered case if the behavior of the 
control surface can be effectively expressed by use of the single node Z,. 
Moreover, it is necessary to design the control system to insure stable opera- 
tion at low airspeeds or on the ground. An analysis of the fully-powered case 
with F£ taken as zero thus yields conservative conclusions as to stability if 
coupling of multiple modes of vibration of control surface and parent surface . 
is negligible. The situation in the power boost case is in general somewhat 
different. As will be shown later in the chapter, it is possible, by a suit- 
able choice of the power boost ratio Ib ,» to adjust frequencies of the open- 
loop transfer function leg % so as to eliminate a stability problem arising 
in ground operation of the power boost system. A consequence of the pronounced 
effect of the power boost ratio on stability is that the aerodynamic force f 
may be either stabilizing or destabilizing, depending pn the value of 6, If 
4 is small, however, the conclusions of Chapter V separding the effect of aero- 


dynamic force on stability are equally valid in the power boost case, and in the 


VI-1 


following analysis it is assumed that the quantity f is zero. 


(VI-3) x, = EE 


(VI-4) E= keg [%, (axe) - x. 


These relationships are illustrated schematically in Figure VI-1l. 





Figure VI-l1. The Power Boost Hydraulic Control System 
with Negligible Airload 


The quantities y , FY , and Y, of Figure VI-1 are obtained from the 


eg < 
determinants associated with (IV-12) and (IV-13). Writing of these results is 


simplified by assuming relative magnitudes of parameters, as was done in Chapter 


¥e 
OB, >> & Ms >> M 
B >> (=), Ms >> Mp 
a-b 
ke >? Re 
VI-2 


i A Se ee ee eee 


re rae 
4 


a 


i 


r= rr" 


wg iow 


ow 


ony 


ew 


a ee a a er ae ee ee ea 


does 


M, M 

he” he >a 
Mi,» M, 

ae SS 2 = 
ky, kp M;, ?? b M, 


kp M;, M, >> 8, &;p Mp 


The numerator of )% takes the form: 
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In a typical system the piston mass Mp is small, and the piston spring Ap 
is several times larger than either the oil spring A, or the cylinder-surface 
coupling spring 4,. Numerator and denominatcr of ¥, then each have a quad- 
ratic factor at the approximate high frequency vkp Mp’ « With these quadratics 
effectively cancelling out, a good approximate factorization of ¥ ie: 
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It should be noted that the expression (VI-7) for ¥, is identical in form 
with Lhe expression for the fully-powered case, (V~25), with f=O. If the power 
boost ratio //b is less than infinite, there is a slirht increase in the gain 
«Jj, and a reduction of the frequencies uw), , @, , and w,. From an examination, 
un the basis of typical parameters, of the frequencies @, and @, , it can be 
stated that the quadratic containing M, can be neglected because it lies at 
very high frequency; if the power boost ratio is high, i.e., if the system 
approaches a fully-powered system, the quadratic at Wy is out of the range of 


interest. The transfer function % may then be rewritten in the fairly simple 


form: 
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where the frequencies and the middle terms of the quadratics are as given in 


(VI-7). 


The effects of the small masses M. and Mz on the quantities eg and Ja, 


can be entirely neglected, and these transfer functions become: 
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The latter quantity is, of course, unimportant in a determination of stability. 
The open-loop transfer function of the positional servo is the product 


Yeg » By examining (VI-8) and (VI-9), it is seen that the control surface 


quadratic at ta), cancels in the product. It would be desirable in order to 


insure stability that the numerator of %, be such as to cancel the denomina-~ 


o 
tor quadratic of  . The numerator of %, can de ficbored to obtain this 


cancellation if the following conditions sre inet: 


kh; + abk, >» ke 
M; B, 
—s () 
(VI-11) 
B. —> 0 


k; (k, + kp) 


bla-b) * ray? 
o Mp 


The determination of 6 in accordance with the last condition may be in 
conflict with other considerations of control system desimn. For a particular 


aircraft, the force to which the pilot or autopilot is thus subjected may be 
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too large or too small for proper operation. 
If the conditions of (VI-11) are satisfied, the block diagram of the servo 
is as shown in Figure VI~2. 





Z, = Mye"+ Bys +(ky + abk,) 


z= Wye" +[By-+ a(a-b)By|s +(Kea + a(a-b)ke| 


Figure VI-2. Power Boost System with Power Boost Ratio 
Adjusted to Stabilize the System 


The valve damping 4. and the coupling spring k, were shown in Chapter V 
to be of value in stabilizing the fully-powered system. In the power boost 
case these effects are masked by the force feedback which arises because 4 is 
greater than zero. It is desirable, then, to keep these parameters as small 
as possible in the power boost system. 

The discussions of the fully-powered system in Chapter V and the power 
boost system in Chapter VI have necessarily been brief, but it is hoped that 
they will have provided the reader with some insight into the behavior of 
hydraulic control systems and into the problems of reconciling a rapid response 


with stability. 
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| CHAPTER VII 
SPECIAL CONSIDERATIONS IN HYDRAULIC CONTROL SYSTEM DESIGN AND ANALYSIS 
SECTION 1 - INTRODUCTION 


The purpose of this chapter is to consider several special topics in the 
design and analysis of hydraulic control systems. The material presented here 
is to be regarded as supplementary to the discussion and analysis of Chapters 
III through VI. 

Section 2 is concerned with forces acting between valve slide and valve 
housing due to inertia of the hydraulic fluid. It is show that the effects of 
such forces are well represented by adjusted values of parameters already con~ 
sidered in Chapter III. 

Section 3 discusses the hydraulic actuator with a three-way valve and how 
it differs from the more conventional four-way arrang: 1ent. The analysis of 
previous chapters is applicable with only small alterations. 

Sections 4 and 5 deal briefly with additional means for stabilizing a 
hydraulic control system. The results are comparable with those of the methods 
discussed in Section 3 vf Chapter V and in Chapter VI. Gain adjustment through 
the use of an auxiliary linkage is the subject of Section 6. In Section 7 the 


recently developed electrically operated transfer valves are described. 


SECTION 2 = LONGITUDINAL FORCES ON CONTROL VALVE SPOOLS* 
In the design of the hydraulic control. valve for a high performance servo 
system, some knowledge of the force required to operate the valve spool is 


necessary. In addition to the inertial and frictional forces which can to a 


*This section is largely an adaptation of Shih-Ying Lee and John T. Black- 
burn, "Axial Forces on Control-Valve Pistons," Meteor Report No. 65, Dynamic 
Analysis and Control Laboratory, Massachusetts Institute of Technology, Can- 
bridge, Massachusetts, June 1950. 
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degree be calculated or which can be estimated on the basis of previous design 
experience, two other axial forces which are functions of flow through the valve 
exist. The more easily recognizable of these two forces is the so-called 
Bernoulli (or centering) force. This force is a nearly linear function of the 
volumetric flow rate through the valve. Hence for a linear flow valve, the force 
is proportional to valve spool displacement. The force displacement gradient 


is thus similar to an equivalent spring constant. The remaining force is pro- 


portional to the rate of change of flow rate through the valve. This is equiva~ 


lent then to an additional valve damping term. The total forces involved are 
those accounted for earlier in Chapter III, specifically the valve damping 4,- 
and the valve spring constant ke . In general, each of these two forces can 
be either positive or negative. The Bernoulli force in the usual valve design 
is positive; i.e., it resists motion of the valve spool in a direction which 
increases the flow rate through the valve. This gives rise to the misnomer 
"centering" force. 

Because these two forces are functions of the geometry of the valve spool 
and sleeve configuration, the designer has some control over these forces. A 
well-designed valve will be of such a configuration that the Bernoulli force is 
minimized and the spool damping effect is positive. 

Three main assumptions are made at this point: 

1. The fluid is considered to be incompressible and non-viscous, 

2. The peripheral width of the orifice is assumed to be large in com 
parison with its length so that the flow can be considered two- 
dimensional. 

3. The flow is considered to be irretational in the adjacent region 


upstream of the orifice, 
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degree be calculated or which can be estimated on the basis of previous design 


experience, two other axial forces which are functions of flow through the valve 


exist. The more easily recognizable of these two forces is the so-called 


Bernoulli (or centering) force. This force is a nearly linear function of the 


volumetric flow rate through the valve. Hence for a linear flow valve, the force 


is proportional to valve spool displacement. The force displacement gradient 


is thus similar to an equivalent spring constant. The remaining force is pro- 


portional to the rate of change of flow rate through the valve. This is equiva~ 


lent then to an additional valve damping term. The total forces involved are 
those accounted for earlier in Chapter III, specifically the valve damping 4, 
and the valve spring constant ks « In general, each of these two forces can 
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is positive; i.e., it resists motion of the valve spool in a direction which 
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dimensional. 
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For small valve spool displacements, the first two assumptions are reason- 
ably accurate. The third assumption is not true throughout most of the hydraulic 
circuit but is fairly accurate in the small region of interest. 

Consider the figure below which shows an element of a control volume made 


up of a square land valve spool and the sleeve in which it operates. 







Sleeve 
(Stationary) 





Figure VII-1. An Element of a Control Volume 


The force existing on the valve spool is the axial component of the net 
rate of efflux of momentum through the control volume boundary a-b-c-d-e-f-g-h- 
idea. The area a-b of the vena contracta of the jet is much smaller than the 
area d-e where the fluid enters the upstream chamber A. Since the velocities 
are inversely proportional to the respective areas involved, the influx of 
momentum through d-e is considered to be negligibly small compared to the efflux 


of momentum at a-b. Therefore, the net axial force is 
(VII-1) f-g = QU0 cosO 


where 
Q is the total rate of flow (in? sec.~1) 
U is the velocity at the vena contracta (in. sec.~+) 


fo is the fluid density (1b. in.~4 sec?) 
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For the configuration shown with sharp edged orifices and very small radial 
clearance, the angle O is 69 degrees. This value has been determined analyti- 
cally and verified experimentally. Since QO is always less than 90° for the con- 
figuration shown, the force f_, of (VII-1) is always positive and tends to center 
the valve spool. This is true regardless of the direction of flow. In either 
case it is the configuration of the downstream chamber that determines the force, 

Because it is the downstream chamber that determines the force, there exists 
the possibility of shaping the downstream chamber of a port in such a way as to 
produce a negative (decentering) force. If such a configuration is used for one 
port of a four-way valve, the negative force can be made to cancel the effect of 
the positive force existing at the other port. 

In the negative force port configuration, the axial component of the efflux 
of momentum from the downstream chamber is made greater than the influx. There- 
fore a negative or decentering force at one of two ports is generated. The con- 


figuration of such a port is shown in the figure below, 


| 






Sleeve 


Figure VII-2. A Negative Ferce Port Configuration 
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In Figure VII-2, the chamber in the spool is shaped somewhat like a turbine 
bucket with the sleeve relieved to form an extension of the downstream chamber. 


The force existing on the spool is 
(VII-2) E-5 = QUe(cos 8 - cos @) 


This force is negative since cos Q is considerably larger than cos ©, . 
In addition, the eddy re-entering the chamber D at the angle Q further increases 


the negative force. A four-way valve with force compensation is shown below. 
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Figure VII-3. A Four-Way Valve With Force Compensation 


The total force on the valve spool (exclusive of inertial and frictional 
effects) is developed in the referenced report (see footnote, page VII-1). The 
total axial force is 
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and 


dQ 2ZAp* dx’ 
dt —@*V 5 dt 


Ap is the pressure drop across the orifice 

x” is the linear opening of the orifice 

ay is the peripheral width of the orifice 

Cz is the orifice discharge coefficient (approximately .65) 
C. is the velocity coefficient (approximately .98). 


Therefore, 


(mist) = ~ (26, C.lcor)urap]x’* [Gur (2Apa" &] SE 


where £= %,-X, of the figure below. 


x 





Figure VII-4. Flow Into the Annular Chamber of the Valve Spool 


For the direction of flow as shown with the jet discharging into the annular 
chamber of the valve spool, /™ X,~- %, and is positive. This results in a nega~ 
tive damping coefficient in the force equation. If the flow is reversed and is 
from the annular chamber out through the orifice, the damping is positive. Ina 
multi-orifice valve, the positive damping must be greater than the negative damp- 


ing to prevent valve instability. 
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SECTION 3 = ACTUATORS WITH THREE-WAY VALVES 
(a) COMPARISON OF ACTUATORS WITH THREE-WAY AND WITH FOUR-WAY VALVES 
In the discussion so far concerning control valves and actuators, only the 
symmetrical four-way valve and balanced area actuator have been considered. The 


typical configuration is repeated below. 
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Figure VII-5. A Four-Way Valve and Balanced Area Actuator 


The operation of this device has been described in more detail elsewhere 
in this volume. In general, there are two equal pressure drops at the two valve 
ports, and assuming an incompressible fluid, there are two equal volumetric flow 


rates into and out of the respective cylinder volumes, 
(VII-5) AP =PB-R =f-R =AR 
and the load induced pressure 


(viI-6-) R= P-R 
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and 
dQ 2Ap’ dx 
at ~ e“Vo 5 at 


Ap is the pressure drop across the orifice 

x is the linear opening of the orifice 

ay is the peripheral width of the orifice 

Ce is the orifice discharge coefficient (approximately .65) 
C, is the velocity coefficient (approximately .98). 


Therefore, 


(VII~d,) Ro=un [2Cg Ca (cos @)urAp |x’+ Ge (2hp0' L| ax 


where £= %,-x, of the figure below. 
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Figure VII-4. Flow Into the Annular Chamber of the Valve Spool 


For the direction of flow as shown with the jet discharging into the annular 
chamber of the valve spool, _™ x,~-%, and is positive. This results in a nega- 
tive damping coefficient in the force equation. If the flow is reversed and is 
from the annular chamber out through the orifice, the damping is positive. Ina 
multi-orifice valve, the positive damping must be greater than the negative damp- 


ing to prevent valve instability. 
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SECTION 3 ~ ACTUATORS WITH THREE-WAY VALVES 
(a) COMPARISON OF ACTUATORS WITH THREE-~WAY AND WITH FOUR-WAY VALVES 
In the discussion so far concerning control valves and actuators, only the 
Symmetrical four-way valve and balanced area actuator have been considered. The 


typical configuration is repeated below. 
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Figure VII-5. A Four-Way Valve and Balanced Area Actuator 


The operation of this device has been described in more detail elsewhere 
in this volume. In general, there are two equal pressure drops at the two valve 
ports, and assuming an incompressible fluid, there are two equal volumetric flow 


rates into and out of the respective cylinder volumes, 


(ViI-5) AP =h-R = R-& = 4B 
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This can be illustrated by a bridge circuit shown in Figure VII-6. 





Figure VII-6. Bridge Circuit Analogous to Hydraulic Servo of Figure VII-5. 


At this point, it appears desirable to describe another type of valve and 
actuator that is somewhat simpler in operation. Because the valve is simpler in 
construction, the device becomes attractive as far as design is concerned. A 
typical configuration of the three-way valve and differential area actuator is 


shovm in Figure VII~-7. 
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Figure VII-7. A Three-Way Valve and Differential Area Actuator 
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Here the actuator piston is constructed so that the areas of the two sides 
of the piston head, subjected to pressure, have a ratio of two to one. To com-= 
pare equivalent systems, suppose that the smaller area is equal to the area 
presented by either side of an equivalent balanced area piston. This smaller 
area chamber communicates directly to the high pressure supply (F, ) « For the 
actuator to be in a state of static equilibrium under no load conditions, the 
pressure ( i ) in the larger area chamber must be maintained at a value of 
one-half the high pressure supply. Therefore, control of the actuator resolves 
itself into controlling the pressure and flow into (or out of) one chamber of 
the actuator. Because of this, only one controllable orifice within the valve 
is required, 

The actuator is operated by controlling the magnitude of a variable force 
against which a constant force operates. For the configuration shown, the same 
maximum force in either direction is available and is equal to the respective 
maximum force attainable in either direction by an equivalent balanced area actu- 
ator. 

Both types require the same amount of supply volume per cycle. However, 
the differential area actuator delivers twice the volume to the return line over 
a half cycle than does the balanced area actuator. Over the other half cycle, 
the differential area actuator delivers zero volume to the return line. In 
addition, the piston diameter of the differential area actuator is larger than 
that of an equivalent balanced area actuator. 

In concluding the discussion of the two types of valves and actuators, it 
should be pointed out that at the mid-point of the stroke, the hydraulic compli- 
ances of the two types are equal. However, in the case of the balanced area 
actuator, controlled by a four-way valve, the effective spring constant of the 


oil within the cylinder increases as the actuator departs from the mid-stroke 
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position. In the case of the differential area actuator controlled by a three= 
way valve, the compliance is proportional only to the volume of oil trapped in 
the larger area chamber. Therefore, the hydraulic compliance of such a device 
is minimum at one extreme of actuator stroke and maxinum at the other. Also, 
in the case of the three-way valve actuator, since the small area side of the 
actuator is directly connected to system supply pressure, pronounced variations 
in system pressure are more apt to affect the actuator output position. Actu- 
ator output motions under load perturbations are also proportional to the actu- 
ator hydraulic compliance, 

In the discussion, a two-to-one area relationship was assumed in the case 
of the differential area actuator. Because it is convenient in the design of 
an actuator to incorporate O-rings for sealing means, it may not always be fea= 
sible to provide a precise two-to-one area relationship. Under such conditions 
the actuator response will not be truly symmetrical about a no-load operating 
point. In general, if the actual area relationship is within 10% of the ideal- 
ized two-to-one relationship, it becomes difficult to detect, any asymmetry in 
the actuator output. 

Similar comments apply in the case of the balanced area actuator where the 
areas at both sides of the piston head are not precisely equal. A typical case 
is where it is not desirable to balance the areas with the actuator blind shaft 
extension. 

A lack of area balance in either type actuator causes a slight. departure 
from symmetrical response, The operation is similar to that which would be 
caused by a moderate load bias, which is a fairly common operating condition 


for any type actuator, 
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(b) ANALYSIS OF ACTUATORS WITH THREE-WAY CONTROL VALVES 

Since the piston areas of this actuator are unbalanced, it might appear 
that analysis is more difficult; however, this is not the case. For this type 
of actuator, the flow is metered across two lands of the valve rather than across 
four lands as in the four-way transfer valve. If the port connecting the supply 
pressure to the small area side of the piston is large in comparison with the 
metering ports of the valve, and the supply pressure is assumed constant, it 
can be seen that the analysis of flow equations need only be concerned with the 
flow in and out of the large area side of the piston, since the small area side 
of the piston can be assumed to do nothing more than supply a constant force. 
If the large area is twice the small area, the expression for the load induced 
pressure / will be: 
(viI-7) R= -£ 
For purposes of computing the oil spring the volume of oil under compression is 
only the oil on the large side of the piston; therefore, the analysis should 
always be carried out for a position of the cylinder where the large area side 
of the piston has its greatest volume, since this will be the point where the 
oil spring is the weakest. This is the extended position for the actuator 
shown in Figure VII-7. With the small changes mentioned above, analysis is 


then carried out by the methods of Chapter 3. 


SECTION 4 ~ EQUALIZATION BY HYDRAULIC MEANS 
Chapter III of this volume describes a typical valve-controlled servo 
element. If the piston rod spring constant Ap is very large and the piston 


mass Mp is very small, (III~44) can be approximated by: 


| A 
(VII-8) X, = A -(Bs +l) (42) 
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Assume that the complete set of performance specifications has been made 
available and that the actuator size has been selected. Therefore, for a prede- 
termined value of system supply pressure ir , the piston area A has been 
established, and the effective volume of fluid under compression vis approxi- 
mately determined. The servo designer has no control over the load induced 
pressure 4 or the fluid bulk modulus VV. 

The only design parameters over which the servo designer has any effective 
control are the slopes C, and Cp of the valve flow curves. Further, because 
of the difficulty of evaluating Cp, the servo designer has only one design 
parameter by which he may optimize his design. This design parameter is Co 
which is essentially a gain term. 

Suppose then that the designer has selected a value of C, which seems to 
be satisfactory. After fabrication of a prototype, certain indications of incip~ 
dent instability become apparent. The designer then can reduce the value of C,. 
However, with the servo gain diminished, the system may not be able to satisfy 
the design requirements. Another solution is to provide such devices as heavy 
valve damping or a bleed across the cylinder parts. However these devices are 
also performance inhibitors and are generally unsatisfactory. 

The problem now resolves itself into one of equalization, of providing 
methods of utilizing elements with essentially fixed characteristics in order 
to achieve a type of behavior that is in no way inherent in the valve controlled 
servo element itself. 

An equalization means which is widely used for electrical networks is the 
RC lead network. This network,which is shown in Figure VII-8,is put in series 
with the control element and effectively modifies the signal to the control ele- 
ment (in this case the servo valve) in a particular manner. The result is that 


the servo can operate at a higher gain in the forward loop (with a higher value 


VII-12 


ieee eerie 


ca PF Fe Re ee eee 


ae ae ee ee ee 


- 


J 


™™S — 


Ee EEO, lol 
JH ES i aS 


WN = 


— 


eee 


Section 4 


of C,) without resulting instability. This is due to the positive phase 


shifting characteristics of the network. 


The basic lead network is shown below 
C 





Figure VII-8. Basic Lead Network 


The relationship of the voltage out @, to the voltage in @, is 


E, (oc Ts + /) 
TT. ol * Re beat In Ss A 
(VII-9) 3 o(Ts*/) 
where 
Oc Rt fe 
Re 


ew eC 
R+R, 
Such a network is of little interest to the hydraulic servo designer unless 
of course the control valve happens to be positioned by some electromagnetic 
device. The remaining portion of this section will discuss means by which hydrau- 


lic and mechanical elements may be used to provide an equivalent effect. 
One possible way of achieving the desired equalization is shown in Figure 


VII-9. 
In Figure VII-9 the main valve and actuator are exactly the same as those 
described in Chapter V. The auxiliary valve, the auxiliary piston, and the float~ 


ing link are added. All hydraulic components are located in a common housing. 
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Figure VII-9. Typical Hydraulic Actuator with Equalization 


Therefore the input signal x; and the actuator output position x, are 
measured relative to a fixed reference. The displacements ~, of the auxiliary 
piston, ~, of the auxiliary valve spool, and € of the main valve spool are 
measured relative to the movable actuator housing. 

First consider the system at rest as shown in the figure. Suppose a step 
input x, displaces to the right point "@ on the linkage. Initially then, the 
link rotates in a clockwise direction about point "d° as an instantaneous center. 
This results in a displacement x, of the auxiliary valve spool and a proportional 
displacement € of the main valve spool. Both of these valve displacements are 


towards the lett. The results are that the auxiliary piston moves to the right 
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relative to the movable housing and the main actuator moves to the right. As 
the main actuator moves to reduce the error to zero, the auxiliary valve spool, 
because of the linkage geometry, crosses its neutral point before the main valve 
spool reaches its respective neutral position. Thus the auxiliary piston, the 
auxiliary valve spool, and the main valve spool all simultaneously are returned 
to their respective neutral positions as the system error reduces to sero. 

The figure below shows the link in two positions. Figure VII-10(a) shows 
the angular position of the link immediately following the introduction of a 
step displacement of X; , before the actuator has changed position in response 
to the signal. Figure VII-10(b) shows the link as the actuator approaches syn- 
chronism with the signal. Note that in (b) the auxiliary valve spool has crossed 
its neutral position and that the auxiliary piston has reversed its direction of 
motion, returning to its neutral position. The heavy arrows show the direction 


of motion and velocity of the points indicated, 





Figure VII-10. Link Displacement 
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From the linkage geometry it is possible to write two equations for the 
link instantaneous angular position O: 


ei TOT ig, hat = Ke 
(VII 10) 8 Z ZL, 
and 
Ap-~ X%e7 Xe Xam E 
(VII-11) G2 ———_ + = —t— 
also 


g 
(VII~12) = a =n 


where A, is the auxiliary system gain. Equation (VII-12) is valid if the aux- 
iliary actuator is unloaded; that is, if 4 and A, are very small. 


(VII-13) *,=-YE 


where 7 is the transfer function of Chapter V. 
The minus signs indicate that valve displacement to the left results in 
actuator motion to the right. Combining (VII-10) and (VII-12) 


(VII=-14) As (x, ~Xo~ Xe) = A, (x, t a) 

from which 

(VII-15) Xq = en. ae 
(4+ te) + 42] 


Equation (VII-11) may be rewritten as: 


(viI-16) (x, - x, 4, +4,) = (4,+ 4,4 4) my -4E 
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Substituting (VII-15) in (VII-16) and solving for E : 


A -3 ae 

(VII-17) Es - (xe- X,) = ceed enemas 

A,tte 4, eel 

4, +, )k, 
or 
€ aa %g (%2- %.) 

Equation (VII-17) may be combined with (VII-13) to form the same type of block 
diagram as that developed in Chapter V which is repeated here. 





Figure VII-11. Block Diagram of Actuator With Hydraulic Equaliser 


It is desirable that the equalization term leg be equivalent to that given 
in (VII-9). Therefore, 


(Let Ay) | + | 


teen ans (< Ts+/) -| A, Ak, 
«(Ts +/) A,+A, ee 
(2+ 2,)k, > 


For (VII-18) to be true, the linkage equation /,=4,+4, met be true; then 


Ltt A 
oo = foci, Riad 2h = ne 
(VII-19) 2, ae & (2+ Lyk, 


The corresponding quantities of the electrical equalizer are: 


(VII-20) Qc = Fit Re 
R, 
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~ ARC 
(VII-21) T= R+R, 

Other methods for providing equalization can be used. For instance, a 
floating valve spool which is not directly connected to the pilot's mechanical 
input motion may be employed. In this case the control valve spool is not posi- 
tioned by the system instantaneous error but rather by some particular function 
of this error. By selection of the particular function of the error by which 
the control valve is positioned, the type of equalization is determined. This 
is by no means limited to an analogy of the electrical R-C network selected in 
the foregoing. By the exercise of some ingenuity in the hydraulic curcuit design, 
a wide variety of equalization characteristics is made possible. 

In general, the proper application of equalization techniques will greatly 
improve the stability characteristics of a given servo system throughout a criti- 
cal frequency region. In many designs the normally critical design parameters 
such as valve gain, represented by the slope of the valve flow curve, becomes 
relatively less important. 

In addition, because most hydraulic equalization means employ an internal 
feedback loop around the contol valve itself, or a functional equivalent of such 
a feedback loop, the normal non-linear action of the valve becomes obscured. 
This also results in the fact that the device instantly responds within itself 
to the effects of load perturbations without error perturbations necessarily 


resulting. Further, the steady state load error can be made zero in most designs. 


SECTION 5 © EQUALIZATION BY STRUCTURAL FEEDBACK 
The principal problem in hydraulic control system design is to find means 


of de-emphasizing the incipient instability arising because of the large mass 
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of the control surface and the low rigidity of the oil in the cylinder. Force 
feedback, operating in conjunction with a fairly weak spring, is remarkably 
effective in this respect, as has been shown in Section 3 of Chapter V and in 
Chapter VI. However, the effectiveness of methods discussed to this point 
depends on the magnitude of flexibility in the input link (cable or pushrod), 
and the design must be tailored to a given control channel. The application of 
such methods often requires numerous time-consuming tests. Merely as a matter 
of speeding design, there is some practical nsed for a means of force feedback 
stabilization in which the dependence on input dynamics is not pronounced. In 
the following there is a brief discussion of a method which has acquired the 
name “structural feedback." # 


Consider the simplified system shown echematically in Figure VII-12. 





Figure VII-12. Fully-Powered System with Structural 
Spring and Multiplying Link 


The quantities B.,k, ,M,, 8: , & , and & are considered zero for purpose 


*# Rahn, J. J., and Kangas, E. W., Improvement of Power Surface Control 
Systems by Structural Deflection Compensation, BuAer Report No. AE~61-5~I1I¢ 
Report of the Second Piloted Aircraft Flight Control System S sium, June 2, 

P| L, and 5, 1952. | 
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of illustration. The pivot point of the input arm is in this case not considered 
fixed, but as being on a "weak bracket" with spring constant k, . At some point 
of the input arm with dimensionless radius a , motion due to both the immediate 
input %,, and the deflection x, of the pivot point is fed to a multiplying link. 
The valve error then becomes: 


(/-a) 


Cc 


(VII-18) € -(2 x) + 





It should be pointed out that derivation of this equation tacitly assumes 
that the input arm does not translate as a whole in response to load, as is 
approximately the case if the input arm is a quadrant with preloaded cable and 
the bracket spring A, is less than infinite. The immediate connection between 
the input system and the input arm of the actuator is now taken as a pushrod, 
in which case the equation (VII-18) applies. 

With no aerodynamic force Ff , the cylinder displacement x, can be expressed 


as a function of the valve error €, 


(VII~-19) m= EE 


°o 


where, for the simplified case shown in Figure VII-.~, 
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A force balance on the cylinder yields the equation 
(VII-21) kK, x, + k, (x,- x,) = O 


The valve error is obtained as a function of x, and ax,/¢ by combining 
(VII-18), the equation about the surface node (IV-1), and the equation express~ 
ing the force balance on the actuator (VII-21). 


(VII-22) E= %, ! %e (224) - a 


where 


ei = 758 
oe 
w, 
k 

wo, = [22 








ie i (I-a) 
Ml? a 


The block diagram of the servo is shown in Figure VII-13, 
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Figure VII-13. Block Diagram of Fully—Powered System with 
Structural Spring and Multiplying Link 


From a comparison of Xe 


numerator quadratic of Neg. is the same as that of the denominator quadratic of 


and % it is seen that if the frequency of the 


% 4 the control surface dynamics cancel and the open loop transfer function 


becomes Kk, /6, S . The cancellation requires that 
/-@ 
hy + Ki —— 
(VII=23) ee ot on ae ee 
ky ke A, kg Kp Ap 


or simplifying, and assuming 4, >>A,, 


ee ee 


(VII-24,) = Kk 


The choice of a and c is restricted by the fact that the ratio a/c enters as 
an initial gain. 


If the multiplying link is not present, c is unity, and (VII-24) becomes: 
(VII-25) -a =— 


The utilization for stabilization of structural feedback without a multiplying 
link thus requires that the input %; be reversed. It is assumed here that the 


spring A, represents added flexibility in every case. 
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It is by no means easy to design a bracket for strength and flexibility at 
the same time. If structural feedback is to be used for purposes of stabili- 
zation, the design process may be expected to be largely cut and try. 


SECTION 6 = MECHANICAL GAIN ADJUSTMENT 
A common fault found with the integral valve-cylinder actuators used as 
illustrative examples earlier in this volume is the difficulty in changing gain, 
since the gain is chiefly the slope of the flow curve. Figure VII-li shows a 
simple actuator where the open-loop gain is dependent upon the mechanical advan- 
tage of the linkage operating the valve. 
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Figure VII-14. Simple Hydraulic Actuator 
The valve error for the actuator shown in Figure VII-14 is defined as: 
(VII-26) E€& = d(x; - x,) 


If (VII-26) is written in the form presented in Chapters V and VI, the form 
would be 
(VII-27) E = Yo (x;- x) 


where 
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Since the transfer function %, Jé = A remains unchanged, the block diagram for 
the system is shown in Figure VII-15. 





Figure VII-15. Block Diagram of a Simple Hydraulic Acutator 


It can be seen from Figure VII-15 that the factor d directly effects the open 
loop gain of the system. A change in gain by a factor of two can be easily real- 
ized by this means, but ratios larger than two may be difficult to obtain due to 
the length of the linkage arms involved. 

It can be seen from Figure VII-15 that the leading effect of valve friction 
on an input mechanism is reduced if @ is smaller than one. This may be desir-= 


able in case an autopilot actuator is used as an input device. 


SECTION 7 = ELECTRICALLY CONTROLLED TRANSFER VALVE 

A type of control valve that has become widely used in recent years with 
hydraulic actuators is the electrically controlled transfer valve, It has not 
as yet been used for primary surface actuators on piloted planes because of its 
relative unreliability compared with manual control valves. However, it is widely 
used on missiles and a few autopilot installations. It is desirable, therefore, 
to have an analysis which will cover this type of valve. 

The valve analysis presented in Chapter III assumes that two variables .? 
and E were sufficient to define the behavior of the valve. The same analysis 
could be extended if a transfer function could be found between valve spool error 
of the electric transfer valve and the electric signal supplied to the valve, 


“or instance the current. However, for most valves used today, the dynamics of 


VII~24, 


l 


r— . 


~-— 


? 


Section 7 


this transfer function occur at such high frequencies that the relation between 


valve error and current may be regarded as a constant. Figures VII-l6a, b, and 


c show the experimental flow curves for one electrically operated four-way trans- 


fer valve in use today, 
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Figure VII-16,. flow Curves for Typical Electric Transfer Valve 
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It can be seen that these curves compare very closely with the curves presented 
in Figures III-4 and III-6 of Chapter III for a manually controlled valve. 
Therefore, the analysis presented in Chapter III can be used if the current to 
the valve is substituted for valve error. The transfer function %, / ¢ may then 
be obtained. The rest of the servo loop,exclusive of equalization, consists of 
feedback potentiometer and amplifier as show in Figure VII-17 and can therefore 


be handled directly by standard servo techniques. 


e4 4 Valve 
“% Amp aa Actuator = 
i Load 


K 
Feedback 
Potentiometer 


Figure VII-17. Block Diagram for Positional Servo Using an 
Electrically Operated Valve 


One effect that has been neglected in this analysis is hysteresis. The 
normal, manually controlled valve has a negligible hysteresis loop; however, 
this may not be the case for the electrically controlled valve. If the hyster- 
esis loop is not negligible, the analysis is more complicated.* Also, if the 
dynamics of the valve fall within the frequency regime of interest, the trans- 
fer function x,/¢ will have to be broken up as previously mentioned and these 
dynamics included, 

Electrically operated valves could be built as a three-way valves however, 


analysis would still be similar. 


% See Chapter VI, "Non-Linearities," Methods of Analysis and Synthesis of 
Piloted Aircraft Flight Control Systems, BuAer Report No. AE-61-4 I,Northrop 
Aircraft, Inc., Hawthorne, California, 1952. 
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CHAPTER VIII 
COMPONENT DESIGN FACTORS 
SECTION 1 = INTRODUCTION 


It is the purpose of this chapter to summarize pertinent component design 
factors and to show how the design of the control surface actuator is integrated 
into the over-all design of the airplane, i.e., how the hydraulic actuator design 
influences and is influenced by the design of other airplane components. Most 
of the discussion is confined to fully-powered actuators incorporating substan- 
tially closed-center valves; however, the presentation includes mention of some 
of the problems associated with other types of systems. The discussion has been 
kept quite general, although at several points specific magnitudes or relative 
magnitudes of parameters are quoted as being typical and/or desirable. It should 
be pointed out that individual design problems may present difficulties which 
are not completely covered in the sections to follow. It is hoped that the 


chapter will at least provide a starting point. 


SECTION 2 = ACTUATOR SELECTION 

A primary consideration in the hydraulic system design is the type of actu- 
ator that will be used. The actuator may be a moving cylinder type with a sepa- 
rately mounted valve or a moving piston type. Most of this volume has been con- 
cerned with the moving cylinder type of actuator, with integrally mounted valve. 
This combination provides inherent servo follow-up, as well as minimization of 
line lags between valve and cylinder; it also results in the simplest operating 
linkage and consequently little backlash in the feedback path. Disadvantages 
of the combination are that the moving cylinder "sweeps out" more airframe struc- 
ture than does the fixed cylinder moving piston type, and the follow-up linkage 
is not as readily adjusted for necessary "fixes" in the later design phascs. 
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It may be desirable to employ two hydraulic actuators rather than a single 
actuator. This is particularly true on multi-engine airplanes where the two 
actuators may be furnished hydraulic power from separate pumps on different 
engines. The dual hydraulic system has obvious advantages from a vulnerability 
and reliability standpoint provided, of course, that the aircraft can still exe- 
cute simple maneuvers with a single actuator operative. Ancther advantage of 
the dual system is that the two actuators may be preloaded against one another, 
thereby eliminating dangerous backlash under otherwise critical surface flutter 
conditions. The preload is accomplished by displacing the valves so that when 
the control surface is faired, the actuator cylinders tend to move in opposite 
directions. 

Many modifications are possible in the actuator configurations. Examples 
are "gain" changers, hydraulic equalization, and complex feedback linkages. 
These complications of design should be used only when necessary. The prime 
effort must be to keep the design as simple as possible, and yet fulfil system 
performance requirements and stay within airframe structural and space limita- 


tions. 


SECTION 3 ~ SYSTEM PRESSURE SELECTION 

Another consideration which must be established early in the hydraulic 
actuator design is the system pressure. To minimize cylinder size and weight, 
it is usually desirable to use as high a system pressure as possible. The 
advantages thus gained must be balanced against the disadvantages of providing 
extra engine power to drive a heavier and more costly pump. A high system pres- 
sure also complicates the problem of leakage from hydraulic lines and couplings. 

At the present time most hydraulic components and hardware are qualified 


for a pressure of 3000 psi. A pressure of 5000 psi is considered the optimum 
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pressure for minimum size and weight of the over-all hydraulic system of a pre=- 
sent-day high speed aircraft, but procurement of complete sets of high-pressure 
components and hardware is not currently possible. Within the next few years 
higher qualification pressures will be realized as airplane speed and hinge 


moment. requirements increase. 


SECTION 4 ~ ACTUATOR-SURFACE GEOMETRY 

The values of certain basic parameters are available to the designer of 
hydraulic control systems prior to the establishment of the control system con=- 
figuration. These data are: 

1. Maximum hinge moment on surface due to air loads, 

2. Surface deflection limits. 

3. Surface rate or maximum velocity. 

Assuming a certain system pressure, it is easily shown that the total cylinder 
oil volume must remain constant, irrespective of the surface horn radius (or 
actuator stroke) selected. To obtain the same maximum hinge moment, the cylin~ 
der area must, of course, be varieds this results in an altered compressibility 
of the cylinder oil even though the volume remains fixed. In many cases the 
dominant flexibility of the combined actuator—surface system is due to the com 
pressible oil, yet the resultant natural frequency of the surface acting through 
this oil spring also remains unchanged for various horn radii when the same hinge 
moment and surface deflection are required. 

Lengthening the surface horn radius has the advantages of relaxing the posi- 
tioning accuracy requirements of a given hydraulic servo and of minimizing the 
effects of backlash at all bearing, or attach, points. A disadvantage is that 
the actuator stroke is correspondingly increased, and consequently the actuator 


assembly may require a large space between extremes of its travel within the 
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parent surface structure. If, further, the servo valve travel limits are kept 
fixed, in addition to hinge moment, surface deflection, and surface rate, the 
hydraulic actuator gain is also increased; that is, the first order time constant 
of the servo is decreased, giving rise to a stability problem. In most instances, 
space limitations severely limit actuator surface geometry, thus making many the- 


oretical considerations purely academic. 


SECTION 5 — VALVE DESIGN 
Most present-day fully-powered actuator designs include modified "closed 
center" servo valves. These valves are constructed in such a manner that a small 
amount of fluid ("neutral leakage") is allowed to flow from the pressure line to 
the return line through the valve ports when the valve is in a neutral, or cen=- 


tered, position as shown in the figure below. 


i Pressure 





Faia The 
— To Cylinder 


With this arrangement, even when the cylinder is not carrying external loads, 
the pressure in the fluid on either side of the piston is very high (approxi- 
mately one-half system pressure). Any air which is entrained within the fluid 
has negligible effect on the compressibility of the oil. The controlled leakage 
feature of the valve also produces an "open-center" servo operation over a very 
er range of valve travel; in this range the pressure differential across the 
cylinder piston is primarily a function of valve displacement. Thus any move= 
ment of the valve will result in a pressure differential across the piston. In 
most valve designs, full system pressure differential is available to the ie: 
inder when the valve has been displaced approximately 5% of its full travel. 


The complete system usually is desijmed in such a manner that full valve travel 
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corresponds to two degrees or less at the control surface. Since a pressure 








differential on the order of 2% of system pressure is normally required to over- 
come cylinder friction, it is evident that any "flat-spot" in the plot cylinder 
velocity as a function of valve displacement is about 0.002 degrees. Care should 
be taken that the neutral leakage flow is relatively small since there may be a 
number of similar valves operating on a single hydraulic power source. The sum 
of the neutral leakage flows may become excessively large for the case of wind 
milling engines, at which time the hydraulic pumps are operating well below their 
rated capacity. 

The slope of the flow curve, i.e. flow valve displacement, isythe mst 
important parameter of a closed center valve, since it, together with the cylin- 
der area, determines the gain of the hydraulic actuator, or the time constant of 
the servo system.* A hydraulic system time constant of 0.05 seconds is usually 
sufficient for present-day fighter aircraft. In the case of an open-center 
valve design, the slope of the flow curve does not necessarily dictate the sys= 
tem time constant. 

The design of the actual valve orifice depends a great deal upon the maxi- 
mum Oil flow required, which in turn is fixed by the required surface speed. On 
many airplanes of the past, the surface speed has been fixed by a rule of thumb: 


"full surface travel in one second." In general this rule leads to surface speeds 


—~ 
— 


~~ 


that are much higher than necessary. The flight condition which usually dictates 
the maximum surface speed is the landing condition, where the hinge moment is 
low but rapid surface corrections are necessary. It is worthwhile to carry out 
an investigation to determine the actual maximum surface velocity to maintain 
good control of the airplane. One approach is to assume a two-degree-of~free- 


dom airframe and solve for the time required to reach a given change in airplane 


# See Chapter V. 
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attitude as a function of surface velocity. The input is assumed to be a chopped 
ramp (Figure VIII-1); i.e., the surface has a constant speed to its limiting 
travel and then stays at that position. The time-velocity curve resulting from 


this approach will have the form shown in Figure VIII-2,. 





Surface Displacement 
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Figure VIII-l. Input to Surface 
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Figure VIII-2. Response of Surface 
It can be seen that beyond the dotted line of Figure VIII-2, increasing surface 


speed does very little toward increasing a response of the airplane. NACA reports 


have shown that satisfactory landings can be made with surface speeds as low as 
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10 degrees per second if the pilot is not aware that the surface is actually 
moving that slowly. Any reduction in surface velocity is reflected as a consid= 
erable weight saving due to the smaller pump required to drive the hydraulic 
system. 

In certain valve designs, the flow of fluid through the valve results in 
pressure drops which produce forces on the valve spool. These forces normally 
tend to return the valve to the neutral or centered position, Poorly designed 
valves can give rise to erratic valve centering forces which are objectionable 
to pilots and might cause system instability. For these reasons, attempts are 
made to minimize the effects by designing the valves in such a manner that the 
flow into and out of the valve is normal to the axis of the valve spool. The 
flow through the valve may be broken into a series of flow patterns where at 
maximum valve displacement the majority of the flow is in a region where dynamic 
centering forces cannot be applied to the valve spool. With this type of design, 
centering forces may be reduced to small values, 

Valve friction, even if rather small in magnitude, is very objectionable 
to a pilot from the standpoint of the character of the force. The coulomb fric~- 
tion component is constant in magnitude, but takes the sign of the spool veloc- 
ity. However, the coulomb effect upon the pilot is most pronounced during 
motion reversals, while, on the other hand, valve "stiction" (or break-away 
friction) appears at the stick only when motion is initiated from a rest posi- 
tion. As a result, both of these valve forces feel to the pilot as though the 
system has high inertia. Also,’a sluggish feel results if the viscous friction 
component of the total valve friction is too high. Maximum acceptable valve 
friction is approximately 2 lbs.; values less than this can only be obtained in 
sealed valve designs by reducing "0" ring squeeze, reducing valve spool diameters, 
etc. Where fluid leakage is no problem, simple lapping of the spool and sleeve 
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(or spool housing) produces relatively low valve friction forces. 

Adequate filters must be provided in as close proximity to the valve as 
possible. Further precautionary filters may be integrated in the valve body it- 
self. The possibility that foreign materials such as filings and chips might 
get into the fluid and jam the servo valve is a real problem. An additional 
safety measure which may be incorporated is a design of the valve and sleeve 
such that, rather than providing one large-area orifice, or slot, which might 
easily be jammed by a particle, a series of small diameter metering holes are 
uncovered in sequence by a sharp—edged valve spool land. The valve spool and 
Sleeve are hardened to Rockwell C-55 to C-58. As a result, large particles can- 
not enter, and the particles which can are small enough to be sheared off by 
forces developed on the spool by the pilot. Certain flight controller servos 
can also develop enough force to provide the necessary shearing action. 

It is desirable in valve design to provide a valve sleeve which floats on 
"O" rings within the valve housing. The radial clearance between the sleeve 
and the housing should be great enough that no binding of the sleeve will occur 
under extreme temperature conditions. The "O" rings also isolate the closely 
lapped sleeve and spool from external stresses which may exist due to mechani- 
cal mounting means. To further eliminate external forces on the sleeve, the 
valve spool travel stops should be on the housing rather than on the sleeve. 
The floating sleeve arrangement also allows easy servicing of the valve, since 
the sleeve and spool may be removed without disturbing the housing and its 


attached plumbing. 


SECTION 6 = ACTUATOR DESIGN 
Several spring constants which may affect system operation are inherent in 


the cylinder and piston design. ‘Two very important sources of flexibility are 
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the piston rod and the extension, or coupling spring, between the cylinder body 
and the surface horn. A reasonable criterion for these two spring constants is 
that they each be at least ten times greater than the spring constant of the oil 
within the cylinder. The axial flexibility of the cylinder body itself may be 
disregarded if a neutral leakage type valve is employed. The reason for this 
is that the sum of the pressures on each side of the piston remains approxi- 
mately constant, even with load variations, thus yielding the equivalent effect 
of a pre-loaded spring. It is common practice to insure long cylinder life by 
inserting a steel cylinder sleeve which is amaller than, but concentric with, 
the cylinder body; the sleeve is supported by appropriate mountings on both ends. 
The sleeve should not be so thin that its radial expansion under pressure will 
materially reduce the effective spring constant of the cylinder oil. The equiv- 
alent spring constant of this sleeve should be at least twenty times as large 
as the oil spring constant. In every case, the piston and cylinder assemblies 
must be designed for rigidity and not merely for strength. 

When the servo valve is mounted separately from the actuator, rigid tubing 


must be employed between valve and cylinder, the use of flexible hose in this 





application dangerously reduces the series combination spring constant of the 
oil and lines. The calculation of the total oi] under compression mst also 
include that in the interconnect lines for the remote valve configuration. 

An additional spring constant, which is avery important system parameter, 
"is the combined series spring constant of the surface horn, the torque tube 
between horn and surface (if such exists), and the inherent flexibility of the 
surface itself. This spring is dependent on the structural design of the sur= 
face and is not under the direct contro] of the hydraulic designer. The struc- 


tural designer should be aware that rigidity of the surface is an extremely 
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important parameter in the design of the hydraulic actuator system. The effec- 
tive spring constant of the surface is difficult to evaluate since the surface 
is a distributed parameter system and as such has many oscillatory modes. It 
is usually sufficient to evaluate the surface spring from the derived or meas— 
ured value of its fundamental torsional frequency as seen from the horn. The 
bending modes are usually relatively unimportant since most surfaces are sup~ 
ported at several points along the hinge line. ‘Parent surface bending if very 
great does not permit the simplification suggested above. Also, on recent air 
plane designs all-movable surfaces have been employed in which the bending modes 
are much more dominant than the torsional modes. 

There is one more design item over which the designer has limited control, 
This parameter is the effective damping coefficient between cylinder and piston 
due to O-ring installation. The O-rings are, of course, used to prevent piston 
rod leakage as well as load induced leakage past the piston from one side of 
the cylinder to the other. The introduction of the necessary O-rings into the 
design entails non-linear damping, the effective viscous equivalent of which is 
difficult to evaluate. Data on O-ring damping for various cylinder pressures 
and degrees of O-ring "squeeze" are contained in Air Force Technical Report 
No. 5997. The use of a single grooved piston sometimes results in O-ring roll 
which resembles backlash in the output motion of the cylinder (or piston) as 
the load pressure reverses sign. This O-ring roll may be minimized by providing 
two O-rings per piston and venting the peripheral area between the two to atmos- 
pheric pressure. In this fashion the force on either ring due to pressure in 


the cylinder.never reverses sign, 


SECTION 7 =- STABILITY ANALYSIS 
A brief review of a typical hydraulic actuator will now be presented to 


illustrate the effect upon system stability of the design parameters that have 
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just been discussed. The case which will be taken is the fully-powered actuator 
of Section 2 of Chapter V. 

In general, stability analyses of hydraulic control systems are made for a 
grounded condition of the aircraft, i.e., when there is no airload on the surface. 
This is usually the most adverse condition from the standpoint of stability since 
airloads entail additional surface damping and an additional restraining spring 
on the surface. It is a well-known fact that although many hydraulic control 
actuactors have been subject ‘e.enoud chatter or buzs, they have proved per= 
fectly stable in flight. Ground chatter also prevents checkout of the manual 
and automatic modes of the flight control system. 

With aerodynamic force on the control surface assumed negligible, the open- 
loop transfer function of the typical example is approximately given by (V-27), 


which is repeated here for convenience: 


2 
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The frequency response plot of this transfer function is shown in Figure VIII-3. 
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Figure VIII-3. Typical Open-Loop Frequency Response 


The most important parameters which place upper limits on the system gain 
are those which are contained in the damping and natural frequency terms of the 
denominator quadratic. It is seen that the frequency @, is dependent upon all 
the springs in the system as seen by the surface mass. It is important to keep 
this frequency as high as possible since it is relatively lightly damped; other- 
wise the system gain term, or zero gain Jine, must be adjusted at the expense 
of optimum performance. The preliminary analysis should indicate whether the 
design will be adequate, i.e., completely stable, yet with sufficient gain to 
yield a satisfactory time constant. 

The design criteria which have been set forth in previous sections of this 
chapter establish the cylinder oil spring and the inherent surface spring as the 
dominant springs in the system. The predetermined hinge moment, surface displac+-- 


ment, and surface rate requirements render the oil spring constant essentially 
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unalterable. Therefore, if the natural frequency of the system is too low, 
practically the only spring that cen be stiffened to the extent of having a 
measurable effect on the first denominator quadratic frequency «, is the sur- 
face spring. This requires coordination with the structures engineers to see 

if this is possible. Another problem which can arise is the possibility of flute 
ter-servo interaction, even though the hydraulic system is in itself completely 
stable. Certain types of flutter coupling modes can exist between the surface 
and the actuator; complete treatment of this subject is presented in Part II of 
this volume, 

Assuming that the analysis proves the hydraulic system to be still unstable 
in spite of alteration of the critical spring constants, various other correc- 
tive means are still available. Surface dampers may be added to reduce the 
height of the @, peak, but, of course, at the expense of reduced system per- 
formance and efficiency. Various types of feedback and equalization schemes are 
also possible with little adverse effect on performance, but resulting in con- 


siderably more design complication.* 


SECTION 8 - AIRFRAME~AUTOPILOT~SURFACE ACTUATOR CONSIDERATIONS 
if the hydraulic surface actuator is to be used in conjunction with an 
autopilot in addition to manual control, several additional problems arise. 
The most important of these are the time constant of the hydraulic system and 
the threshold. The time constant is not of primary importance on a manually 
controlled airplane since the pilot is capable of rate judgment and can make 


any necessary corrections within limits for excessive phase lag of the hydraulic 


_ system airframe combination. This means he can introduce a large amount of lead 


since the frequencies involved are usually low. The value for time constant 


previously given should be sufficient. However, the amount of lead that. can be 


* See Chapter VII, Sections 4, 5, and 6. 
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introduced through an autopilot is limiled; therefore, it becomes important that 
the combined phase curve of the autopilot, hydraulic system, and airframe allow 
the gain to be adjusted to give good control of the airplane. Since the airframe 
frequency response is dictated by the designed configuration, it is unalterable, 
and since the phase lead introduced by the autopilot is limited, the hydraulic 
actuator must be good enough to make the combined system satisfactory. There= 
fore, if an autopilot is to be used in conjunction with an airframe, it is essen- 
tial that the autopilot designers be consulted to determine what time constant 

is necessary for the hydraulic system. 

The second problem, low threshold, is important since the autopilot actu- 
ator will be required to make a great many small corrections because it cannot 
distinguish between phugoid and short period oscillations of the airframe. If 
the threshold is high, the autopilot hydraulic actuator airframe system will tend 
to have a limit cycle oscillation that the pilot will find objectionable. In 
conjunction with the low threshold requirement, it is necessary to maintain a 
satisfactory time constant for very small input amplitudes. Many hydraulic sys- 
tems tend to have a higher time constant for smaller amplitudes because suffi- 
cient care was not taken in the valve design. This effect will again tend to 
produce limit cycle oscillations in the combined system. 

Other hydraulic system parameters that affect autopilot design are loads 
imposed upon the autopilot actuator by the hydraulic system and velocity limit- 
ing of the hydraulic cylinder. For a fully—powered system, the only load imposed 
upon the autopilot actuator by the hydraulic system is the valve operating force. 
Assuming that Bernoulli forces have been designed to be small, the larger part 
of this force is coulomb friction. Since most autopilot actuators are electric 


motors whose torque error curve may have a low slope, this coulomb friction can 
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produce a flat spot in the autopilot actuator output. It canbe sccn that the 
valve operating force should be as small as possible. For a power boost system, 
the valve operating force is not nearly so important since it is a smaller part 
of the force required to actuate the surface. The velocity limiting of the 
hydraulic system is important only if the autopilot actuator is capable of bot- 
toming the valve, i.e., if it has a higher velocity. 

It can be seen from the foregoing discussion that the design of the hydraulic 
servo actuator involves information from the aerodynamics section, structures 
section, and autopilot section. It is important that the designer coordinate 
closely with all these groups to produce an integrated design. 
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PART IT 
HYDRAULIC ACTUATORS IN FLUTTER SYSTEMS 
CHAPTER IX 
THE INFLUENCE OF SERVOMECHANISMS ON THE FLUTTER OF SERVO CONTROLLED AIRCRAFT # 
SECTION 1 = INTRODUCTION {O SERVO-FLUTTER INTERACTION 


The extensive application of servomechanisms in the control of aircraft 

has been both cause and result of some of the most important recent advances in 
aircraft design. To a greater degree than was formerly the case it has become 
necessary to examine the controlled aircraft system as an integrated whole. 
In particular, the introduction in high speed aircraft of complex control sys- 
tems embracing several servomechanisms is making it a necessary part of design 
to study in some detail the dynamic interaction between the control system and 
the aeroelastic airplane. 

This chapter is concerned with the modification of classical flutter theory 
to account for the effect of a complex control system. Since control systems 
are highly variable in content and layout, a considerable portion of the chapter 
is devoted to control system elements which are likely to be present in practice. 
The material presented here is not to be regarded as a ready-made expression of 
any problem in servo~flutter interaction which may arise, but rather as a system- 
atic organization of the considerations which the design engineer must take into 
account in order to set down his problem. 

The central consideration in the study of dynamic interaction between the 
aeroelastic airplane and a complex control system including, in general, sev~ 
eral servomechanisms is that there are present in the control system energy 
sources which may either promote or retard a tendency toward flutter of the 
overall system. Most of present, day airplanes have powered controls, at least 
"4 Portions of the material contained in this chapter appeared previously 


in Air Force Technical Report #6287. Special permission for its use here has 
been granted by Wright Air Development Center. 
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in some of the control channels; so in a majority of cases there is at least 
one control system power source, besides the pilot, within the airframe. This 
power source is that of the surface actuator, which is most often a hydraulic 
servo actuator, but which may in some case be powered electrically or pneu- 
matically. Since the surface actuator has an energy source and is the direct 
coupling member between the control surface and the control system in the over~ 
all controlled aeroelastic airplane system, it occupies a position of primary 
importance in the present study. 

Surface actuators are in general not the only signal channels between the 
aeroelastic airplane and the control system in the over-all flutter system. 
Non=-rigid body motions of the airframe may, through structural deflections or 
sensor excitation, cause actuator inputs, which appear to the actuator or 
actuators as normal control inputs. The resultant control surface motions may 
promote or retard a tendency toward flutter of the over-all system. By far the 
most difficult portion of the modified flutter analysis is the determination, 
either experimentally or analytically, of the characteristics of elements of 
the additional signal channels. 

This chapter is not directly concerned with methods for solving the prob- 
lem of flutter of the controlled aeroelastic airplane system; it is rather 
devoted to an examination and phrasing of the problem. Available methods of | 
solution are described in Chapter X. 

The emphasis in this chapter and also in Chapter X is placed on the longi- 
tudinal control axis. However, all of the concepts and many of the developed 


details are directly applicable where the dominant control axis is lateral. 


SECTION 2 - MODIFICATION OF THE CLASSICAL FLUTTER SYSTEM 
(a) THE CLASSICAL FLUTTER SYSTEM 


A typical classical "flutter system" may be considered as a mechanical 
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system having three dogrees of freedom (torsion, bending, and surface rotation) 
and acted upon by two-dimensional aerodynamic flow forces. Each degree of — 
dom of the mechanical system is characterized by self~inertias, dampings and 
springs, and by coupling terms (via inertia, damping, or springs), with the other 
degrees of freedom. These self and coupling terms are all mechanical reaction 
forces, and the mathematical models describing their behavior (with zero applied 
aerodynamic forces) are linear differential equations with constant coefficients. 

The mathematical model of the applied aerodynamic forces cannot be expressed 
in the same way. Basically, the forces are considered as the effects of two= 
- dimensional flows and are usually described in terms of the rosults obtained by 
solving Laplace's equation for a particular two-dimensional configuration. The 
mathematical.models of the aerodynamic flow forces are therefore functions of 
the mechanical degrees of freedom and of time. To put the aerodynamic forces 
into a usable form (or even to solve Laplace's equation), it is necessary to 
assume’ a steady state oscillation of the three mechanical deyrees of freedom. 
Since the forces also depend upon the speed of the airflow about the particular 
two-dimensional configuration, the oscillation frequency is "non-dimensicnaliced," 
and referred toas the "reduced frequency," ber. On the basis of thic assump 
tion, the aerodynamic forces are expressed as a linear combination of functions 
of the three mechanical degrees of freedom but with the coefficicnts as functions 
of the reduced frequency. It is this complication, with the attendant special 
methods of solution, that has necessitated the creation of a speeial field for 
considering the dynamics of flutter. 

With this background, the mathematical model describing a typical classical 
flutter system is given in (IX-1), (IX--2), and (IX-3), where the constact eegtfie 


cicnt mechanical reaction forces are on the left. and the frequeney dependent 
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aerodynamic forces are on the right. The equations are written for the two 


dimensional case with zero structural damping and no tab. 


Surface Rotation Equation: 


(Ix-1) 


[bf 6A] +[f, +h(e-2) 5] & 45h = 7" 


Torsion Equation: 


(Ix-2) 


he + b(e-a) 58 + [4 K+, a] + Su 4A =M' 


Bending Equation: 


(Ix-3) 


> —— 
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Spf + 5nc & + (Mp + Ch) #h 


is the coordinate of axis of rotation 

is the instantaneous value of local torsional angular displacement of 
airfoil 

is the semi-chord of entire airfoil 

is the instantaneous value of control surface rotational displacement 
is the arbitrary damping coefficient 

is the torsional stiffness of parent surface about elastic axis per 
unit span 

is the torsional stiffness of control surface restraint about control 
surface hinge axis per unit span 

is the bending stiffness of parent surface about elastic axis per unit 
span 


is the instantaneous value of local bending displacement 
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is the moment of inertia of parent surface per unit span about, elastic 

axis 

I, is the moment of inertia about hinge line of control surface per unit 
length 

L is the total aerodynamic lift on entire airfoil per unit span 

M is the total oscillatory aerodynamic moment acting on parent surface 

about its quarter-chord point per unit span | 

is the total aerodynamic moment acting on entire airfoil about its 

elastic axis per unit span 

5. is the static moment of entire airfoil per uit span about elastic axis 

Ss isthe static moment per unit length of control surface about hinge 

line 

T” 48 the total aerodynamic torque acting on control surface per unit 

length about hinge line. 

These classical equations are complete in themselves if no further energy 
sources are present and if the description of the reaction forces is accurate, 
However, when a fully-powered or power-boosted surface control system is added, 
two different types of modifications to the classical equations are in general 
required. 

The first of these types of modifications deals with one degree of freedom 
vibration and is made necessary by the fact that the mechanical reaction forces 
described by the flutter equations no longer describe the physical situation. 

In this instance, the self-reaction forces in (IX-l), LA+GBA, representing 
effectively a spring (for example, torsion springs and a pushrod) and the surface 
inertia, require modification to describe the presence of the hydraulic actuator 


having no input valve motion. The nature of this entire modified self-reaction 


ope owen 
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force may be found by looking back into the surface actuator system. Since the 
first approximation is still an inertia and a spring term, with the inertia 
essentially /, , the entire system may be considered almost the same as the con- 
ventional system for sufficiently low frequencies. However, for high frequencies 
the spring term must be considered complex. This type of modification is con- 
Sidered later in detail. 

The second type of modification to the classical flutter system deals with 
systems of two or more degrees of freedom. In this case, one of the flutter sys- 
tem motions x, © , or hb, or a function of one of these motions, causes an input 
to the actuator and hence additional coupling of © with torsion or bending. As 
examples of this type of coupling, considers 

1. The case where rotation of the surface bends the fuselage, stretching 

one portion of a cable system more than another, resulting in valve 
motion which changes the surface rotation, etc., or 

2 The case where the torsional displacement, ©, causes a change in 

moment on the aircraft, in turn resulting in a change in load factor 
which is picked up by an accelerometer in a stability augmenter (or by 
a bobweight in a feel system), which then gives rise to a stability 
augmenter motion, a valve motion, and consequent surface rotation. 

The changes to the flutter system model in Case 2 above are obviously much 
more complicated, depending upon the characteristics of the stability augmenter 
and upon the characteristics of the feedback of a flutter system motion to the 
stability augmenter sensor. In the following portions of this chapter, stability 
augmenter characteristics are reviewed, and a method of handling this sort of 
feedback is presented. The types of feedback from a flutter system motion to 
various stability augmenter sensors are discussed, but nothing definitive can 


be said about exact forms. 
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In summary, the presence of the surface actuator makes it necessary to con- 
sider two types of modification to the flutter system: 

1. A fairly simple change in effective surface dynamics (with no input to 

the actuator). 

2. Much more complicated couplings due to the fact that the subsidiary 

systems cause surface actuator inputs. 

It should be stated at this point that a special instance of this second 
modification, which may be of great value in future designs, is the case where 
a function of a flutter system motion is picked up by a sensor, specifically 
installed for this purpose, which then moves the control surface in such &@ way 
that the influence of undesirable couplings is reduced; such a sensor is a flut- 
ter system stability augmenter. This sort of device could become quite impor= 
tant in future designs because it would make possible the use of certain air- 
plane configurations now considered unsatisfactory because of flutter. It 
would also be particularly helpful as a remedial measure in alleviating trouble 
caused by unexpected flutter and as a safety device for flutter testing. 

(b) MODIFICATION OF THE CLASSICAL SYSTEM 

The purpose of this subsection is to indicate the general form taken by 
the flutter equations when a surface actuator is added to the system. The md- 
ifying terms are expressed as transfer functions which are operators relating 
input and output quantities.* The procedure will be to assume that all possible 
feedbacks to the surface exist, and then to show how the general flutter equa- 
tions are modified. Later sections discuss in detail the actual forms of the 
general feedbacks assumed here. The reason for including this section at this 


a all 


# See: Methods of Analysis and Synthesis of Piloted Aircraft Flight 
Control Systems, BuAer Report No. AL-61-41, Northrop Aircraft, Inc., 


Hawthorne, California, 1952. 


IX-7 


Section 2 


point is to supply a better background for certain of the more detailed analyses 
given later. 


A general surface actuator has the defining equation: * 
(IX~4) M. = YAt Yx; 


where M, is the net moment due to the combined effect of (and its derivatives) 
and the effect of any input ~; to the servo device. The moment component 7.x, 
implies the existence of surface rotation through feedback coupling involving 

cc, , or fy which produce ~;. In terms of the flutter equations, 


(Ix-5) M, = T’~ [Ie + b(e-a)5,]% - Soh 


Note that in the classical flutter case the actuator system is not considered, 
and Y/F =I,8+C,8. Equation (IX-5) then becomes identical to (IX-1). 
Physically, %/O is the moment component due to surface rotation, and 
is the moment component due to an actuator input, 
If the input to the actuator, %,, is a function of <,, or + (due to 
the presence of coupling stability augmenters, autopilots, structural feedbacks 


directly in the actuator, etc.), the form of x, becomes 
(Ix-6) X= Yo + Eh + kA 


where Y is the total transfer function from the torsion angle ° , to actu- 
ator input “7, etc. 

In general, the X , % » % transfer functions are combinations of various 
feedback paths and contain portions which are functions of ba)/u as well as por- 
tions which are functions only of &@. Since the transfer function concept allows 


~~ % It is assumed that the surface actuator applies no forces that directly 
cause torsion or bending. 
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the analyst to effectively lump degrees of freedom, complex feedback paths can 
be handled readily in this fashion if the feedbacks are functions only of @ . 
When portions of the transfer functions become properties of bay , it is desir- 
able to consider certain of the degrees of freedom as additional ones to the 
three used thus far as examples, without lumping all the feedback degrees of 
freedom as transfer functions. This case is discussed in Section 3. In this 
section, the examples used are considered to be functions only of @ and are 
selected to illustrate the fact that feedback terms modify the conventional 
flutter terms. The equations used in this section would be general in nature 
if the Yle))’s are replaced by Y(bayv) "Ss . This 4s undesirable as far as actu 
ally solving the general flutter system problem is concerned. 

Combining (IX-4), (IX-5), and (IX~6), 


(K-7) YM + Ye + Er Ka) =T- [ht dc-alG,]& - 5h 


The modified surface rotation equation then becomes (for sinusoidal «<, G@ , and 
h)s 


(mr-8) Le) * Kd) (6 + {Le * 6E-a) SG ]Cw')+ % (d)% (w) Po 
+ {-a*So+ Fl ¥lW)}h = 7" 
In general, %(W) = {-J,@) + Gy + Ce (w)} 


(m9) EIy a) +6, + Gold) Ele E14 + {[0, + b(c-alSg|- a") 
+ GUNG} + Go + LN KPH = 7° 


The extent to which the flutter system is affected by any modifying system 


is measured by comparison of the values of the additive terms relative to the 
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corresponding classical terms (e.g., ¥(@)) (a) relative to -«) 5 )- In this 
instance, if ~a) 5, =e Fla) ¥ la) at the frequency of interest, the contribution 
of the additional bending coupling will be negligible. 

~ In a formal fashion, the entire problem of servo systems and flutter is 
described by (IX-2), (IX-3), and (IX-9) if the Y(.))3 are replaced by Y(bay7)’s. 
In the practical case, the forms of the additive quantities and detailed methods 
of analytically handling the additions are important. The sections that follow 
consider these aspects in detail. 

Section 3 discusses the forms of )} and 7; due to hydraulic surface actu- 
ating systems. A simple fully-powered case is taken as an example, While only 
hydraulic surface actuators are considered in detail, the use of electrical 
actuators leads to essentially the same results. Section 4 discusses the por- 
tions of © , % » and % due to stability augmenter systems. Section 5 dis- 
cusses the portions of Y , J; , and % due to cable systems and bobweights. 
Section 6 discusses the portions of 7, % , and % due to feedback of a flut- 
ter motion to a stability augmenter sensor, a bobweight, etc. With these "com 
ponent" characteristics established, the combination of individual characteris- 
tics established, the combination of individual characteristics into the total 


modifying termsG, ,7X% »% % » and YY is discussed in Section 7. 


SECTION 3 ~ SURFACE ACTUATOR SYSTEMS 
Hydraulic surface actuators are in common use in current piloted aircraft. 
In the case of missiles, the surface load and speed requirements and the missile 
airframe response requirements permit the use of a greater variety of surface 
actuator types, e.g., electrically powered units, pneumatic units, and electro- 
hydraulic units, as well as hydraulic actuators. The dynamical characteristics 


of all surface actuators are of essentially the same form, i.e., expressible by 
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similar generic transfer functions with interchange of certain basic system 
parameters. 

With the abové facts in mind, the system chosen for illustrative treatment 
is the fully~powered type formed by simplifying the system shown in Figure ITI-5. 
The power boost ratio Wb is taken as infinite, and the input mass M; » input 
damping 4 , valve damping 4, , and input coupling spring constant ke are each 
taken as zero. The system is then truly fully-powered; i.e., the input member 
is subjected to no force at any time. 
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Figure IX-l1. The Fully-Powered Hydraulic Servo Actuator 


If 2 is the surface horn radius in Figure IX-1] andd is the surface span 
served by the hydraulic actuator, there exists a simple relationship between 
the force £ of Chapter IV and the moment M, of (IX~5): 


(IX~10) M, = = F 
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The displacement x, and parameters 4, and 4, of Chapter IV are given in rota- 


tional terms: 


x, = Af 


d 
(IxX-11) Mz = Wis 


A, = 7G 
where and /, are the same as in (IX-1). The parameter 6, is a damping of 
control surface rotation which‘is seldom used in flutter analysis, since in this 
case it does not simlate any part of the aerodynamic force, and any remaining 
damping is usually negligible. 
Equations (IV-7), (IV-8), (III-58), and (III-61) are easily specialised to 
apply to the truly fully-powered system under discussions 
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The 73 given above are polynomials in § which reduce to unity if s equals 


SOro es 
Equations (IX~10), (IX-11), (IX-12), (IX-13) are combined to form: 


(IX-14) M, = 5a t ye 


where 
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The expansion of Ge 5G y is simplified by taking the flow 


roraL » and 
damping 4, as large in comparison with the cylinder damping 4. , as is always 
the case. In the typical case, the piston mass Mp is quite small, and the 
piston spring A, is several times as large as either the oil spring A, or the 
cylinder-surface coupling spring A, . The polynomials'g ,9 » 9s / k, , and 

9, then each have high frequency quadratic factors at the approximate frequency 
kp/Mp. These factors effectively cancel out of C, , er ae ay » and % ,» which 


may be written: 
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Frequency Response of Corrective Term -G 


Figure IX~2, 
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Figure IX-2 shows a generic frequency response of the corrective term ac ° 
This quantity acts to alter the effective rotational stiffness of the control sur- 
face from Cg=(4/d)k, , which would be the stiffness if the actuator cylinder 
did not move under load. The actual stiffness at any given frequency is the com 
plex quantity obtained by adding G, and Cy. It should be noted that the servo 
action is effective at zero frequency, in which eae a =-(L Hd x k: /k.+ k, ) 
The dynamics in Figure IX-2 are limited to a very low frequency range and a very 
high frequency range. For the simple case used as an example, thereis a broad 
intermediate region where there is no phase correction. This is the range of 
usual interest in flutter. However, the correction is not as simple with more 


complicated actuators. 
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Figure IX-3. Typical Frequency Response of Total Complex Spring Ge,,_,, 
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In Figure IX-3 it is asaumed that jw, and Wy, are practically equivalent; 
iee., (i/k,) =(I/k,) + (1/kp) « This is true if the oil spring A, is suffi- 
ciently high or the flow spring A, is sufficiently low. The figure graphically 
demonstrates the fact that when the two frequencies ay and &, approach coinci- 
dence (which is the case for 4,>>A. ), the complex spring becomes for all 
frequencies a pure spring, its value being very nearly that of the coupling 
spring Cz = (td) k, » Although the plot shows the value of the complex spring 
at zero frequency to be C, = (e'/d) A. , this is not precisely true. If will be 
noted that the complete analytical expression for the complex stiffness at zero 


frequency is 


L* kek; 


(IX~18) Os roras |w--O : d (k. ig k,) 


Then at sufficiently low frequencies, i.e., those approaching static or steady- 

state load conditions orn becomes a pure spring formed by the series link- 

age of the coupling spring (e/a ke and the effective hydraulic servo spring 

(2’/d) k, « The approximate relationship Coreras L el ha holds true if 

k, >>Ks , which is the case for a very high gain, closed-center hydraulic servo. 
The range of frequencies higher than the frequency peak indicated by wy 


on Figure IX-3 may also be examined. The exact value of Co in this range 


TOTAL 
cannot be obtained in general terms. However, numerical calculations for typical 
systems have substantiated the physical reasoning that the equivalent spring is 
essentially the series spring combination that would be obtained if the hydraulic 
system were viewed as a passive network. In other words, at sufficiently high 


frequencies the servo action of the system is greatly suppressed, and the effec- 


tive spring at frequencies above wW, is 
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It should be emphasized that normal static tests (e.g., loading the surface and 
measuring its deflections) will not provide a satisfactory means for obtaining 


the above spring constant. 
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Figure IX-4. Frequency Response of ~ 7; 


The typical frequency response of ~7, is shown in Figure IX-4. The moment 
component 7 x; (see IX-14 and IX-5) due to an actuator input requires additional 
consideration. The transfer function 7* cannot be classified as a complex 
spring in the same sense as Gy 7,574, , even though it has the dimensions of a 
spring. The servo input must first be knowm as a function of <,/, or /A 
before the complex stiffness concept can be extended to include / . 

The expressions (IX-15), (IX-16), and (IX-17) were obtained by considera- 
tion of a very simple type of hydraulic actuator. Many hydraulic servos fit 
into the pattern of Figure III-5, and the results of Chapters III and IV are 


available for use in investigation of more complicated cases, 
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The validity of the linear analysis in the high frequency region near WwW, 
and “y in Figures IX-3, IX-4, and IX=-5 is to be questioned, and if unusually 
high flutter frequencies are encountered a purely analytical study of the mod- 
ified flutter system is probably useless. The linear analysis fails at high 
frequency because of the increased importance of distribution of mass and damp- 
ing. 

The discussion concerning the incorporation of hydraulic system data in 
flutter calculation can be summarized as follows: 

1. Both general and approximate relationships have been presented as 
transfer functions, or portions thereof, which relate surface dynamics, 
hydraulic servo dynamics, coupling moments due to S& and h » and the 
total aerodynamic moment on the surface. 

2. A complex spring stiffness concept has also been developed, thus pro- 
viding a convenient means for analytical or graphical inclusion in the 
normal flutter solution techniques. 

3- The reductions of the complex stiffness term and the moment-input 
transfer function to effectively pure springs in certain frequency 
regimes have been presented merely to give the reader a keener physi- 
cal insight into the servo-flutter problem. The frequency variant 


expressions must be employed in the usual analysis. 


SECTION 4 — STABILITY AUGMENTER SYSTEMS 
A stability augmenter may be defined as a device which augments or creates 
desirable airplane stability characteristics by applying proper control surface 
motions and hence proper stabilizing forces and moments to the aircraft. Exam 


ples of physical stability augmenters are the pilot-artificial feel~airframe 
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system, yaw damper systems, sideslip stability augmenter system, and automatic 
pilots. 

For a stability augmenter to function in the desired way, i.e., to control 
a "controlled element," several distinct elements are required. 

The first is the "controlled element," which in this instance will always 
be the airframe. The controlled element is characterized by inputs to obtain 
control (via control deflections) and outputs generated by the input (such as 
changes in speed, angle of attack, and rate of pitch). Since the controlled 
element portion of a stability augmenter does not figure prominently at this 
point of the discussion, it will not receive detailed treatment here. 

The second element required in a stability augmenter is a device to supply 
& surface rotation. This is usually a positional servomechanism, either driv— 
ing the surface directly or, more commonly at present, driving the input member 
of a surface actuating servomechanism. The latter element is discussed in Sec= 
tion 3. 

The third important element required in a stability augmenter is the sensor, 
This is a device capable of detecting a particular motion of the controlled 
element and of transmitting a signal proportional to this motion. The sensors 
which are important in flight control system work and which may also be impor- 
tant in the modified flutter problem are: 

1. Rate gyros 

2. Accelerometer or force pick-ups 

3. Stabilized gyros 

The characteristics of these elements are summarized in later portions of 
this section. 

The fourth important element in a stability augmenter is the equalization. 


Equalization includes ail the means required (including physical connections) 
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to modify the performance of any of the system elements and of the over-all sys- 
tem to achieve satisfactory system operation. This definition is rather general, 
but, from a dynamic standpoint at least, equalization in stability augmenters 
usually takes the form of a rate circuit and an amplification. This element will 
also be discussed in this section. 

With the four prime elements of a stability augmenter defined, it is possi- 
ble to establish a fairly general pictorial representation of a stability aug~ 


menter system. 
‘ee Motion 

















Surface 
Actuator 


Positional 
Servo 






Equalization 





Airframe Motion 
~~ 


Figure IX-5. Block Diagram of a General Stability Augmenter 


From the standpoint of modifying classical flutter analysis, the presence 
of the stability augmenter creates the possibility that the sensor may pick up 
some function of a flutter motion <,/9, or 4 in addition to the "rigid" air 
frame motion it is expected to sense. In these instances, the flutter motion 
signal can result in motion of the surface, thereby creating the coupling of the 
second type mentioned in Section 2(b). 

With this information on stability augmenters as a background, it can now 
be stated that the purpose of this section is to consider in detail the mathe~- 
matical model of that portion of the stability augmenter system from sensor 
input to positional servo output. These input-output relationships will be 
developed in terms of the transfer functions of each block, with the end result 


being expressions containing positional servo output motions and sensor inputs. 


IX-21 


Section 4 


(a) SENSORS 
RATE GYROS 
A rate gyro is a gyroscopic element restrained to move in only one degree 
of freedom. If the wheel is spinning at an angular velocity «© and has a 
moment of inertia about the spin axis, J, , an angular velocity 1, about the 
x axis results in a precessing torque about the # axis as shown in Figure IX-6. 


This torque is given by the equation: 


(IxX-19) Le = L, wy IL, 


Ny 


7) 


Figure IxX-6. Rate Gyro Relationships 


Lz Zz 


For small precession angles, where the gyro is restrained about the Zz axis 


with a spring and with coulomb and viscous dampers, 
(x-20) 4, = (La,)n, = 14+ B,0+407 7 son O 


The symbol O is not to be confused with the O used in later portions of this 
report as an airframe angular rotation. 


The gyro then has a transfer function 
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(IX-21) a = ne! a 
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and 
Lawn I, 
If some sort of pickoff is then attached to give a linear electrical output 


signal, % , for a given precession angle 0, 


k 
.. (ih 2,) 


(IX-22) 
A, (3. , 2h, 
(fr+ 2%, ) 


a, 


It will be noted that a minimm detectable input angular rate due to the 


gyro construction itself exists, and is given by 


(IX-23) n, . tlhe 


MIN. I, al, 

Resolution of potentiometers or thresholds of other types of pickoffs must 
be added to this minimum signal to obtain the total minimum detectable signal. 

Since a rate gyro gives an output proportional in the steady state to 
angular rates, the airframe motions roll rate, p , pitch rate, G » and yaw rate, 
r , can be directly measured by properly mounted rate gyros. 

It must be mentioned at this point that considerable portions of the effec- 
tive angular velocity presented to the rate gyro could conceivably be due to an 
elastic motion of the airframe itself in creating local angular velocities dif- 
ferent from those of the "rigid airframe." This possibility is usually mini- 
mized by proper mounting of a rate gyro. It is obvious that the location of 


the airplane and also the frequency of any disturbing influences are of the 
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utmost importance in considering any modification to flutter analysis. This 
same conclusion is valid for the remaining stability augmenter sensors consid- 
ered. It should also be noted that a certain minimum amount of feedback is 
required to actuate the sensor. 

ACCSLEROMETERS OR FORCE PICKUPS 

Because most sensing elements used in aircraft stability work require an 
unbalanced force or moment to change some motion or other quantity of the sens- 
ing mechanism, they are in reality force or moment pickups. However, the unbal- 
anced force on an object may be a function of a velocity or displacement, etc., 
Or may be made proportional to a displacement or velocity, etc., over some fre~- 
quency range. As an example, consider the simple second order system shown sche~- 
matically in Figure IX-7. (Since most accelerometers have equivalent network 


diagrams, this diagram will suffice for all the accelerometers considered here.) 


xo f 
4. Applied oving Viscous Coulomb 
Oo orce Mass Damping k > Spring Fell Friction 
y ame ae ar ae aa 


Figure IX-7. Network Diagram of Accelerometer 


The mechanical system shown has the equation: 


(IX-24,) fF = Mx, + Bx, + kx, + F. SEN X 


3° 
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giving a transfer function when £ =O, 


2 % J 
ary F = 24, /) 
K [se + —“s+/ 
w, a), 


At various forcing frequencies, the system acts as follows: 





= ~ Xo ~ wn 
(a) G) >> wp, e F ya a)" ~ x" 
o— - Xe Ong 
(TX-27) (b) eu Fa ™ 23,K 
e Xo ~ / 
(c) a) << Gh, Z ra = ra 


It is clear from (IX-27) that the displacement x, is a measure of displacement, 


velocity, and acceleration; that is, for 
a) >> wp, » the system is a displacement pickup. 


@) = wu) 


e ? the system is a velocity pickup. 


a) << @)p, » the system is an acceleration pickup. 


Since the primary concern at present relates to the use of an accelerometer, 
it is evident that the approximation of (Ik~27c) will be of most concern; that 
is, the natural frequency should be higher than frequencies of interest in the 
system. 

With a suitable pickup for the motion of the mass, the accelerometer trans- 


fer function becomes: 








“lo” oy, 
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where ¥, is the accelerometer output signal and a; is the input acceleration. 


The minimum detectable ecceleration is given bys 


(IX-29) as ~/% . 2 


Mimn M ¢ MIN. 


To this minimum detectable value must be added the resolution or threshold of 
the pickup. 


Here again, the sensor is capable of picking up only the local input, which 


can easily be a function of a flutter motion. 
STABILIZED GYROS 
Stabilized gyros, as normally used in aircraft controllers, establish ref- 


erence planes fixed in inertial space. In the case of a vertical gyro, this is 


achieved by utilizing the gyro-torqueing system as an extremely low natural fre~- 


quency filter for a vertical sensor (basically a two degree of freedom pendulum 
or accelerometer). Gyro-torqueing is employed as a similar filter for a mag- 
netic field (or other reference) detector in the case of a directional gyro. 
From a vertical-direction gyro system, the airplane roll, pitch, and yaw angles 
can be measured directly when airframe displacements from essentially straight 
and level flight occur. 

For all practical purposes, a properly designed gyro system has a constant 


for a transfer function for fairly low airframe speeds: 


Vsyro ae 
(1X~30) » VW, oO Ks 


As mentioned above, the vertical gyro is actually nothing more than a filter 
for a pendulum, mercury switch, or other force detector, and since phugoid 


periods in seconds are about one-fifth the airspeed in mph, the dynamics of 
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the gyro-erecting system may become important. Therefore, for this case, 


/ 
Sta) 
(rr-31) Mereo _ 7 | Wag, 


= I 
TOTAL sr + s+! 
"gy Ongy 





or 
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Some instances of possible gyro-erection system coupling with the phugoid are 
already on record for high speed aircraft exhibiting phugoid periods of the order 
of 120 seconds. For all flutter modification considerations, however, the trans- 
fer function is a constant, 

As discussed above, stabilized gyros normally measure roll, pitch, and 
yaw angles. These angles are always the local angles at the gyro mounting, and 
hence the non-rigid airplane motions could affect the output of the gyro. Mounte 
ing is, again, all important as a detail consideration. 
(b) POSITIONAL SERVO SYSTEMS 

Two types of positional servo systems are important in stability augnenting 
systems. These are the electrohydraulic and the electromechanical. The electro~ 
hydraulic servo system utilizes an electrical input signal in conjunction with 
torque motors or hydraulic amplification or both to position a hydraulic valve, 
The hydraulic valve meters flow to a cylinder or to a rotary hydraulic motor, 
An electrical follow-up unit is added to provide feedback, The device operates 
in a fashion similar to the fully~powered surface actuator discussed in Chapter 
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V and has a transfer function from input signal to output motion which is very 
similar to the low frequency approximation of the fully-powered surface actuators 
K, 


Xo ae PS 
(IX-33) X; ra Usr/ 


The second type is the conventional electric motor with follow-up. As is 
well known, this type of system is approximated by: 


A 
Xo 23 

rm ee 
Ops ps 





(c) EQUALIZATION 
The equalization normally required in any stability augmenter, even in the 
most complex case, is usually derived from an electrical or mechanical rate 


circuit and amplification. The transfer function of such a device is given bys 


1; 7 s#/ 
(IX~35) ae = Ma ‘oar 


(d) THE SENSOR~EQUALIZATION-POSITIONAL SERVO COMBINATION 
Combining the characteristics of the sensor, equalization, and positional 


servo of a stability augmenting system gives a total open loop transfer function 


of the form: 
6 eee 
(1X-36) fDL,a@g ore | 
a dp 2 Sseave + / 
| Keease | ka (Gist 0 seeve — 
or oe any a, eae 
Ksauson (Es . !) Kuro 
Sy PSeameoe 5 4] Tranves? / 


2 
SENSOR Wsenson 


a ce tl 


SENSOR EQUALIZATION POSITION SERVO 


IX-28 


Pe 


= Fe Free. Fe eee ES UE ES PO Oe Oe ES 
a ell alll aol aot ell eel ol ot 
* Ph é 6 a” a a er ee ae” ee ne ee” (ce, “oe 


‘ : in 
% 


Section 4 


If ja) is substituted for s , the transfer function becomes a frequency 
response, The entire transfer function can then be changed into a real and an 


imaginary part: 


a aa! 
(IX-37) Alame = Rarenanten * J£avemenen % YR +I" omit 


In servo work, the "amplitude ratio," (R*+I rye » is usually converted into 
decibels for plotting purposes by taking 20O1oe,, (R* +I *) e The phase, 
TAN' I/R_—_, 48 used directly. The plot of typical cases of (IX-37) would look 
like Figure IX-8. The real and imaginary parts, K and J , can be easily obtained 


from such a Bode plot, 
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Figure IX-8, Frequency Response of Sensor-Equalization-Positional 
Servo Combination 
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SECTION 5 - UNBALANCED CONTROL SYSTiuM INERTIAS AND CABLE SYSTEMS 

The control systems of missiles and airplanes (particularly piloted air- 
craft) contain certain masses or inertias which are subject to forces as a result 
of accelerations. These forces usually arise from the rigid body dynamics of 
the airframe, say, in the longitudinal case, the short period oscillations of 
the airplane. A bobweight installed as a part of an artificial feel device is 
usually the most effective element for producing surface actuator input signals 
via a flexible cable system. The stick assembly inertia is another important 
element in the control system. In any event, the combination of rigid=-body 
acceleration sensitive masses, ‘system frictions, and cable or feel system flexi- 


bilities may be lumped together into a single second order system, 


y, Me. 
8 
ee 
sais . ate + | + 26-2 | 
we sas 50 
where 


X; is the surface actuator servo input displacement. 

m is the local normal acceleration. 

M, is the unbalanced mass subject to the acceleration. 

kK, is the mechanical advantage between the unbalanced mass displacement 

and the servo input displacement. 

k is the effective centering spring for the unbalanced mass. 

fo is the effective system damping ratio. 

@).= (k /M, , where A7% is the total effective mass of the control system. 

Another case where a surface actuator input can occur via a control system 
element is when rotation of the surface bends or twists the fuselage, stretching 
one portion of the cable system more than another. The transfer function relat- 


ing the surface rotation and the resultant displacement of some portion of 
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the cable system will be discussed later, 

The particular transfer function characteristics under consideration at 
this point are those caused by cable system dynamics when excited by a local 
cable input X; « A functional relationship may be developed between the dis- 
placement xe and the displacement of the mass or inertia center, X,, of the 
complete control (cable) system. This transfer function also may be represented 


as a second order expression: 
’ 
Xe 


K 
PEs) i Re ee ae REND in eee eRe eee ee 
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at 7 "@, 


T 





Since a mechanical advantage X- usually exists between the cable system and 
the surface actuator input point, the final desired transfer function is: 





(m40) =f = pF 
X, 


SECTION 6 = FEEDBACK OF A FLUTTER MOTION TO A STABILITY AUGMENTER 
SENSOR, BOBWEIGHT, ETC. 
The previous discussion of stability augmenters, bobweights, etc., has 

treated one phase of description of possible flutter couplings due to a flute 

ter motion detected by a sensor and actuating the surface. Before anything 
definite can be said about these couplings, the possibility that sensors might 
pick up some function of the flutter motion must be explored. As previously 
pointed out, the analytical formulation of this link is dependent almost entirely 
upon the installation of the sensor in the airframe, and upon the flutter motion 
transmitting characteristics of the elastic airframe. Therefore, this section 


cannot take all possible couplings into account. The possibility that any one 
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feedback may occur must be a subject for individual study on each aircraft con= 
sidered, but in many cases, proper sensor mounting could obviate the necessity 
of considering some of these links. (Proper sensor mounting would be defined, 
for this consideration, as the process of placing a sensor where it picks up 
only those quantities desired for the particular device). 

The auantities required are expressions for sensor input in terms of flute 
ter motions, i.e., 


n,a, 9 


IX= 
ve Ee 


KR, * TL, 


This expression can then be combined with the transfer characteristics of the 
stability augmenter to form a complex number representation of the particular 


coupling: 
(1x42) (R, + LVR, + 71,) = R+ GI 


The result can then be examined by the methods outlined in Section 7 of this 
chapter to determine whether a particular coupling is important. 
(a) COUPLING WITH SURFACE ROTATION 

All stability augmenters are coupled with surface rotation through the 
"rigid" body motions of the airframe. The airframe transfer functions for such 
cases are very well known* and need only be repeated here. Normally, the fre- 
quency response of the augmenter—airframe combination is concentrated in such 
a low frequency region that the effect on all but the lowest frequency flutter 
is negligible, 

Of more direct concern at present, is the coupling of the surface rotation 


and the "elastic" motions of the airframe throngh the stability augmenter. 


% See, for example, "Aircraft Controller Design," Northrop Report SMD=-31A, 
Northrop Aircraft, Inc., Hawthorne, California, October, 1952. 
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For example, consider the case where an elevator rotation, by increasing 
the horizontal tail load, bends the fuselage of the aircraft. This phenomenon 
can easily move one portion of a cable system, produce a local component of 
effective angular velocity sensed by a rate gyro, create a load force on an 
accelerometer, or cause other similar effects. Any of these effects can ulti- 
mately result in input signals to the surface actuator system and thereby affect 
surface motion. The transfer characteristic would have a first approximation as 
@ second order system, with the natural frequency being that of the first fuse~ 
lage bending mode.* In any situation where fuselage bending is an important 
mode, the possibility of such a coupling siould not be overlooked. 

Summarizing the transfer functions of importance: 

l. "Rigid Body" Airframe (Longitudinal Axis only) 



























(1-43) 2 ke (ia) CF, (70) * Miz, (7a))+ AL-%, Fw) + I] 
A =) + ede ve pas 2 S.p,. 
e we (ja) +/ we, + we, (jo) 
(IX-44,) 7 or £8 - Ke Ue (Ja) + ae i} 
—— SP 
er (id) + te 20) +1) 
2. Elastic Airframe (Fuselage Bending) 
a OR = oR O is 
(IX-45) a ; 
hd -w 2d 
(“Set Star) 


*# The basic assumption here is that the aerodynamic forces are of a quasi- 
steady nature; i.e., they are not functions of beYy as such. 


IX-33 


Section 6 


where k,, is the static amount either of upward displacement (in the case of 
the accelerometer or bobweight) or angular displacement (in the case of the 
rate gyro or vertical gyro), etc., per unit of surface rotation. 

(b) COUPLING WITH TORSION 

A change in the twist, %, of a tail surface causes a change in the local 
angle of attack and hence in pitching moment and tail load. The tail load 
change is nearly instantaneous, and consequently fuselage bending effects anal- 
ogous to those considered above are possible. The tail load and moment change 
can also excite the sensors through the airplane rigid body motions which, of 
course, are of fairly low frequency. 

(c) COUPLING WITH BENDING 

An excellent example of possible bending coupling is given in the case of 
an elevon controlled airplane with a normal force stability augmenter. The 
first wing bending mode can be of very low frequency (particularly with pods or 
stores), and is almost invariably picked up by the normal accelerometer or bob- 
weight. The transfer characteristic is again approximated by a second order 
system, with the natural frequency in this case given by that of the first wing 
bending mode. This type of coupling could be quite Serious in cases of elevon 
flutter. 

Another example is that of asymmetrical bending of a horizontal stabilizer 
giving rise to fuselage torsion, which is picked up by either a roll rate or 
roll gyro. In some cases this coupling is probably not serious, especially 
where torsional frequencies are very high, and the sensor will usually be used 
to control the ailerons. However, the downwash is changed by aileron rotation, 
so a very complex coupling of asymmetrical bending-roll gyro-aileron-downwash= 


stabilizer could conceivably occur. 
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No summary of possible bending-sensor transfer functions will be given here, 

since the form depends very intimately upon the actual configuration. 
SECTION 7 ~ TYPICAL FLUTTER SYSTEM MODIFYING TERMS 

In this section, the previously developed component characteristics are 
combined into the total flutter-servo system modifying terms for a typical case. 

For illustrative purposes, assume that the total system is composed of: 

1. The conventional flutter system. 

2. A pitch rate stability augmenter using an electrical servo and no 

equalization other than a gain. 

3. A fully-powered hydraulic surface actuator. 

In addition, assume that the only feedbacks are: 

l. A rigid body feedback due to (, 

2. An elastic feedback due to fuselage bending caused by/S . 

3. A rigid body feedback due to&, 

4. An elastic feedback due to fuselage bending caused by x, 

Physically, & and & produce airloads on the parent surface and on the 
control surface seauestively: These loads cause the airframe to execute rigid 
body oscillations and fuselage bending, both of which are detected and trans~ 
mitted by the pitch rate sensor. 

The airframe augmenter system under discussion can be represented as shown 
in Figure IX-9. The block labeled "airframe" will be different depending on 
whether the pitch rate is caused by “or (9 and whether the mode of motion under 


consideration is a rigid body mode or a fuselage bending mode. 
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Figure IX-9. Block Diagram of a Typical Stability Augmenter 


For the case where the rigid body feedback is due to /O , the modifying 


term has the form: 


(Ix-46)  %% 


] 


Kk, Kno K stavo 

(IX-47) YG LE) s* 2h, s oe 
_ 6 Cle ion? 

"ne n6 SfRVO 8fAaV0 


kK. (les +/XZ,s+/)s 
4 


sé 2¢ s§ 2&, ) 
—e. / + ——H2 5s + / 
a ap. 6 AW,” ah 


=~ %%% 


By substituting 7a) for Ss , (IX~47) becomes a frequency response, and can 


be represented as: 


(m8) Xa) %(w) = Rio) + 72(W) 


A Bode plot of (IX-47)is shown in Figure IX-10. 
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At frequencies of concern in flutter analysis, the phugoid is not important. 


The rigid body feedback rpresentation can then be approximated by Figure IX-ll. 
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Figure IX-l1. Simplified Frequency Response of Modifying Term Y; % 


i] 


For the case of the fuselage bending feedback due tof , the airframe 


portion of (IX-47) can be replaced by the fuselage bending transfer function. 
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As in the previous case, the substitution of 7a) far 3 allows the transfer 
function to be separated into its real and imaginary parts. A Bode plot of 
(IX-49) is shown in Figure IX=-12. 
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Figure IX-12. Frequency Response of Modifying Term /, Ia, 
The Bode plots shown in Figures IX-ll and IX-l2 are complete except for 


the zerodb line which establishes the system gain. The over-all system gain 
is a direct function of the gain of each of the component transfer functions. 
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These gains are of course functions of the dynamical characteristics of the 
particulaP-systen being investigated. In the general case, the gain can only 
be approximated, 

In the typical case under consideration in this section, the Bode plots for 
the rigid body feedbacks due to & and 2) differ only by the location of the 
zero db line; that is, the motion caused by a given deflection of either the 
parent surface or the control surface will be the same except for the amplitude, 
This is also true for the case of the fuselage bending feedback. Thus 
¥; (w) Yo, (w) » %@) Me, (Ww) , Ye (ad) % (a) , and % (@) a ()) can be evaluated 
for a particular frequency from Figures IX-l] and IX-12 with the appropriate 
zero db lines. Once the real and imaginary parts of these terms have been deter- 
mined, they can be included in (IX-9). The magnitude of these modifying terms 
compared to the classical flutter terms determines the influence of the stabil~ 
ity augmenter system on the flutter system. 

The generic forms of the Bode plots for stability augmenter-flutter systems 
using other types of sensors can be inferred by inspection of Figures IX-10 and 
IX-11, For instance, a system using a stabilized gyro rather than a pitch rate 
gyro would sense pitch angle instead of pitch rate. Inthe case of either the 
rigid body or the fuselage bending feedback, the plots in Figures IX-10 and IX-l1 
would be altered by the addition of a slope of ~-6 db per octave to each of. the 
asymptotes, and the entire phase curve would be shifted -90°. For example, the 
Bode plot of Figure IX-12 altered to pertain to a system using a stabilized 


gyro would appear as shown in Figure IX~13,. 
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Figure IX-13. Frequency Response of Modifying Term / b, 
for a System Using a Stabilized Gyro 


Expressions that account for other types of feedbacks, such as hydraulic 
valve motion due to cable system movement, can be developed in the manner used 
in deriving (IX-9). 

The technique of modifying the classical flutter equations to account for 
the presence of complex servomechanism systems has been presented. Chapter X 
will develop the mechanics of solving these modified flutter equations, 1i.ée., 
the methods of handling Re))+ 7) ()) modifying terms. The cases where feedback 
path transfer functions become functions of bwi/v » thereby requiring additional 


degrees of freedom, will also be discussed, 
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CHAPTER X 
METHODS OF ANALYSIS OF SERVO-FLUTTER INTERACTION* 


SECTION 1 ~- INTRODUCTION 


This chapter considers the theoretical determination of the flutter charac- 
teristics of a coupled servo-airframe system. The method of attack is the 
successive consideration of three increasingly complex systems. This approach 
is chosen as being of particular value in gaining an insight into the manner 
in which the servomechaniam affects the characteristics of a flutter system. It 
is believed that the method to be used here yields a practical approach to this 
problem, without departing too radically from conventional methods of flutter 
analysis. By progressing from a simple system to ones in which more complex 
problems of interaction exist, a physical feel for the matual influences of the 
servo and the flutter system may be obtained. 

The three type of systems to be considered are: 

1. Single degree of freedom, no structural feedback to servo input. 

2. More than one degree of freedom, no structural feedback to servo 
input. 

3. More than one degree of freedom, with structural feedback to 
servo input. 

In (1), the only degree of freedom to be considered is the rotation of the 
control surface; no deformations of the control surface, the parent surface, or 
other structures are taken into account. 

In (2), the system is to be regarded as having, besides a degree of freedom 
associated with control surface rotation, additional degrees of freedom which 
express any combination of the following: Parent surface bending and torsion, 


fuselage bending and torsion, and control surface torsion. None of the 


*Portions of the material contained in this chapter appeared previously 
in Air Force Technical Report #6287. Special permission for its use here has 
been granted by Wright Air Development Center. 
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additional degrees of freedom are to cause any appreciable input to the servo. 
In (3), the system is of the type mentioned under (2), but modified by 
structural feedback which provides an additional input to the servo, with no 


specification as to how or where the feedback originates. 


SECTION 2 - BASIC CONCEPTS 

The flutter characteristics of a coupled servo-airframe system may differ 
considerably from those of a conventional flutter system. A physical under- 
standing of these differences is facilitated by the introduction of three 
relatively new concepts: 

1. Frequency Dependent Restraint 

2. "Complex Stiffness" 

3. “Feedback Uncoupling" 
In addition to their usefulness as aids to intuitive thinking, these concepts 
have a basic utility in analysis. It is worthwhile, then, to begin the develop- 
ment of methods of analysis with a discussion of these very important concepts. 
(a) FREQUENCY DEPENDENT RESTRAINT 

In the conventional flutter system, the resistance offered by the control 
system to control surface rotation is considered to be that provided by a spring 
of fixed stiffness. The definition of the control system is somewhat arbitrary. 
It has, in general, several inputs, one of which is the control surface rotation. 
The determination, either experimentally or by analysis, of a spring constant to 
be used in the flutter analysis involves some sort of specification of all control 
system inputs except the control surface rotation. Once the spring constant is 
established, it is considered to be applicable for any frequency at which flutter 
may occur. ‘The control surface thus has only one natural frequency of rotation, 
and this natural frequency is invariant with the frequency of a forced oscillation 


of the control surface. 
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For the more complex servo-actuator control system, this simplification is 
generally not valid; in this case the effective stiffness becomes a complex 
function of frequency. Why this is true can be seen as follows: Consider that 
the control surface is removed and that the actuator is grasped at its connec- 
tion to the control surface horn and its output shaft oscillated with some 
fixed amplitude and at a fixed frequency. The servo will attempt to hold its 
output position fixed against the action of this load disturbance. However, 
due to the nature of the servo system components (electric motors, hydraulic 
fluid, distributed and concentrated masses and springs), and to the fact that 
the servo is an error-operated, closed-loop device, the force output of the 
servo will differ from the disturbance both in amplitude and phase. 

In general, the amplitude and phase of the servo response under the above 
conditions will be functions of the frequency with which the load disturbance 
oscillates. This indicates that the restraint applied by the servo system to 
the control surface cannot be adequately represented by the fixed spring of 
the conventional flutter analysis. 

As an example of the possibly serious errors which might be introduced 
into a flutter analysis when the frequency dependence of restraint is not 
taken into account, consider the following: 

Suppose that a single value of rotational frequency ie measured for a 
given servo actuator--control surface combination, and that this is psed in 
a conventional flutter analysis. Then, unless the resultant flutter fre- 
quency is very close to the measured natural frequency, the restraint pro- 
vided by the servo system at the flutter frequency will not be of the correct 
amplitude nor in the correct phase to afford adequate control of the surface 
-at this frequency. It may actually occur that although calculations based on 


this single measured frequency indicate stable flutter characteristics, the 
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system becomes unstable in practice at some other frequency. This will be the 

case if the control system restraint of the surface at this other frequency is 
sufficiently reduced, or if its phase is markedly changed, by the frequency depend- 
ence. 

It will be noted that both amplitude and phase changes as a result of varying 
frequency were mentioned in the preceding discussion. If only the amplitude were 
affected, the servo-actuator system could be regarded as merely a variable spring 
constant. However, the mention of variation in phase with load disturbance fre- 
quency implies that damping (positive or negative) occurs; that is, that the 
system contains an energy dissipating or an energy producing device. The magnitude 
of the energy absorbed, or supplied, will also in general vary with frequency. 

(b) “COMPLEX STIFFNESS" 

The second ccncept to be introduced is that of the "complex stiffness" of the 
control system. This is a quantity which can be used to describe the properties 
of the control system from its input (or inputs) to the point at which the actuator 
is connected to the control surface. It may be defined as the steady-state reac- 
tion on the control surface, exerted by the control system, due to control surface 
motion at fixed amplitude and at a fixed frequency. It is assumed in this defini- 
tion that there is no input to the control system other than that due to the sur- 
face load disturbance. Since the servo itself is an error-actuated, closed-loop 
device, 1t will have an output to counteract this load disturbance, but this will 
in general differ in amplitude and phase from the disturbance. The complex stiff- 
ness thus has a real and an imaginary part, each of which is a function of the 
frequency with which the surface is oscillated. 

The complex stiffness of the control system can be experimentally determined 


on the ground with the control system actuated only by the control surface motion. 
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Control surface inertia and aerodynamic forces play no part in determining the 
complex stiffness. By using the concept of complex stiffness, analyses of 
conventional and of servo control systems can be performed in the same general 
way. The stiffness characteristics of the control system can be evaluated in 
& routine manner. With conventional control systems, the natural frequency 
can be measured. With the servo type system, a single measurement of natural 
frequency is not sufficient, but the characteristic for this type can be 
determined by a frequency response test. Calculated or measured stiffness 
characteristics can then be supplied to the flutter specialist in much the 
same way that basic weight or stiffness data are now supplied. 

With the approach discussed above, the servo engineer need not become a 
flutter expert, and the flutter engineer need not become a specialist in cal- 

ulating the effective stiffness of a control system of either conventional 
of servo type. Each can concentrate on his own phase of the theoretical 
uivestigation and use his experience in his own field efficiently. Their 
activities correspond to regions of the physical system which join at the 
connection of actuator and control surface. 

(c) "FEEDBACK UNCOUPLING" 

The third concept concerns a method of attack for the case where feedback 
arising from other degrees of freedom may, through the servo actuator, result 
in control surface response. The input to the actuator may result from elastic 
deformation of cables or bellcranks, from inertia forces due to concentrated 
masses, or from the response of sensing devices and servo transducers. Actuator 
inputs may thus attend an elastic or rigid body degree of freedom in flutter. 
Hence, if a certain actuator input results from a given degree of frecdom, 
whether it is wing first bending or airplane yaw, this particular input is 


associated with a fixed shape of the deflection distribution over the airframe, 
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which in general does not include control surface rotation. The result is that 
the control surface, because of its fixed position (i.e.,2=O ), is subjected 
to actuator imposed forces (or moments) as a consequence of feedback. The forces 
which thus arise from a given degree of freedom (such as wing torsion) are then 
capable of doing work in the control surface degree of freedom. This is analo= 
gous to structural coupling between two generalized coordinates, 

To achieve what is in effect feedback uncoupling, it is necessary to con- 
sider the previously mentioned wing torsion degree of freedom as including con- 
trol surface rotation of such relative magnitude and phase as would be deter- 
mined by the condition that the actuator force applied to the control surface 
is zero. 

Figure X=-1 shows simple systems which exhibit these characteristics. Motion 
is restricted to control surface rotation and parent surface bending (or trans- 


lation). 






——Undeflected 


—Undeflected 


——_ 


_Deflected 


Figure X-1.. Simple Systems 


Figure X-1(a) represents the general case of a more conventional system, or, 


more explicitly, the case where the control surface is not subjected to internal 
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forces (or moments) due to bending of the parent surface. In this degree of 
freedom, there is no rotation. It should be noted that this system is also 
representative of the case of a servo actuator with no feedback. 

Figure X-1(b) illustrates the effect of feedback. As a normal consequence 
of the conventional type of bending mode / (no rotation ), the control sur- 
face is subjected to an internal moment due to the stretching of the attached 
spring. This is similar to the case where bending motion results in actuator 
inputs, and hence actuator imposed moments on the surface. 

Figure X-1(c) illustrates the case shown in Figure X-1(b) except that the 
bending degree of freedom 4 is altered to include a proportionate amount of 
rotation 4: such that the attached spring does not stretch. Hence, no spring 
(or feedback) coupling occurs in the adjusted bending degree of freedom. 

In the simple case shown, there is no phase angle between the bending 
motion / and the rotation 4; for the directions indicated as. positive in 
Figure X-1. Also, the ratio of 6 to / is a constant and is independent of 
the frequency of oscillation. However, in the case of a practical servo system, 
due to the nature of the mechanical, electrical, hydraulic, and structural ele- 
ments between the source of the feedback and the control surface proper, both 
the amplitude ratio and the phase angle between the degree of freedom which 
causes feedback and the "induced" control surface rotation (; are in general 
functions of the frequency of oscillation. 

In general, then, associated with a given degree of freedom 7; is a con- 
trol surface rotation #6; , which is determined by the condition that the 
internal force applied to the control surface is zero. Hence, Pet = FR; Gi 
where Ay is 4 complex function of frequency. 

It should be noted that, in accordance with current techniques of flutter 


analysis, the Se derree of frcedom (e.g., wing second symuetrical torsi on) is 
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fixed insofar as the shape of its deflection distribution is concerned and is 
assumed to contribute to the flutter mode, regardless of the flutter frequency, 
without changing its shape. Therefore, use of this concept implies the need for 
determining the relationship between the rotation (; and the deflection ¢- 

for a fixed mode shape and any required frequency. This is not a difficult prob- 
lem if the major contribution to the servo input is from a single sensing device, 
such as an accelerometer, because the result of acceleration at the sensor loca~ 
tion can be either calculated or measured. However, when the feedback is the 
result of combined inertia, elastic, and electrical effects and is dependent on 
both the deflection distribution and the frequency, experimental determination 
of the physical phenomenon represented by the symbol 7; is difficult. 

As an example, consider the case where A; due to an elastic degree of free- 
dom, say wing bending, is to be determined in a ground test. The control surface 
is removed to meet the necessary condition of zero force at the horn. The servo 
actuator is then connected to a light weight link which does not restrict motion 
of the actuator extension. With proper instrumentation, the actuator extension 
Be and the bending motion at a reference point on the wing can be measured. How- 
ever, except at the natural frequency of the wing bending mode, it is not possible 
to excite the structure at frequencies much different from this natural frequency 
and still retain the same wing bending mode shape. Yet, this is required if 
compatibility is to be achieved between the resultant flutter frequency and the 
feedback contribution at the flutter frequency. 

The difficulty mentioned above, however, does not necessarily restrict the 
use of this method of attack insofar as the application of experimental results 
is concerned, Experimental measurements can be made at and near resonant fre- 
quencies, and a combination of calculations and engincering judgment can be used 
to arrive at the character of A; at frequencies above and below the measured 


frequency. 
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The conditions of fixing the mode shape (if the feedback system is dependent 
on it) and varying the frequency of oscillation can probably be realized in a 
completely analytical approach to the determination of RK, as a function of fre- 
quency. 

In the event that there are pressure sensing devices which influence the 
servo action, and hence the amount of feedback resulting from a given degree of 
freedom, it may be necessary to determine the function 4; as a function of 
dynamic pressure, in addition to frequency. In this case, there is the additional 
complication that compatibility must be achieved between the dynamic pressure at 
flutter and the feedback associated with the flutter conditions. To determine 
R: experimentally in a ground test, it may be necessary to measure 4; at each 


frequency for each of several typical pressures applied to the sensing device. 


SECTION 3 - ANALYTICAL PROCEDURES 

The basic concepts discussed in Section 2 are used in the following analyti- 
cal treatments of the three systems described in Section l. 

(a) ONE DEGREE OF FREEDOM FLUTTER (NO FEEDBACK) 

The one degree of freedom flutter case involving only rigid control surface 
rotation will be considered. Structural twist of the surface is most conveniently 
handled by considering it as a separate degree of freedom. This is discussed in 
Section 3(b). 

In a manner similar to that employed in current flutter practice, a fic- 
titious damping will be assumed to exist at the control surface horn. The 
amount of this damping required to maintain steady state oscillations will be 
determined. A negative value of this damping corresponds to an excitation, i.e., 
to the presence of a fictitious "energy providing" device. When negutive values 
result, it means that an exciter is necessary to maintain steady state motion. 
Hence, lacking this exciter, motion would decay. If positive values result, a 
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damper would be required. Hence, without this required damper, divergence will 
occur. A zero value of this fictitious damping indicates that the system is 
capable of constant amplitude flutter without the benefit of either of these 
fictitious devices. 


Consider the case shown in Figure X-2: 


”\ Aerodynamic Hinge Moment 
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~_ Ts w Torque Due to Control 
System Restraint 


~—_lg Torque Due to 
Fictitious Damper 


Figure X-2. One Degree of Freedom Case with Fictitious Damper 


The torque due to the restraint of the servo control system as a result of 
constant amplitude oscillatory displacements /S can be expressed (in the linear 


range) as 
(X-1) i, = Kf 


It is important to note that iL, is the reaction supplied by the control 
system as a result of motion imposed at the control surface. Any servo input 
that may result would occur due to reflection through the input system as a 
result of control surface driving motion. The influence of any feedback that 
may result due to control surface rotation only, for example, feedback through 
a sensing element attached to control surface, would be integrated within the 


complex stiffness discussed in Section 2(b). 
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The complex stiffness K * is in general made up of real and imaginary 


parts. Hence, 


I, 
(X-2) oe = K,+ Jk, 
where 
K, = Kk, (@) 
(X=3) 
K, = K, (a) 


The torque due to the fictitious damper can be considered proportional to 
the amplitude @ and, in the sense shown, leading it by 90°. Hence, 


(X-4) Lt =7Ae 


where c is some arbitrary constant. 

If it is desired. to consider the damping action as viscous, and the torque, 
therefore, proportional to velocity, it should be noted that for steady state 
harmonic motion, 7c is replaced by jwAc. Thus, C= wc’ . 

Inasmuch as the object is to find the value of c required for steady 
state oscillatory motion, the resultant flutter frequency will allow a value of 
c” to be readily obtained. In addition, here, the realistic case is considered 
to hold for c= O; hence, also for c=O. Therefore, the concept of a damping 
force proportional to amplitude rather than velocity (whether it is structural 
damping or viscous damping), as employed in classical flutter analysis, need not 
be considered restrictive. It can be thought of as a convenient artifice 
employed to obtain equations of a form which can be solved more readily without 


necessarily sacrificing physical significance. 


*For the hydraulic actuator treated in Section 3 of Chapter 9, this complex 


stiffness is denoted by Os ° Ce * Osrorac 
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The aerodynamic hinge moment due to/S8 motion only is obtained from page 34 
of Reference*, Integrated over the span, for constant f distribution, the total 


aerodynamic hinge moment is 
(X-5) L= ftw Q 
where 
Q = no [b* [To - (Po + Te le-e') +R (c-e) Jax 


The equation for the one degree flutter case is obtained by the use of 


Lagrange's equations, 
at 
where ° =/% e ‘ 
| se 
Kinetic energy = > fa ve [, ax 


Potential energy=0O (Due to the absence of structural deformation in this 
degree of freedom) 


Work done by external forces due to virtual displacement,5, , is 
(Ty - Ty - Ta) So 


The resultant equation is 
2 
(X-6) -w SfI,dx- T+ Ly = O 
or, substituting (X-1), (X=-2), (X-4), and (X=-5)into (X-6), 
(X-7) -w/8 [I,dx + 3(Kp+7K;) t 7 fie - 30 Q= O 


A different measure of the fictitious damping will be introduced for con-= 


venience. The relationship between this new damping g, and c is 
Cc 
(X-8) hL= Ke 


* Smilg, B., and Wasserman, L. S., Application of Three-Dimensional Flutter 
theory to Aircraft Structures, Air Corps Technical Report No. 4798, War Department, 


Air Corps, iiateriel Division, Dayton, Ohio, July 9, 19h2. 
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It should be noted that for K,>O, the quantity g has the same physical mean- 
ing as has c., 
Substituting (X-8) into (X-7) and dividing by q)/3 , the equation for the 


One degree case becomes 


(X=9) to - “5 (1 +7g%)+Q 20 


where 
Je, * [Jo dx 


gq = Erg, 
Re 

Equation (X-9) is, of course, essentially no different from the more famil- 
iar control surface one degree of freedom equation (see page 63 of Air Corps 
Technical Report No. 4798), inasmuch as a concentrated real spring of stiffness 
K, is equivalent to writing we 1, . Also, for g=O, K,/Kp is equivalent 
to the available structural damping. 

A solution of (X-9) is accomplished by solving for values of k, Iw and 
G that are required at each value of v/bw . From a knowledge of the char- 
acter of A, as a function of @ , the compatible value of Ww results. From 
this value of @), gd can be obtained by subtracting the available amount of 
damping (i.e., A;/A, ) from the required value Je + The corresponding air~ 
speed is obtained at each u/hw by using the compatible value of @. By taking 
enough values of y/h@, a plot of Jo VB~ V can be made in the usual manner. 

The flutter speed will occur when g*O . The detailed procedure is 
described below: 

1. Obtain plots of k, Jay vs. @ and K,/Kp, vs. @ by calculation, meas- 


urement, or any other means desired. 


X~13 


Section 3 


26 
36 


De 


Te 


8. 


Je 
10. 


Select a value of y/bu) « This fixes the complex number Q=Q, +7Qr. 
Compute the required value of K, /w* from (X-9). This value is equal 
to Jy + Q, « Compute required value of 9, from (X-9). 

From plot of K, /w* vs. @ , pick off value of @ at which the required 
value of X,/ w is provided. (If more than one solution exists, obtain 
each value of @. In general, there will probably be an odd number of 
solutions.) 

Using this value of w@ , read the value of K,/Kep that is available 
from the plot of K; /Ke VS. @). 

Using this value of X, / K, and the required value of J » compute gy . 
Obtain corresponding value of V from V/ba , 4, and@. 

Plot the values of V and g on a plot of g vs.V. If step (4) 
resulted in 3 solutions, there will be three points on this plot. 
Repeat steps (2) through (8) for several values of u/b«). 

Pick off flutter speed from plot of g vs. V at point corresponding 
to g.= O in the usual manner, At speeds corresponding to negative 
values of J, » the system is stable. At speeds corresponding to posi- 


tive values of J, » divergent flutter will result. 


(b) TWO DEGREE OF FREEDOM FLUTTER (NO FEEDBACK) 


The one degree of freedom flutter case discussed in Section 3(a) will be 


extended to include any other elastic degree of freedom. (Note: The addition 


of a rigid body degree of freedom, e.g., airplane pitching, does not necessitate 


a change in the character of the solution for the one degree case.) 


The method of solution of the flutter case involving mre than two elastic 


degrees of freedom (no feedback) is discussed in Section 3(c). 
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The two degree case treated here will consist of: 

a. Control surface rotation (as a rigid body) due to control system 

elasticity (i.e., servo elasticity). | 

b. Any other elastic degree of freedom, such as parent surface (or 

fuselage) bending or torsion or control surface structural torsion. 
It is assumed that deformations in this degree of freedom do not give 
rise to servo inputs. The feedback case is discussed in Section 3(d). 

It should be noted that attempting to account for control surface structural 
torsion by combining it with control surface rotation as a coupled degree of 
freedom is not recommended for the following reason: 

The implicit assumption made in using the coupled degree of freedom is 

that the ratio of the magnitude of the structural torsion degree of free- 

dom to the magnitude of the rotational degree of freedom is constant for 

all (flutter) frequencies, This assumption does not hold for the general 
case and could lead to more erroneous results for the servo case than for 
the conventional case, depending on the nature of the frequency dependent 
stiffness characteristic of the control system. 

However, sufficient experience with practical systems may lead to the gen=- 
eral conclusion that a coupled degree of freedom can be employed ‘in an analysis 
involving a servo with a tolerable error in results. Until such time, however, 
it appears that a theoretical investigation of a flutter mode which may involve 
a variation in/9 along the span of the aerodynamic control surface, arising 
from structural torsion of the control surface,* is more likely to yield real- 
istic results if the structural effects are included as a separate degree of 
freedom (i.e., rigid control system). 

"4 This is to be distinguished from the well known phenomenon of a variation 


in @ resulting from structural torsion of the parent surface alone. This effect 
can be included in the usual manner. 
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A parent surface bending degree of freedom will be added to the one degree 
case of Section 3(a) with the understanding that the analytical procedure that 
follows is general for any additional flexible degree of freedom. Then, from 
page 61 of Air Corps Technical Report No. 47968 the following two equations 
result: | 

Ah, * C 3,07 = O 
Gh, Ip, e’™= 0 


where, from page 63 of referenced report, for n=O , and £,2/ , 


2 ; . i 
(X~10) A - [ = (Drsg)| [asta + molt bh dx 
- “4 te ‘ 
(xn)  € = [Soha ¢ nf [L, - Le (c-e)]b’Kax 
—~ ~_'¢ ov) ’ 
(X-12) G = { Shdx+ nef [h- Rlc-el|b hd 


and, for the case considered here, I is given by (X~9): 
(X~13) I=, - Fe (+3 )* Q 

= T a)? I Ge 
For the determinant to vanish, it is required that 
(X-14) I = S& 


Now, for a given airframe, and at a fixed value of y/bw » the following 


terms can be represented by fixed vectors. Then, letting 


(X-15) y =4, + Q 
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a es / 
(116) 8 Fagan 
< 
ee U4 l, as 
(X-17) V [ t+ wef 


and substituting (X-10), (X-11), (X-12), (X-13), (X-15), (X-16), and (X-17) in 
(X=14), 


K, 2. i Vv, 
(X-18) ee e = V, 


- aN + 79,) + Vs 

The solution to (X-18) is obtained by finding the frequency @, at which 
compatibility is achieved between the required value of K, / «)* » represented by 
the real part of the vector V and the available value of K,/w)* e Then, in 
the same manner as for the one degree of freedom flutter case discussed in Sec- 
tion 3(a), the required value of J, is compared to the available value of 
K,/ K, at this compatible frequency. The difference between 9, and K; /Ke 
represents a measure of the stability or instability at this value of u/ba) 3 
and, hence, at the associated airspeed. 

Now, for a given value of w, , and a selected value of g, , (both being 
fixed characteristics of the known airframe) the right hand side of (X-18) (i.e., 
the vector ¥, ) would be a function of). Also, the available value of A, , 
in general, would be a function of @). Therefore, it is clear that, unless 
kK, can be conveniently expressed in terms of @, the required solution will 


involve either a trial and error procedure or a graphical procedure. 
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Furthermore, if G), is not definitely established (e.g., during the early 
design stages or for a general study), or if the nature of K, changes as the 
design progresses, it may be necessary to repeat a major portion of the work 
for each change of G@, or K,. With this in mind, it would be desirable to 
obtain results which are essentially independent of the magnitude of 4, or 
of the nature of K, vs. @, and which, for a minimm of additional labor, could 
be translated into results corresponding toa specified value of @), or a speci- 
fied character of XK, . 

The following procedure, which consists of three steps, is presented as 
one manner in which the above might be achieved. 

Step l 

For a given value of u/bé) , and at an assumed or specified value of 

structural damping, 9, , calculate the vector V, for several values of 

a), /w) from (X-18). This could be accomplished in more than one way. A. 

graphical method could employ the Bleakney circle which is discussed in 

Air Corps Technical Report No. 4798. The real part of V, represents the 

required value of XK, a » and the imaginary part of V, divided by the 

real part of represents the required value of J, necessary for con~ 
stant amplitude flutter to occur at each value of a), /.) . Repeat for as 
many values of x / bw) as deemed necessary. 

These results are then plotted as shown in Figure X-3 where curves repre=- 
senting results for two values of u / ba) are shown. 

These plots of required values of Kp |) and q, » once obtained, are 
fixed, 
The second step is independent of the first step and is limited to the 


determination of the complex stiffness as a function of @). 
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Figure X-3, 


Step 2 


Determine X, /w)* and A, /K p 8 functions of @). These characteristics 


of the control system might be determined by analysis, #* by experimental 


means, estimated, assumed, etc. In any event, these functions, which are 


% Numerical subscripts denote particular values of y/ bw). 


*## The expression for the complex stiffness for a typical hydraulic actuator 
is developed in Section 3 of Chapter IX. 
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then plotted as shown in Figure X-4, are intended to represent the available 
values of Ky | and K,/K, « 


Kp 
“we 
a) —— » 
Ky 
Kp 
QQ ——-» 
Step 3 


In the third step, the plots resulting from Step 1 and those obtained in 

Step 2 are combined in one specias chart. This chart allows for the nec- 
essary compatibility to be realized between the required value of K, /a', 
the available value of X, /u)* , the associated value of @) , and any speci- 
fied value of @, , The corresponding values of q, Ke /K eg», and V are 


read directly from this chart. 
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This chart is made up of four basic types of curves as illustrated in 
Figure X-5. These are as follows: 

Type (1) curves are shown as solid lines. They are dependent only on the 

ordinate scale, and consist of: 

1. A single hyperbolic curve which is a plot of @ vs. //) . This is 
shown in the left center part of the chart. 

2. Straight lines, each of which represents a fixed value of @, . These 
are located in the right center part of the plot. (Only one is shown 
in Figure X-5 but any number may be drawn). The tangent of the angle 
between this line and the vertical is equal to 4, multiplied by a 
constant which is determined by the relative scales of the @), / ws) and 
//) ordinates. 

Type (2) curves are shown as solid lines. These are results of solving 

(x-18) (4.e., plots obtained from Step 1). One of each of the following 3 

curves for each value of v/bW . (Only one set is shown in Figure X-5, 

but curves from all values of v/ba) may be drawn on the same chart). 

1. Curve of (K, / uw") VS. (a), /w) .» These are located in the upper right 
hand part of the chart. 

2. Curve of I, vs. (w/a). These are located in the lower right hand 
part of the chart. 

3. Straight line of V vs. @. These are located in the left center part 
of the chart. 

Type (3) curves are shown as dashed lines. These are characteristics of 

the control system (i.e., plots obtained from Step 2): 

1. Curve of ( K, /@)) vs. @ . This is shown in the upper left hand part of 


the chart. 
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2. Curve of K,/K, ve. 4). This is shown in the lower left hand part of 
the chart. 
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Figure X-5. Chart Showing Procedure for Constructing "Compatibility Curves" 
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Type (4) curves are shown as dashed lines. These are "compatibility curves", 
one for each value of /ba) , and are located in the right center part of 
chart. (Only one is shown): 

1. Curve of (a), jw) vs. ( '/w). A point, the locus of which describes this 
curve, is located by straight line projections which form an imaginary 
rectangle between the (K, /ay ) required curve, the ( K, Jw) available 
curve, the hyperbolic curve, and the point in question. An example of 
the procedure followed in locating any one point on this curve is shown 
in Figure X-5. 

The procedure which, for a given value of a, , yields values of the airspeed 

V, the required damping 7, , the available damping K, /Kp » and the associated 
frequency a) , for each value of v/ba) for which compatibility is achieved, is 
as follows: (Note: There may be more than one solution as mentioned in Section 
3(a).) 

For curves corresponding to /ba) , (see Figure X-6). 

1. First proceed along line corresponding to desired value of a) to 
intersection with "compatibility curve'’ shown as point "a" on Figure 
X-6. 

2. Project vertically down to intersection of J, curve (point "b'"). Read 
value of Jo ° 

3. From point "a", project horizontally left to intersection of hyperbolic 
curve (point "c!"). From point "c't up, read @. 

i, From point "c" down to K, /K, curve (point "d'?), Read available 
K./ Kp value. 

5. From point "c" down to airspeed (straight) line (point "e""). Read 


value of V. 
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Chart Showing Procedure for Obtaining Results 
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Repeat steps (1) thru (5) for curves corresponding to (u/ ba)), ; (v/bx)), ’ 

etc. 

Now, for each value of uv/bw , corresponding values of Vg, , KX, Kp» 
and @) have been obtained. As previously mentioned, more than one set of these 
values are possible for any single value of /s. 

The fictitious damping, g , is found from (X-9) to be 

R 

At this point a plot of V vs. g (for fixed value of @, ) can be made 


as shown below. 
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Figure X-7, 


It should be noted that: 
1. An arbitrary change in the value of 6, requires a relatively small 


amount of labor to obtain a corresponding V vs. g, curve. 
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Repeat steps (1) thru (5) for curves corresponding to (u/ ba)), F (v/bw), ? 

etc. 

Now, for each value of u/hbe) , corresponding values of lV, Je 9 Ky [Kes 
and @ have been obtained. As previously mentioned, more than one set of these 
values are possible for any single value of U/sw. 

The fictitious damping, g , is found from (X-9) to be 

R 

At this point a plot of V vs. g (for fixed value of 4) ) can be made 


as shown below. 


| 
| 
bo | 
Speed Range Were Negative | 
Additional Damping (i.e., excitation) ,|—| Speed Range 
is Required for Flutter When Additional-~~ 
Damping is Required 
i 
- 
- Vv a 
Speed at Which Flutter 
Would Occur for Existing 
System 
Figure X-7, 


It should be noted that: 
1. An arbitrary change in the value of @, requires a relatively small 


amount of labor to obtain a corresponding V vs. J, curve. 
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2. A change in the nature of the (available) curve sequites the con- 
struction of new "compatibility" curves, but the curves resulting from 
the analyses remain unchanged. In this way, several prints of the 
chart in which Type (3) and Type (4) curves are excluded could be made. 
Then the effect of changes in the nature of the complex stiffness (Lec, 
A, and A, ) could be obtained graphically from these charts with a 
minimum of effort. 

(c) THREE (OR MORE) DEGREES OF FREEDOM (NO FEEDBACK) 

The three degree of freedom flutter case is a normal extension of the two 
degree case discussed in Section 3(b). The additional labor required to include 
the third degree of freedom here is comparable to that required in adding a 
third degree of freedom to a conventional Bleakney circle analysis, | 

A parent surface torsional degree of freedom will be added to the two 
degree case of Section 3(b) with the understanding that any flexible degree of 
freedom may be added and the same procedure would still apply. 

Page 61 of Air Corps Technical Report No. 4798 gives the following three 


equations: 
Ah, ¢ Be, et + CA,0% =O 
Dh, + Ex, d* + F 3,2" =O 
Gh, + Hac +T Qed = O 


=m * 


Page 63 of the same reference gives the equations for A ’ B -Cy D 


IM; 


3 
F ,G, and H . These will not be repeated here. The equation for J is given 


by (X~-13). 


— 
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For the determinant to vanish, it is required that, 


(x-19) I= 





This is a more complex form of (X~-l4). However, by speci fying values for 
Jy 9» So 9 and (4), /a,.)s (X-19) can be written as follows: 


Ke | 
(X=20) wal $ i9,] = V. 


where, for a fixed value of(V/bw), the vector V- becomes a function of (a), /a)) 
or (a) /@)), whichever is desired. Then from this point on, the procedure is 
identical to that outlined in Section 3(b) for the two degree case. The essen~ 
tial difference is that the evaluation of V- is more laborious. One graphical 
method of computing 4 is by the extension of the Blealmey circle method to 
include the third degree of freedom whereby the term (see (X-19)) 


(o- Ae) 


ya 


becomes a vector which is added to the vector obtained from the Blealmey circle 
plot for a given value of (4), /uw)) or(«)./a)). This is a well know method of 
solution for a case of this type .and will not be discussed in detail here, 

The extension to more than three degrees of freedom is, of course, even 
more laborious, but the same procedure will still apply. 

The foregoing has been concerned with the case where specified magnitudes 
of the structural frequencies were not required until the major portion of the 
work was completed. In this section, however, the ratios of structural fre- 


quencies relative to each other are necessary. In the event that.all of the 


X~27 


Section 3 


structural frequencies are known prior to any analyses, the chart would not be 
required, and the solution would become somewhat more direct. 
(d) THE FLUTTER SYSTEM WITH FEEDBACK 

The equations for the system in which motion in a degree of freedom other 
than the control surface degree of freedom results in feedback coupling through 
the control system are developed by employing the concept of "feedback uncoupling" 
discussed in Section 2(c). 

The feedback terms, i.e., R; » employed in this treatment are taken as those 
which do not involve aerodynamic couplings; hence there is no functional depend- 
ence on the reduced frequency u/beo . AS an example, consider the two cases 
described in Section 1 of Chapter IX. In the second of these cases, torsion, <, 
causes a change in the aerodynamic moment on the aircraft and hence accelerations 
which are detected by a sensor and in turn result in actuator inputs. If the < 
degree of freedom involves no motion of the fuselage (e.g., cantilever torsion of 
wing or stablilizer) and if the accelerometer or bobweight is located within the 
fuselage, then no actuator input is directly associated with the « degree of 
freedom. However, if the « mode is a free-free symmetric mode involving fuselage 
pitching or a coupled mode involving any fuselage motion, there is a direct 
relationship between accelerations at the sensor and the torsion <. Then, any 
actuator input arising from sensing elements within the fuselage is associated 
with zero airspeed modes involving fuselage motion (i.e., airplune rigid body 
pitching, translation, or roll, fuselage flexible modes). also, if a sensing 
element is in the wing, any actuator input resulting from sensor excitation is. 


associated with only those degrees of freedom involving wing motion. 
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Associated with some degree of freedom (e.g., fuselage side bending) is a 
control surface rotation which is determined by the condition that zero 
torque is transmitted to the control surface by the actuator. 

The magnitude of (5, due to the feedback which arises from this fuselage 
side bending mode is not influenced by control surface inertia effects. (Notes 
In actual practice, the feedback terms could as readily be obtained, either 
analytically or experimentally, with the control surface attached. However it 
would then be necessary to exercise care in order to eliminate certain of the 
inertia coupling terms in the final flutter equations.) With this in mind, a 
procedure for analytically determining the A; functions from the transfer func~ 
tions discussed in Section 2 of Chapter IX would be as follows: 

Equation (IX-l,) is 


(X-21) Mm," TA + EX, 

Equation (IX-9) states that, in general, 
(x22) lw) = [-1,0°+G + Glo) 

Equation (IX~6) is written here for only the 4 degree of freedom, 
(X-23) X= GA 


In accordance with the requirements for A; (i.e., M,2O and no surface 


inertia effects) the above three equations, when combined, yield the following 


expression: 
ia: en a 


where, the A function required to obtain the A, function is not dependent 
upon uv /bu) « 
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(e) WO DEGREE OF FRELDOM FLUTTER WITH FEEDBACK 
The two degree case considered in Section 3(b) will be used as an example. 
The additional complication due to feedback is that the bending degree of free- 


dom is allowed to include some control surface rotation 4, , such that 


(Note: 4 is a constant along the span of the control surface.) 


The total angle, @ , due to these two degrees of freedom becomes 
(X=26) PP, i. R, h, +3, 


Now, in lieu of proceeding through the mathematical derivation of the 
flutter equations here, they will be obtained by inductive means. This proce- 
dure will often allow equations (or flutter determinant elements) for systems 
with more degrees of freedom to be obtained more readily. 

Referring to (X-10), (X-11), (X-12) and (X-13) and writing them in deter- 


minant form as follows: 


Work Done by: / A Degree of Freedom 


Work Done in: 


[>| 


(X-27) h 
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Degree of fA 
Freedom 


a Open ee age eee ee 


It can be seen that each element represents the total work done in a 


partic»lar degree of freedom by a particular degree of freedom (e.g., C repre- 


sents total work done in the 4 degree of freedom due to forces arising as a 
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result of @ motion.) 
Now, I represents the work done by 
1. Inertia torque 
2. Aerodynamic torque 
3. Actuator restraining torque (i.e., strain energy) 


By separating (1) and (2) from (3), J can be written as, 


A 


(X-28) I=I-§ 

where from (X-13), 

(X-29) I=iI+#Q 

and 

(X-30) S = Kal, +7g) 
ae "74 


Hence, J “represents the work done by the inertia and aerodynamic torques 
only. 

Now, inasmuch as the "adjusted" bending mode due to feedback consists of 

1. The bending mode as represented in (X~27) 
plus 

2. A motion equal to Ah. 
and, inasmuch as this @ motion does not include any actuator restraining 
torque (due to required condition that this torque be zero), it follows that 
the work done in the "adjusted" bending degree of freedom by the "adjusted" 
bending degree of freedom would include a term, R. I Z This term arises from 
work done by tho (3 component of Iray in the f portions of the Pad’ dogree 


of freedom. Similarly, the work done by the (@ component of lady in the A, 
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portion of Aas » would be represented by R, C (see X~27). Also, wrk done by 
the 4, component of / i in the @ portion of Pacty is expressed by 4G, 
etc. 

Extension of this reasoning allows the determinant elements for the case 


represented by (X-27) with the addition of feedback to be written as 


A+R (C+G)eRI | €+RI° 


For 4,=O , the above determinant becomes identical with the "no feedback" 
case represented by (X-27). 

It is worthyof note that the determinant elements for the case with feed- 
back can be written by combining the elements for the "no feedback" case as 
shown in (X-27) and (X-31). This is advantageous in that if the determinant 
elements (for a given value of y/, ba) ) have been calculated for a "no feedback 
case", the additional work necessitated by the inclusion of the "induced" 7 , 
due to feedback, is not prohibitive (for a specified value of #, ). 

Consideration will be given to three special cases of a system represented 
by (X=-31), each of which requires a different procedure for solution. These 
are as follows: 

1. A, does not vary with frequency, @) . 

2. The general case where f, is a function of frequency, @). 

3. The case where flutter of the two degree of freedom system represented 

by (X-27) exists, and a "beneficial feedback," sufficient to stabilize 


this flutter mode, is desired, 
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For the case where A, (a complex number, in general) does not vary with 
frequency, a solution can be obtained by following the procedure of Section 
3(b). The essential difference is that the determinant elements on the right 
hand side of (X-1,) would be replaced by those shown in (X-31). For a known 
value of #,, the additional terms can be readily evaluated. 

This case can be thought of as a control surface which is geared to parent 
surface bending by a linkage with a (complex) gear ratio = A, . The chart dis- 
cussed in Section 3(b) would be applicable here. 

The general case in which the magnitude and/or phase angle of the vector 
A, may vary with the frequency of oscillation requires that an additional con~ 
dition be satisfied; namely, that compatibility be achieved between the fre- 
quency associated with the value of &, used in an analysis and the frequency 
resulting from said analysis, 

Three slightly different procedures for solution are presented. They are 
all similar in that they require more computations than do any of the previously 
discussed systems. Each of these proposed methods requires that different com 
binations of three of the following five characteristics of the structure be 
known prior to solving the flutter determinant for the remaining two character- 
istics, These five characteristics are: 

1. The real part of the available complex stiffness (i.e., Kp ) asa 

function of @. 
2. The imaginary part of the available complex stiffness (i.e., K,) as 
a function of @. 

3- Both the real and imaginary parts (or the magnitude and phase angle) 
of A, as a function of & , 

4. The value of the bending frequency, a, . 


de The value of the available structural damping coefficient, J, e 
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Method 1. 


Here, it is required that characteristics (3), (4), and (5) be known. 


The object is to solve (X-31) for the value of S (i.e., K,/ w)* and Gy ) 


required to sustain constant amplitude flutter. The procedure is as follows: 


l. 


26 


30 


4. 


De 


Youth, 


Select a value of u/bw . Compute aerodynamic terms. 

Select a value of a). (Note: This fixes A and %, .) 

The only unknown left in (X-31) is §. Solving for S yields the 
value of X,/s)" and the value of g, required. 

For the same value of u/bs), repeat Steps 2 and 3 for several values 
Of @) « 

Repeat Steps 1 thru 4 for several values of u/bw. 


Plot the results as curves of KX, / w)” v8.4) und Jy V8- @ for each 


value of u/ba) (see Figure X-8), 
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7 When characteristics (1) and (2) are obtained, either by analysis, 
measurement, estimation, etc., plot them as showm in Figure X-§ 
(shown as dashed curves). Intersection of the curve of available 
Kk, Je) with the curves required k,/w)* determine points at which com 
patibility exists. Associated values of GY are obtained from corre- 
sponding curves of J, required. 
8, From the resultant values of a. V), y required, and A, /K, available 
a plot of g, (see (X-9)) vse. V can be made in the usual manner. 
One advantage of this method over the following two methods is that the effect 
on the flutter characteristics due to arbitrary variations of A, (available) 
and/or A, (available) can be found with a small amount of additional work. 
Method 2, 
Here it is required that characteristics (1), (2), and (3) be know. The 
object is to solve for the values of a), and g required for flutter to occur. 
For this method, the fictitious damping, Q , is not allowed to exist. Hence, 


% = (K, /Ka) or 


K kK, 
ae —_£ eee 
2 SA(i 7 Z} 


The procedure follows: 

1. Select a value of v/ha) . Compute aerodynamic terms. 

2. Select a value of &. (Notes This fixes 5 and R, )e 

3. The only unknowns left in (X-31) are G@), and g, . Find values of 4), 
and g, required for the determinant to vanish. 

1, For the same value of u/ba) , repeat Steps 2 and 3 for several values 
of w. 

5. Repeat Steps 1 thru 4 for several values of u/ba@. 

6. Plot the results as curves of 4), vs. @ and g, vs. @ for each value 
of u/bhe). (Seo Figure X-9), 
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7. For a fixed value of @, (shown by single dashed line in figure shown 
below), the resulting intersections with the curves of required @, 
yield compatible values of @). Associated values of g are obtained 
from corresponding curves of 9g, required. 

8. From the resultant values of a) (\. V) and g, required, a plot of g V8. 
V can be made. Intersection of this curve with the available value 


of Y, indicates the flutter speed in the usual manner. 





Figure X~9, 
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One advantage of Method 2 over Method 1 is that the variation of flutter 
speed with bending frequency (for the same bending mode shape) can be found 
without too much effort. In practice, however, this is not likely to prove 
very useful for this type of problem unless a preliminary study is being made 
wherein /, , A, , and A, functions are estimated. 

Method 3. 

Here, it is required that characteristic (3) be known. In addition, the 

ratio of characteristic (1) to characteristic (4) is necessary (i.e., K, /a\. 


as a function of @), Also G% and Y are each represented by 9, « Let 


(X-32) 4 = Z = Mm 
(X-33) 5 = (Bs \-Sa) (7 + 79,) 


The object is to solve for the values of a), /a) and g, required for flutter 
to occur. (This is similar to the Material Center Method of Air Corps Technical 
Report No. 4798.) The procedure is as follows: 

1. Select a value of u/ha). Compute aerodynamic terms, 

2. Select a value of G. (Note: this fixes XK, /w), and ) | 

3. Now, the terms 4 and S (see (X-10), (X-31), (X~32), and (X~33)) each 

contain the unknown vector, ( @),, /u))'( /+79,) . Solve for the values 
of a)/a) andg,. 

4, For the same value of v/ba), repeat Steps 2 and 3 for several values 

of &. 

5. Repeat Steps 1 thru 4 for several values of v/ba) . 

6. Plot the results as curves of w/w, vs. @) and J, VB- @) for each 


value of y/bw. (See Figure X~10), 
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Figure X--10,. 


7. A straight line drawn from the origin of co-ordinates at an angle, Z, 
from the vertical (see above figure) represents a given value of a), 
(note: @} = TAN QD )e The intersections of this line (shown as a dashed 
line in above figure) with the curves obtained from Step 6 represent 
solutions for this value of a). At these intersections, values of @ 
(and, therefore V ) are obtuined. Associated values of g, are 


obtained from the corresponding curves of g. VS. @) 
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8. Plot V vs. g, « The intersection of this curve with g, equal to 
the minimum available value of X,/A, or g, equal to the available 
value of g, , whichever is less, should yield a conservative flutter 
speed. It should be noted that if the available value of K,/K, , for 
all practical frequencies, is not approximately equal to the available 
structural damping, 9, , then either of the first two methods may be 
more satisfactory insofar as obtaining results that are not unduly 
conservative. 

One advantage of this method over the first two methods is that, for a 
fixed character of the:X, /0), vs. G@) curve, the effect of scaling both X, 
and a), (by the same factors), up or down, can be found in a few minutes. 

In Section 1 of Chapter IX, mention is made of a flutter system stability 
augmenter. This augmenter would be "designed into" the airplane s0 as to sta-~ 
bilize an otherwise unstable flutter mode by providing, in effect, a source of 
"beneficial feedback". The design of a device which would be capable of pro~ 
viding a given amount of feedback that is required is another problem. 

One possible method of determining the required amount of "beneficial 
feedback" will be discussed briefly. 

Suppose a two degree of freedom system (represented by (X-27)) exhibits 
divergent flutter for a certain airspeed range and at associated flutter fre- 
quencies. The required value of feedback could be found as follows: 

1. Allow the existence of some feedback, 4, . 

2. Equation (X-27) would then appear as (X~-31). 

3. From the divergent speed range and associated flutter frequencies, 

select a value of V and the corresponding value of ag). This fixes 
the value of v/ba) . 
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4. let g, and g, be sufficicntly negative. (At most, less than the 
available values), 

5. The only unknown then left in (X-31) is A. Solve for A, . This is 
the amount of feedback which, for constant amplitude flutter to occur, 
now requires more negative amounts of 9, and g than are available. 
Therefore, the available dumping would cause motion to decay. 

6. Using this amount of feedback, establish that this mode, and all other 
possible flutter modes, are sufficiently stable by using one of the 
methods associated with Figures X-8, X-9, or X-10. It should be noted 
that determination of #, vs. @) will be necessary before this step 
can be satisfactorily accomplished. This requires a knowledge of the 
"beneficial feedback" path from sensing element to control surface, 

If this "beneficial feedback" has an adverse influence on a different, and 
otherwise stable, flutter mode, it would be necessary to relocate the "built 
in" source of this feedback (e.g., a sensor) in order to obtain an optimum con - 
dition. 

Experimental verification of this desired feedback might later be realized 
on the ground by measuring #, . 

(f) THREE DEGREE OF FREEDOM FLUTTER WITH FEEDBACK 

For the three degree of freedom flutter case wherein each of the two struc- 
tural degrees of freedom has a certain amount of feedback associated with it, 
the equations for the determinant elements will be obtained by an extension of 
the method employed in Section 3(e). 

The three degree case will consist of: 

l. Any degree of frecdom, P .« 

26 Any other degree of freedom, GF - 


3. A control surface rotational degree of freedom, @ 
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Now. for this system without feedback, the flutter determinant would have 
the following form: 


| 1 
P | G | ° 
i | 
aed 
(X~34) P A ,e , ¢ 
a eee eee ee eee eae 
ty 
g D es | F x O 
sae aed alae oie = 
@ [Ss l4 ‘i 
| 


where, the equations for each of the determinant elements, A thru // would be 
obtained in conventional fashion. J -S remains unchanged. 

The existence of feedback due to > and JY results in control surface 
rotations A and A, » such that 


(X~35) fo ™ Rop 

(X=36) A = kg 

Then the total rotational angle, Pr » becomes 

(X~37) = mpthgtfF& 

The "adjusted" degrees of freedom can be considered as 


Pray ~ P* Pp 
Fahy = %* Py 


(X-38) 
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Hence, each of the "adjusted" degrees of freedom is composed of motion 
represented by the original mode shapes (i.e., and g ) and some motion. 
Now, by considering each "adjusted" degree of freedom to be composed of two 
parts (see (X-38)), the determinant elements for the feedback case can be written 
by properly combining the elements for the "no feedback" system given in (X-34). 

As an example, the correction terms to the 8 element of (X-34,) will be 
obtained. The desired determinant element represents work done in the 2 ely 
degree of freedom due to forces and moments arising from motion in the 7 act 
degree of freedom. This can be divided into four parts as follows (see (X-38): 

1. Work done in ~ due to g motion. 

2. Work done in p due tof motion. 

3. Work done in Pr due to g motion. 

4, Work done in & due to Ky motion. 

From (X-34), it can be seen that 

Work done in ~ due to & motion = C (p8) 

Work done in @ due to g motion = H(g¢/Z) 

Work done in /Z due to @ motion = ([7SYA 4) 

Therefore, inasmuch as 2 , /7y, » and |» motions are identical except for 


amplitude, the above four portions of the "adjusted" element become 
B(pg) = 8(pg) = B(pg) 
C (74) = Clpkg) = CR, (2g) 
H(¢ho) = HG Rp) = HRp(pg) 
I(bp f= I"(Rpp Rg) = I’RpRe (79) 


* Notes S does not appear in this element. 
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and the couplete detemminant elenenrt is 
; ? 

[B+ R,C + Rp * RR I | (pg) 

Similar corrections to cach of the other affected elements of (X~34) result 
in the following flutter determinant for the feedback case, 


A+ Rp(C+G)+ Rel” (B+ R,C+ Ry H*RRI | C+ Rp” 
| 


A 


i 
(x-39) D+ RyF + R.G+ BRINE + RF +H) + RET | FRI” |= 0 
—— — ee YH —» eee ass a aaw aase ame mei as —e a= a= ame aes ae oo 
G*RoI° | H+eRI | 275 


The same reasoning can be applied to cases with more than three degrees 
of freedom. 

It should be noted that for a case in which ao system (possibly not a true 
servo) exhibits irreversible characteristics (i.e., complex stiffness approaches 
an infinite value) and yet is, by some device, able to sense motion elsewhere 
in the structure and convert it into a physical motion of the control surface, 

Equation (X-39) would include only the four elements corresponding to the 
p and g degrees of freedom. However, the influence of the induced ( motion 
in each degree of freedom is introduced by the correction terms which are added 
to the A ,8,C,0D, and — elements. 

The methods presented in Section 3(e) for solving the two degree of free- 
dom case with feedback are applicable to the case represented by (X-39) with 
the following minor changes (note that both A, and AQ as functions of a) 
must be knowm). 

Method 1. 

Hero it is required to know the natural frequencies and coefficients of 

structural dumping, for cach of the » and g dogreos of ircedom (1.0., @),,» 


a > Ty? and Dy must be fisted). 


a h3 
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Method 2. | 
Here, let Gn = Ig = Y), i 
Fix, @,/W, 
’ 
Solve for w, (or a))and 9, « i 
Method 3. ¢ 
Here, let In > Ig ~ Js * Ip p, 
Fix, ap/al as well as X,/ wy ( 
bey 
Solve for @p/4) and g, . 
6: 
The procedure given in Section 3(e) for solving for a "beneficial feedback" L 
can be applied to the three degree of freedom case. Here, either an Ap or Ay 
f 
(not both) would be obtained in order to stabilize an otherwise unstable flutter LL 
mode. ra 
bes 
Gy 
= 
= 
es 
bes 
- 
= 
= 
i 
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